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| Untersuchung des Wellenwiderstandes von Schiffen 
auf flachemWasser bei gleichförmig scherender Grundströmung” 
Von F. Kolberg 


Die Theorie des Wellenwiderstandes von Schiffen auf flachem Wasser bei gleichförmig scherender Grund- 
strömung wird entwickelt unter der Voraussetzung, daß sowohl die zu Grunde gelegten Differentialgleichungen 
als auch die Randbedingungen an der freien Wasseroberfläche wie an der Schiffsoberfläche linearisiert 
werden. Explizite Formeln zur Berechnung des Wellenwiderstandes und des Wellenfeldes weit hinter dem 
Schiff werden angegeben. Die durch die gleichförmig scherende Grundströmung hervorgerufenen Änderungen 
in Wellenwiderstand und Wellenbild gegenüber den entsprechenden Größen in stehendem Wasser werden er- 
örtert. Es ergibt sich, daß im unterkritischen Bereich die maximale Wellenlänge in Schiffslängsrichtung, im 
überkritischen Bereich der Öffnungswinkel des Wellenfeldes in erster Näherung die für die Änderung von 
Wellenwiderstand und Wellenfeld durch die scherende Grundströmung maßgebenden Größen sind. 


The theory of the wave resistance of ships in shallow water with uniformly shearing flow is developed on 
the condition that the basic differential equations as well as the boundary conditions both at the free water surface 
and at the ship’s surface are linearized. Explicite formulae are obtained for computing the wave resistance 
and the characteristics of the wave surface far behind the ship. The modifications in the wave resistance 
and the wave surface due to the uniformly shearing flow as compared with stationary water are discussed. 
In the first approximation the quantities determining these modifications are found to be: (a) in the under- 
eritical region the maximum wave length in the direction of the ship, and (b) in the region above the critical 
velocity the angular aperture of the wave field. 


Pourvu que les equations differentielles prises pour base ainsi que les conditions marginales & la super- 
ficie libre de l’eau aussi bien qu’a la surface du bateau soient linearisees, on developpe la theorie de la resi- 
stance des ondes de bateaux sur l’eau peu profonde avec courant fondamental embardant regulierement. Des 
formules explicites pour la calculation de la resistance des ondes et du champ des ondes dans une grande 
distance derriere le bateau sont indiquees. On discute les changements dans la resistance des ondes et dans 
le champ des ondes provoques par le courant fondamental embardant regulierement compares aux donnees 
correspondantes dans l’eau stagnante. Il en resulte que des donnees determinantes pour le changement de la 
resistance des ondes et du champ des ondes cause par le courant fondamental et embardant sont dans le 
domaine sous-critique la longueur maximale des ondes dans la direction longitudinale du bateau et dans le 
domaine super-critique l’angle d’ouverture du champ des ondes de la premiere approximation. 


PasBuBaetca TeopUA BOJIHOBOTO COHPOTUBJIEHNSI CYIOB HA MEJIKUX BOMOEMAX C PABHOMep- 
HbIM HONepe4HbIM OCHOBHBIM TeYeHHEeM B IIPEeAIIOJIOFKEHHH, YTO IIONIOFKEHHBIE B OCHOBY ANdhe- 
PeHIMaIbHbIe YpaBHeHuA A Tak;ke KpaeBble YCAOBHA HA CBOÖONHOH IIOBePXHOCTH BONbI H HA 
‚HOBEPXHOCTU CyAHa IuHeapn3aupyIoTcH. JlawTcH ABHBIe POPMyJIbI ANA BRIYMCHEHHA BOJIHOBOTO 
COIPOTUBJIEHNUA U BOSIHOBOTO II0OAA HA Ö0JIbBIIOM PAccTOAHHN OT KOpMbI cyaHa. PaccMaTpnH- 
BAI0OTCA U3MeHEeHMA BOAHOBOTO CONPOTUBIEHNUA MH BUNA BOJIH, IIO CPABHEHHIO C COOTBETCBYIO- 
INMMM BeJIMYUHAMU B CTOAUeil Bole, Bbl3bIBAeMbIe IIONePeYHBIM OCHOBHEIM TeyeHHeM. OKA3bI- 
BAeTCH, YTO B HEKPHTNYECKOÜ O0NACTU MAKCHMAJIbHAA NIMHA BOJIHbI B HANPABIIEHHH IPONONB- 
HOüÜ oOcH Cy1Ha, a B CBEPXKPHTNYeCKON O0JTACTH YTON PacTBOPpa BOJIHOBOTO HOJIA B IIePBOM IIPH- 
ÖJIM;KeHUN ABIIAIOTCH BEeJIHYUHAMN, OTNPEeNEeJIAIOINMMM M3MeHEeHNE BOJIHOBOTO CONPOTUBNAECHUA 
uU BOJIHOBOTO IOJIA, BbI3bIBAeMOE TIOMEPEYHLIM OCHOBHEIM TeyeHHeM. 


. 1. Einleitung 

Die Theorie des Wellenwiderstandes von Schiffen wurde von Michell, Havelock, Hog- 
ner, Wigley, Weinblum, Sretensky u.a. zu einem Werkzeug entwickelt, das sowohl in 
qualitativer wie auch in quantitativer Hinsicht in der Lage ist, ganz entscheidende Aussagen zu 
machen über das Verhalten von Schiffen bei gleichförmiger wie auch bei beschleunigter Bewegung 
in ruhiger See und im Seegang. Eine Zusammenfassung dieser Theorie, soweit sie die Bewegung 
in ruhiger See betrifft, gibt Lunde [5]. Bezeichnend für alle diese Untersuchungen ist, daß die 
betrachteten Strömungen als Potentialströmungen angenommen werden. Für unsere Betrach- 
tung, bei der die Grundströmung als wirbelbehaftet angesehen wird, ist es nun nach einem 
Satz von G. I. Taylor [1] nicht möglich, die zu untersuchende dreidimensionale, durch den 
Schiffskörper hervorgerufene Störbewegung als Potentialströmung zu behandeln. Richtung- 
weisend für unsere Darstellung ist die Arbeit von Kärmän-Tsien [2], welche die Prandtlsche 
Tragflügeltheorie auf den Fall einer ungleichförmigen Anströmung verallgemeinert, wobei die 
Linearisierung der auftretenden Differentialgleichungen eine einschneidende Voraussetzung ist. 
Da unsere Theorie sich in wesentlichen Punkten sowohl von der Kärmän-Tsienschen wie von 
der Potentialtheorie unterscheidet, sind die zugrunde liegenden Voraussetzungen neu zu for- 
mulieren. 


*) Mitteilung der Versuchsanstalt für Binnenschiffbau in Duisburg. 
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3, Differentialgleichung des Problems und Darstellung der Lösung 


Betrachtet werde die gleichförmige Bewegung eines Schiffes an der Wasseroberfläche mit 
der Geschwindigkeit U,. Wir beziehen uns auf ein schiffsfestes kartesisches Koordinatensystem, 
die x-Achse weise entgegen der Fahrtrichtung, die y-Achse sei parallel zur Wasseroberfläche und 
senkrecht zur x-Achse und schließlich weise die positive z-Achse nach oben und sei senkrecht zur 
Wasseroberfläche. Sind dann p’ der Druck, o die Dichte und u’, v’, w’ die Geschwindigkeits- 
komponenten in x, y, z-Richtung, so lauten die Bewegungsgleichungen für unsere stationäre 
Bewegung in einer reibungsfreien inkompressiblen Flüssigkeit En 


ee 
u’ rt En a ER 
sen tu --—, I 
Hinzu kommt die Kontinuitätsbedingung 
rt. N N N WE u e 


Ist U,(z) die Geschwindigkeit der Grundströmung und setzen wir 
U) DT u EFT 
so haben wir folgende Darstellung für die Geschwindigkeitskomponenten und den Druck 
vuU+u,r vers twWen, pe BER r2 7 Dr 


u, v, w, p sind hier die durch das Vorhandensein des Schiffskörpers hervorgerufenen Störgrößen, 
die für einen hinreichend schlanken Schiffskörper als von erster Ordnung kleine Größen ange- 
nommen werden, und U ist die allein von z abhängige Geschwindigkeitskomponente der Haupt- 
strömung. Wird (2.4) in (2.1) und (2.2) eingeführt, so erhalten wir das folgende System (2.5) von 
vier simultanen Differentialgleichungen für die vier gesuchten Funktionen u, v, w, p: 


wrote 
Urda +) +u =. 

(U-+u) = + 02 + = = 3. ae. 
ae En 


Wie in der Störungstheorie üblich, nehmen wir an, daß für die überlagerte Strömung die in 
den Störungsgeschwindigkeiten quadratischen Glieder vernachlässigt werden können. Bei Be- 
rücksichtigung dieser Vernachlässigungen haben wir dann zur Bestimmung der u, v, w, p das 
Differentialgleichungssystem: j per 


ou dU 1 op 
DB: E = — 
x yü dz o 0x" 
ov 1 op 
U zn ee 
dr Day (2.6) 
U eR 1 De. op 
OX o 0 
und 
ou ov ow 
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Elimination des Druckes p aus den Gleichungen (2.6) durch kreuzweise Differentiation ergibt: 


0 |ow ov dU ov 
ie ee a u 
x 2 02 dz 0x “ 
oe ou ow dU-9v.. AU | 
Je Se RT ee (2.8). 
0 | ou dU öw 
er ee y ” | 


Eliminiert man in der zweiten Gleichung des Systems (2.8) die Funktion u mittels der Kon- 
tinuitätsbedingung (2.7), so erhält man | 


uw Aw O°Vv ger 00 edel 
ae ee dz &y da“ 


Aus der ersten Gleichung von (2.8) folgt nun durch Integration nach x zwischen den Grenzen 
— oo undx, wobei zu beachten ist, daß für > — ® u,v, w sowie deren Ableitungen ver- 


schwinden: 
ö 
en 0 . (2.10). 


Differenziert man diese Gleichung nach y, so erhält man 


U u N, om, dU 0 
Pr a 


Führt man dies noch in (2.9) ein, so erhält man als partielle Differentialgleichung für die Ge- 
schwindigkeitskomponente w: 
d=U 
Ga dz? 


— = PABBR 
a a a eh 


Nehmen wir nun an, daß w aus der partiellen Differentialgleichung (2.11) bestimmt ist, so ergibt 
sich v durch Integration der Differentialgleichung (2.10), welche auch in der Form 


du 
ov dz ow 

ee 1 
en 212 


geschrieben und somit als gewöhnliche Differentialgleichung erster Ordnung für v aufgefaßt 
werden kann. Nachdem w und v bestimmt sind, erhält man dann u aus der 'Kontinuitäts- 
bedingung (2.7) zu 


x 


/v cw 
Be i _ 2 A s e 5 5 A r . . . . . 213 . 
. (2 £ a) Fr. 


Bemerkenswert ist nun, daß die Differentialgleichung (2.11) für die Geschwindigkeitskomponente w 
sich im Falle eines linearen Geschwindigkeitsprofils der Grundströmung auf die Potential- 
gleichung reduziert, wodurch dieser Fall, der, wie sich noch zeigen wird, die durch die Strömung 
bedingte Veränderung im Wellenbild und Wellenwiderstand prinzipiell richtig trifft, sich durch 
eine leichte Zugänglichkeit auszeichnet. Im folgenden wollen wir uns deshalb ausschließlich mit 
diesem Fall des linearen Geschwindigkeitsprofils der Grundströmung beschäftigen. Hier muß 
also w der Potentialgleichung 

uw. uw. em 


a. 0. oe ee 
oa Ip 02° 2 en 


Ale 


genügen. Die Bestimmung der Geschwindigkeitskomponenten v und u geht bei bekanntem w 
dann wieder über die Gleichungen (2.12) und (2.13). Geht man in der beschriebenen Weise vor, 
19* 
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so erhält man unter Beachtung der Bedingung, daß u, v, w, sowie deren Ableitungen für — 0 
verschwinden und v, w für 2 + 2 > o, —h<z<( beschränkt sind, schließlich: 


2 du rer 


w de w dw 
uw, y,2)=— f [e-5/ fa® dz ee dx, 
BA. (2.15), 
un (Fe I dz, 


Ss 


w(x, y, 2) 


wo w eine Lösung der Potentialgleichung (2.14), U(z) = v* + bz ist und berücksichtigt wurde, 
daß die u,v, w ein Lösungssystem des Differentialgleichungssystems (2.8) sein sollen, also ins- 
besondere der dritten Gleichung des Systems (2.8) genügen müssen. Der zugehörige Druck 
P(z, y, z) ergibt sich aus der zweiten Gleichung des Systems (2.6) unter Beachtung der zweiten 


Gleichung des Systems (2.15) zu: 


z dU z z 
cw dz cw 
p(&,y,2)=—oU S or A ee (2.16). 


Man verifiziert leicht, daß (2.15) und (2.16) dem System (2.6) und (2.7) genügen. 


Stellt man jetzt noch w(z, y, z) in der Form 


dar, wo 
ed 88 a Zu 


0x2 oy? 022 . 


ist, so erhält man für die gesuchte Lösung der Systeme (2.6) und (2.7) unter der Voraussetzung 
U(z) = v* + bz schließlich: 


U: : 
ur, U, ee. [ee 
au h 
re Zu | 
Een (2.17), 
wir, 4,2) = — 8 = 
dU ;: 


DIT, U, 2) = Dr ud 
P(@&, y, 2) =o0 = 77 = 


wo 
ed AD od 
DT aa ta = 


ist und die Randbedingung am Tankboden 


w, _ m 


er er rn 


berücksichtigt wurde. 


P. 
4 
A 
2 
re 
4 
F 
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3. Die Randbedingungen an der freien Oberfläche 
Die in (2.17) vorkommende Potentialfunktion ist erst durch die Randbedingungen des 
Problems eindeutig bestimmt. Diese Randbedingungen sind nun zu formulieren und zwar inter- 
essiert zunächst die Randbedingung an der freien Oberfläche. 


Betrachtet werde ein auf die Fläche S der Wasseroberfläche wirkendes Drucksystem p(z, y). 
Außerhalb der Fläche S ist der Druck an der Wasseroberfläche z = £(z, y) stets gleich dem kon- 


' stanten Atmosphärendruck p,, den wir ohne Einschränkung der Allgemeinheit als p, = 0 an- 


s 
nehmen können. Bezeichnen wir nach Stokes mit 2 die substantielle Ableitung, so gelten an 


Dt 
der Wasseroberfläche z = {(z, y) die Randbedingungen: 
Der a © 
ur An a 02 eo a) 
und | für z = &(@,y) 
DREEL a Pr 
rn a) ee le 
0 0 : 5 a 
Werden {(z, y), ar &y ebenfalls als von erster Ordnung kleine Größen vorausgesetzt, so erhält 
man durch Linearisierung der Gleichungen (3.1) und (3.2) unter Beachtung von (2.4) und (2,6) 
öp _P 
a et se en 0 (3.3) 
und fürz—=0 
& 
Z 21 9 A Er... (3.4) 


Ersetzt man in (3.3) — durch die mit —e multiplizierte linke Seite von (2.6), so erhalten wir 


nach einigen Umformungen als Randbedingung an der freien Oberfläche 


du 

ou dz 9° 1 op 2 
IE | ER 5). 
en Dr 2) 


Nach Einführung der Potentialfunktion ® mittels der Beziehungen (2.17) geht diese Bedingung 
über in 


25 (&® g\0® u) e=— 
nn Non m. — . . . . 3.6), 
E = | oy? 2 +(e 2. de Bude 0X 5 a) 
= £ U=v*+bz 


wo D(z, y) das auf die Wasseroberfläche wirkende Hognersche Drucksystem ist. Für p@, y) = 0 
geht (3.6) über in die für das Michellsche Schiff geltende Randbedingung an der freien Ober- 
fläche 

b 


20. (20 g \ 00 h e=— 

ut BI 170 ee f = * a R3>2). 

ER +e| äyi e+(e + =) = 0 fürz=0 v (3.7) 
—_H U=v*-+b.z 


4. Die Gleichung für die freie Oberfläche 
Setzen wir den Ausdruck für w aus (3.4) in (3.3) ein und kürzen durch U(0) = v*, so erhalten 
wir 
op 0 0p # 
_— 0:0 — = - N 4.1). 
u ® OX 0X er) 
Diese Gleichung nach x integriert zwischen den Grenzen — oo und x ergibt: 


ia, y)= — r - z ar A LEER 2 EISERBUTT (4.2), 
z=0 
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wobei &(z, y) > 0 für &— — w beachtet wurde. Wird hier die aus der ersten Gleichung (2.6) 
durch Integration nach x zwischen den Grenzen — oo und © folgende Beziehung 


z 


a Be 
Pau wde für z=0 Radler 
0 dz | 
eingeführt, bei deren Herleitung wieder u> 0 und p> 0 fürz=0 und? > — © zu beachten 
ist, so erhalten wir . Br 
1 du Dez, y) 
y)=—-—|U-u+— Iwde+—=l.. 2.242.» . (4.3). 
Diese Gleichung geht bei der Einführung von ® mittels (2.17) nach einigen Umformungen über in 
e b 
* = — 
ee ee a e—y* EISEN 
g | x 0x o v* Ba 
2 U@) =v* +bz 


Hier ist wieder p(x, y) das auf die Wasseroberfläche wirkende Hognersche Drucksystem. Beim 
Michellschen Schiff ist D(x, y) = 0 und wir haben als Gleichung der Wasseroberfläche 


z b 
v* [oo od e=— 
BEIDE SER. 9 ee 
ee, ö 4) 
S ei Ul) =v* +bz 


5. Darstellung der Lösung und Berechnung des Wellenwiderstandes 


a) Michellsches Schiff 


Hier ersetzen wir unser Schiff durch eine in der x-z-Ebene angebrachte Quellsenken- 
belegung. Gegenüber der später zu besprechenden Ersetzung des Schiffes durch ein auf die 
Wasseroberfläche wirkendes Drucksystem hat diese Methode den Vorteil, daß einmal der Zu- 
sammenhang zwischen Schiffsoberfläche und äquivalentem Quellsenkensystem zumindest nähe- 
rungsweise bekannt ist und zum anderen das Schiff als Deplacementfahrzeug durch diese Methode 
besser approximiert wird. 

In Anlehnung an Lunde [5] gehen wir aus vom Potential einer Quellsenkenbelegung 
der x-z-Ebene, deren Dichte o*(&,£) sei. Um die Randbedingung am Tankboden 


oD| 
w =— =( 
z=—h 02 z=—h 
zu erfüllen, addieren wir bereits das Potential der an der Ebene z = — h gespiegelten Quellsenken- 


belegung, so daß wir haben 


= |[@9-l, + „dat, | 


n=a@—HP+P+@—); h=a—P+pP+@+2h+l%) 


Für die weitere Überlegung ist es zweckmäßig, dieses Potential ®,, das also der Potential- 
i : oD, | N 
gleichung wie auch der Randbedingung Fr — (0 genügt, in einer anderen Form darzustellen. 
”" ke—h 
Es gilt nämlich für z—£>0 und z+2h+t£>0: 


6.1). 


eur Rn a 
D, ie jo" (&,0)de a [ao [or [K* (i®@*—z)]{exp [K*£] +exp [— K* (2 h-+L)]}dK*,) 
; 32 sr r 9.2 E} 
wo | 
ar = —H cos Ott Hysin Ok en (9.3) 


gesetzt ist. Dieses Potential erfüllt jetzt noch nicht die Randbedingung an der freien Oberfläche. 


!) Von den hier und im nachfolgenden Text auftretenden komplexwertigen Potentialen haben nur deren 
Realteile physikalische Bedeutung. Deshalb sind bei der Berechnung der hydrodynamischen Größen &(%, Y), 


ee den Formeln der Abschnitte 2 und 4 stets die Realteile der jew eiligen komplexwertigen Potentiale 
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Hierzu addieren wir ein weiteres Potential 


2.— | a0» | rer, x) Co [K* (z + h)]exp [i K*- ö*] dK*,, En 


o* = x: cos O* + y- sin O* 


das ebenfalls der Randbedingung z. = (0 genügt. Die Funktion F(©*, K*) ist dann so 
2=—h 


zu bestimmen, daß 
= 18= j dO* [ Bell o*(&, 2) dE dd exp [— K* (i E cos O* + z)]{exp [K* - 2] 
—n N) s 


+ exp [—- K* (2h +{)]} + F(O*, K*) Co) [K* (z + Dj expili K*:o*]dK* . ..(0.5) 
der Randbedingung (3.7) genügt. Geht man mit (5.5) in (3.7) hinein, so folgt: 


2Or,.KR)= -— Ije® &) exp [-iK*.£&cos O*] 


Co] [K* (h +) {K*-+[K$ —e- sin? O*] sec? O*} . exp [— K*h] dede 
X &of [K* ]{K* — [Kf —e. sin? O*] : sec® O* .Tg [K* Al} 


Ki=-g+e, R=— 


va 


Damit erhalten wir das gesuchte Potential: 


rent ee 


] 
| 
Sehens var [00 [et Co) [K* (k+ D{K* + [K} — € - sin? O*] sec? O*} | 
I 


Co) [K* h]{K*— [KF —e- sin? O*] - se O*-Tg[K*h]} | (5,7). 


x Co] [K* (e + h)] exp [K* (i ©* — h)] dK* 


H=@- HH P+@—95 1 =@—H+P+@+2h+dH 
o* = (C—£) cos O* + y-sin O*. 


Wir bestimmen mit diesem Potential ®’ das Wellenbild weit vor bzw. hinter dem Schiff. Nach 
R (4.5) wird 
(at 


“un 5 [- [mes 


2 
£ [= Walss- ] ca +, | [060 | cos 03740? 
Helm mu ER 


hack; (5.8). 


63 


"Sof [K*(h+d)] el [Kg —e: sin? O*] sec? A | [K*h]} exp[—K*h] 
— [Kf# — e- sin? O*] sec? O*. Tg [K* - h] 


0 


x sin [K* (@ — £) cos O*] cos [K* y sin ©*] - are] 


Für 2&— + o strebt hier der erste Term gegen Null. Das Wellenbild weit vor bzw. weit hinter 
dem Schiff wird also bestimmt durch den zweiten Term in (5.8). Dieser fortan allein interessierende 
Bestandteil hat eine reelle Polstelle für 


K* = [Kt —e-si O*]-sea O*Tg [K*h] . ....... (5.9). 
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Um die Lösungen dieser Gleichung zu diskutieren, geben wir ihr folgende Form: 


K'h Ren (5.10). 
Ki: h-sec?O*. I _ - sin? or) 


0 
Diese Gleichung hat als Funktion K* - h betrachtet, dann und nur dann eine reelle Lösung, wenn 


K*.h- see 9*11— 2. - sin O*\ > 1 
0 Kr 
d * 
u K*h> cos? O* + (eh) - sin? O* 
ist. 
Mit sin? 9* = 1 — cos? O* wird aus der obigen Ungleichung: 


K’:h—e-h>(1—e:.h)- c0o80* 


oder A 
Bee v*2 1 Br Let. 
cos? O* < — — ee Ar 7 We * 
1—e-h 1—ceh Fi(li—e-h Fi = 4 
Setzen wir e* B Vg-h 
g-h 
Y Ei : 5.11 
0% = arc cos umge cos! Fi ar 1 ee (5.11), 


so bedeutet die obige Ungleichung also, daß für F}-yl—e-h>1 nur solche Werte von ©* 
in Frage kommen, die größer als © sind: 


v* 
+ 205 — Te cos er A en neh (5.12), 
Fiyl—e-h 
wogegen für F}-Yl—e-h<1 die Größe ©* keinen zusätzlichen Einschränkungen unter- 
worfen ist. Lassen wir den Term, der mit wachsendem Abstand vom Schiff abklingt, fort, so 
haben wir zu untersuchen 


vr 4 Ar 
ee | o*(E,C) dE dd | cos O* dO* 


„ TSorim* (a + O1R* + IKE —e: sin? O4] sector} {Ar —e Xp [R* - Al}exp [— Ken] 
K* — [Kf— € - sin? O*] - sec? O*. Tg [K* h] 


0 
x sin [K* (© — £&) cos 9*] cos [K* y - sin O*] dK*. 
Für das Auftreten von reellen Polstellen des Integranden war die Gleichung (5.9) maßgebend, 
welche wir vorhin diskutiert haben. Alle anderen Polstellen des Integranden sind rein imaginär, 
wie man sieht, wenn £ 
Kir 
eingeführt wird, womit die Gleichung (5.9) übergeht in 
K* = [K} — e sin? O*].. sec? O* .tg [K* h] 
Kr un 
[K$ — e- sin? O*]. h- sec? O* 
und daraus ist alles weitere ersichtlich. (Siehe auch Lambs [6] Diskussion des ebenen Problems.) 


Beachtet man nun noch, daß man den interessierenden Term von £*(x, y) schreiben kann in der 
Form: 


oder 


— tg [K* A], 


n/2 

RER be 

E*z, y) = Re r — |) o*(E, L) de ar | cos» dQ* 
0 


s 
oo 


"Coj [K* (h + DI{R* + [Kg —e- sin? O*] seo? O*} {K* — e: Tg [K* h]}exp [— K* h] 
2 R*— [Kg —e- sim Orlsea OK 


0 


1 i f 
3; [exp [EK*(@—£) cos O* + ysin O*}] + exp[i K* (© — &) cos O* — y.sin o*} j) dK* (5.13), 
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so läßt eh die Integration bezüglich. K* durchführen, wenn man die in Lunde [5] angegebenen 
Konturen benutzt, bei denen jetzt auch noch die auf der imaginären Achse gelegenen Polstellen 
durch kleine Halbkreise zu umgehen sind. Die zu verwendenden Konturen richten sich wie bei 
Lunde danach, ob (x — £&) cos O* + y - sin O* > 0 oder (e — &) - cos O* + y sin O* < 0 ist. 
Dies liefert dann für ö*(x, y) weit vor bzw. weit hinter dem Schiff die Formel 


Era = LEN E j 0*(£, 6) de dt 


e TC [R*(hLOJK* + [Kt —e- „sin? *] sec? O*} {K*—e-Tg [K*-h]texp[— K*-h] 
1— [K$ —e- sin? O*] h- sec? O* . Se [K* - h] 


x cos [K* (e — £&) cos O*] cos [K*- y . sin O*] cos O* dO* (5.14), 
K* = [KL — e- sin? O*]- se O*. Tg [K*h]; 
0 für Fafl e-hel 
für Ffyl—e-h>1 


* 


fi] Gt 
arc cosl.- yı u) 


wo das obere Vorzeichen zu den Wellen vor dem Schiff, das untere zu den Wellen hinter dem 
Schiff gehört. 


Zu beachten ist, daß bei der Herleitung der obigen Formel 
1. der Beitrag der imaginären Achse, 


2. der Beitrag der auf der imaginären Achse gelegenen Polstellen fortgelassen wurde. [Lunde 
S.15, 16 u. 30.] 


Für diese Beiträge läßt sich aber wie bei Lamb [6] nachweisen, daß sie mit wachsendem 
Abstand vom Schiff gegen Null streben. Ferner ist zu beachten, daß der Ausdruck (5.14) für 
C*(x, y) abgesehen von der eben erwähnten Vernachlässigung nur exakt für y = 0, d.h. auf der 


x-Achse ist. Setzt man O* = arc tg — ‚ wobei wir voraussetzen wollen, daß @ — £ und y 


das gleiche Vorzeichen haben, so würde man | genau genommen den folgenden Ausdruck für 
C*(x, y) erhalten: 


Avr ([' 
aN-F- | J o*(&, 2) de dk 


8 


ra Co [K*(h+Z)]{K* + [K$ —e- sin? O*]- sec? O*} {K* —e-Tg [K* h]}exp[—K* -h] 
x | 1— [K4 —e- sin? O*]-h-se® O*:Gea [K*-h]| 


or 


x [cos K* (x — £) cos O*] cos [K* y - sin ©*] - cos O* doO* 


je (&,&) dE dd 


"Col [K*(h + &)]{K* + [Kt —e- sin? O*] sec? O*} {K*— eTg [K* h]} exp [— K*-h] 
= [ 1— [K* —e- sin? O*] - h- sec? O* . See [K* - h] = 


(5.15), 


oO*r 


x sin [K* (e—£) cos O*] sin [K* y sin O*] cos O* dO*, 
Kr=IR: 2: sin! @r]|see 9*. Te |[K*-.h], 
0 N 


Ru 


= 1 Er 
z are 0 | für Fr.fl1—e-h>1 
h 


wo das obere Vorzeichen zu den Wellen vor dem Schiff, das untere zu den Wellen hinter dem 
Schiff gehört. Daraus ist ersichtlich, daß dieser Ausdruck für große Abstände vom Schiff gegen 
den früheren (5.14) konvergiert. 
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- Der Ausdruck (5.14) zeigt nun, daß, wenn man vom Potential ®’ ausgeht, auch unendlich - 
weit vor dem Schiff Wellen vorhanden wären, was physikalisch unmöglich ist. Man sieht nun 
leicht, daß 


li. ©) dede 


ESEL EI EEE sin? O*] sec? O*} exp [— K*-h]- Co [K*@+ h)] 
. fe  "&of [K* - h]{1 — [Kt —e - sin? O*] h- sec? O* - Sec: [K*- h]} 


x sin [K* (© — &) cos O*] cos [K* ysin O*]d@* . . . (5.16), 
K* = [K* — e: sin? 0*] sec? @*. To [K* h] 
0 für FFfi—e-h<i 
EEE ER rn 
Fr/1—e-h 
eine Potentialfunktion ist, die selbst allen Bedingungen unseres Problems genügt und die darüber 


hinaus die Eigenschaft besitzt, daß die Wellen weit vor dem Schiff aufgehoben werden. Das 
vollständige Potential ist somit 


®-8,+0,+16,- [[eeol, +, te« Bil 


ei a [K$ —e- sin? O*] sec? @*}Coj [K* (z+ h)] 


9 = 
arc COS 


Co) [K*-h KR*— [Kt —e: sin? O*] sec? O* Tg [K* - h]} 


x exp [K* (iO* —h)]dK* +4 N | 0*(&,Daedt ugs 


es [K*(h+O)K*+ [RK —e- sin? O*]sec?O*} exp [—K*-h]Coj [K*(z+h)] 
01 [K* - IT — [KF —e : sin? O*] . h. sec! O*. See (K*. h)} 


x sin [K* (e — £) cos O*] cos [K* y - sin O*] dO*, 
n=(@—P’+YP+(l@—i); HSaR—H+P+@+2h+D), \ 
0* = (0— £&) cos O* + y-sin O* 


wo im letzten Integral K* = [Kf — e sin? O*] - sec? O*. Tg [K*. h] und die untere Inte- 
grationsgrenze 


0 für, Ps VL achel 
= = a 
arc cos Fr 7 = a für FF Y1—e-h>1 


ist. Die zu unserem ER D=G®, +0®,+ ©, gehörige Wellenoberfläche 
wird dann in einiger Entfernung hinter dem Schiff dargestellt durch 


*a,y)= fee &)dEdE 


8 


fe [K*(h+OR* + IK —e- sin? O*]sec? O*) {K*—e-Tg[K*-h]} exp[—K*- h] 
6 1— [K* —e- sin? O*]. Rh. see? O*. See [K*. A] 
or 
x 008 [K* (8 — £) cos O*] cos [K* » y - sin @*] cos O* dO*, ein 


K* = [Kg — e- sin? O*] sec? O*. Tg [K* : A], 
0 für-AR ser el 


RR 1 een 
; I en 5 für FA.Yl—e-h>1 
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‚Hieraus ist die Änderung des Wellenbildes bei Fahrt in une NE was uns wi noch 


beschäftigen wird. 


Zu bestimmen ist nun noch der Widerstand, den wir ee Druckint 
Schiffsoberfläche erhalten nach der Formel BE egration, über “ 


R=2 | N p@,0,0)- “n dEBae. 619, 


Es wird es angenommen, daß die Gleichung der Schiffsoberfläche Sehen ist durch 


m 1 RO FR ae 1.02.20) 
und das Schiff zur Längsschiffsebene symmetrisch ist. Nur in diesem Fall können wir ja auch 
das Schiff durch eine Quellsenkenbelegung auf der Längsschiffsebene ersetzen. 

Ist nun y=y(&,z) die Gleichung der Schiffsoberfläche, so bestimmt sich die Dichte 
0*(x, z) der Quellsenkenbelegung nach der Formel 
U@) RZ 
a an de 


Mit er re und dem Potential (5.17) bestimmt sich jetzt en Druck 
p(&, 0, d auf der Längsschiffsebene nach (2. 17) zu 


gr, 2) — 


| dU % 
ER; 08 dz 
P@,0, d)= nn x [mu Me er a 2 


u 0% 
=0- ae: 1 a) 


ıl 
x |e— = : Sp 
| ve —9 +2 +(€0— 0% Ve — 9°: + PP +(d +2 h+0)%° 


a ei 3 E+2h+tL 2) 
er Ze > 2 er in dEd 
U a VE—S+E rl | En: 


rı/2 
Fe | j 0*(&, dee [ cos O* dO* fe [R* (hHO)]{R* + [K* —e- sin? O*] sec? O*}exp [—K*h] 
Fa ö 


= RZEl 


Cof [K* h]{K* — [K# — e- sin? O*] sec? O*-Tg [K*. h]} 


x I Co [K* (€ +] — — . Sin [K* (© N] sin [K* (€ — &) cos O*] dK* 


Cof [K* h] {1 — [Kt —e sin? O*] h- sec? O* . See: [K* h]} 


+ Alle (&,9 «e] Cof [K* (h +) {K* + [K$ — e - sin? O*] sec? O*}exp [— K*h] 


Die 
x 5 Co) [K* C+Hh-7 Sin [K*(d’ + DIE [K* (€ —E) cos O*]cosO*d0*|. . . (5.23). 
Hierbei ist U = U(l'), o*(&,&) = ar a und im letzten Integral sind O* und K* ver- 


knüpft durch 
K* = [K$ —e- sin? 9*] sec 9*.Ta[K*.h]. 


Die untere Integrationsgrenze des letzten Integrals ist wieder 

0 fir eeh<] 
er Beta 2a für FFVi- e:h>1 
a Fry1—e:h] ur „yl—-ece.h> . 


ei) Der. hier auftretende Integrand wird für &' = & singulär. Diese Gerade ist durch einen Cauchyschen 
Hauptwert auszuschließen, da nur “der Grenzwert von p(x, y, 2) für y— 0 gesucht wird. Entsprechendes gilt 
auch für Formel (5.24). 
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Damit erhält man für den Widerstand, wenn man die Bestandteile fortläßt, die bei zwei- 
facher Integration über S den Beitrag Null ergeben: 


| | | vie) Mas az 
EB de ir ale 
zn — = Ve -+e—D YyEe—-p+e+2h+E” 


I U); las az j= wa [m O* a0 [».0m [K* (© +h)] 
o = 


Cof [K*(h+L){K*+ [Kt —e - sin?O*]sec? O*} exp [—K*-h] er Ber 
X Sf [R* AK IK ee OH Re TE RTER 
a/2 
+4 | | Ui) Dasdt | | 21 ar at | [U(@)K*.Cof[K*(C +] —b Sin [K*@+H)]} 
Cof[K* (h+&)]{K* + [Kt —e- sin?O*]sec2@*}exp[— K*-h] 
Co] [K* h]{1—[K} —e- sin? O*]h- sec? O* Sec [K* . A]} 


cos [K* (€’—£) cos 9*] cos 9* d@* 


(5.24), 
wo wieder beim vorletzten Integral 
K* = [K} —e - sin?-0*] sec? 9%. Tg [K*-h]; 
0 für F}-Yl—e-h<1 
9: = ee! 
arc cos BET ech fürF}-yl—e-h>1 
ist. 
Für ein Elementarschiff, d.h. ein Schiff dessen Oberfläche durch 
gay) 2. (3.25) 


gegeben ist, wie auch für ein zum Hauptspant symmetrisches Schiff vereinfacht sich die Formel 
(5.24) wesentlich. Hier verschwinden nämlich die beiden ersten Glieder, und wir haben: 


r-2([vo; Zur; at 


x | {U €) K* Co] [K* (€ + M]—b- Sin [K*(& + M)]) 


 &oj [IK*(h+L)]exp [—K*h]{K* + IK’ —e: sin? O*]se® O*| , (5.26). 
Cof [K* h]{L — [K} — e sin? O*] h. sec! Ö*. Sec [K* h]) 


x cos [K* (£' — £) cos O*] cos 9* dO*, 
K* = [Kö —e - sin! @*] sec! 9* Fo [K*h], 


0 für Rfizern<i 
9 = 


1 NE 
are cos Ir EB: nr Fx-yl—e-h>1 


Unseren Überlegungen hatten wir zugrunde gelegt, daß die Dichte der Quellsenkenbelegung 
sich aus der Gleichung der Schiffsoberfläche bestimmt nach der Formel: 


U(z) 0y 
2 


Dee) EU U) =v* +bz. 


EN, AN 


an U aus su anna ir u 
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Diese Formel müssen wir nun noch begründen. Hierzu haben wir genau wie Lunde den Grenzwert 
der zur Ebene y—= 0 normalen Geschwindigkeitskomponente v für y—0 in einem le 
der Quellsenkenbelegung zu bestimmen. Nun ist nach (2.17) 


= 
od 
DE U 2) = ut dz, 
Pr] 
also das zu bestimmende v: 
lim »60, 9,9 = — 1m 2 + tim # m | a 
y—0 Umso 
wo jetzt (x, 0, z) ein Pt der ee ist. wie bei Lunde [5] folgt nun 
0) 
im * 
RT al Ze a Be (3.28), 


wo das negative Vorzeichen für die positive, das positive Vorzeichen für die negative Seite der 
x-z-Ebene gilt. Ferner läßt sich mit den in der Theorie der Potentiale von Flächenbelegun- 
gen üblichen Hilfsmitteln zeigen, daß 


im (= PD EU 28, 


ist. Damit wird dann also nach (5.27) 
NL E ZU A cha A De u (5.30). 
Beachten wir nun noch, daß nach der kinematischen Bedingung an der Schiffsoberfläche 


°y 0 _ 
ge Bl ee ee (5.31) 
sein muß für y = y(z, z), was linearisiert 
oy 
v2, 0,2) = U ORDER TEL ErTE (3.32) 
lautet, so erhält man zusammen mit (5.30) die Dichte der Quellsenkenbelegung zu 
ee A Day 
or(x,2) Im a Me (3.33), 
oder 
U) &y 
+ = S=ADt Zr ae See . (9.34). 
a2, 2) RE U(l)=v*+b (9.34) 


b) Hognersches Drucksystem 


Wir wollen hier ein über die Oberfläche verteiltes Drucksystem betrachten, wie es Hog- 
ner [8] für den Fall stehenden Wassers getan hat. Der Zusammenhang zwischen Bodenform 
und äquivalenter Druckverteilung bleibt hierbei noch unbekannt. 


Zur Darstellung des vorgegebenen Druckes ist hier das Fouriersche Integraltheorem für 
eine Funktion zweier Veränderlicher anzuwenden, d.h. man geht aus von der Darstellung des 
Druckes in der Form 


Fame ll j da dß j 1 ep lea dr eu nildidn 68). 


Bei Einführung von Polarkoordinaten gemäß 
RC @* B= K*.sin O* 
da dß = K* dK* dO* 


haben wir für den vorgegebenen Druck die Darstellung 


PD, y) = | Pe n) dEdn je [or [i K* {(e — &) cos ©* + (y—n) sin O*}] K* dK* 
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oder x = x a 
P(&, y) = - | | P(&, n) de dn | dQ* j exp [| K*@*] K* dK* . . . . (6.36), 
IT 
8 —n 0 
wo 
or = (2 —9)cs 9% +y—n)sinO* ....... 0. (5.37) 
ist. 


Mit S wird hier das Gebiet bezeichnet, auf das das Drucksystem wirkt. Für das Potential 
wird man dann auf den Ansatz ' 


D, | |P® n).dE | or | res K*) exp [i K* ©*] Co] [K* @+ h)] dK* (5.38) 


geführt. 
®, genügt der Potentialgleichung sowie der Randbedingung Er N y- 0. Geht man 
z=Hh 
mit (5.36) und 5.38) in die Randbedingung (3.6) hinein, so folgt 
a K* .sec O* f 
AUbrISR- eo v* Cof [K* h]{K* — [K$ — e - sin? O*] sec? O* Tg [K* h]} 
wo 
ch \ . v* 
RENT Yg-h 
gesetzt ist. Somit gilt also 
era + dO* 
Ds N d£dn [sec do 
2 exp [i K* ©*] Eof [K* (2 + h)]K* 3 
Coj [K*- A] {K* — [K, — € - sin? O*] sec? O* Tg [K* - re j 


0 
oder 
n/2 


1 Je 
D, | |P6 n) d& in [ee O* dO* 
0 


| sin [K* (® — 8) cos ©*] cos [K* (y—n) sin O*] Co] [K* (+ h] 


ren, — * r% s 
Coj [K* - AI {K* — [K,— e - sin? O*] sec O* Tg [K*- hl} K* dK* (5.40). 


0 


Bestimmt man mit diesem ®, das Wellenbild weit vor bzw. hinter dem Drucksystem nach (4. 4) 
[bier ist p(z, y) = 0], so wird: 


n/2 
DEE BET RES NN 
* * 
eg SE | renea | 62 fi * — [K} — e- sin? O*] sec? O* . Tg [K* h] 


x cos [K* (0 — £) cos O*] cos [K* (y —n) sin O*] dK*. 
Für eine reelle Polstelle des Integranden muß 
K* = [K5 — e: sin? O*] sec? O* Tg [K* A] 
sein. Es folgt dann, daß das Wellenprofil in hinreichender Entfernung vom Drucksystem dar- 


gestellt ie, durch 
2 


K* {K* — 2: %g [K* h]} 
‚| Dan |, _IKE 


Sr, y) = 
— & - sin? O*] h. sec? O* . Sec? [K* h] 


7.0. 


S 


Ss or 


x sin [K* (© — £) cos O*] cos [K* (y — n) sin O*] dO* , 


e 
E 
er 
e 
F 
£ 
E 
E 
4 
> 


wi 


A TREE 


F. Kolberg, Wellenwiderstand von Schiffen auf flachem Wasser ... 267 


wo ©*, K* durch die Beziehung 
K* = [Kf —e- sin? @*] sec? O* Tg [K* h] ; 


DE era lt,-eene) 

9 = T 
arc COS I — — ———— ij * i 

Be yı a für Fi Vi ern >> 


verknüpft sind und das obere Vorzeichen für die Wellen vor dem Störkörper, das untere für die 
Wellen hinter dem Störkörper gilt. Mithin lautet das Potential d,, das die Wellen vor dem Stör- 
körper aufhebt: 


er ll P& m) dedn 


RR sec +30 [74 (2 u nl 


* | Co [K*A]IL— [Kt —e sin? O#]h. sc O*. ©e@ [K*h]) 
or ß 
x cos [K* (x — &) cos O*] cos [K* (y—.n) sin O*] dO*, [ oo) 
K*=|R, -esin: OF] seckQ*+,- 2. NER], 
0 für FEVI—e-h<1 
rg arc cos ne für Ft\1—e-h 
erh ur Fy 1—e-:-h>1 


Die zu unserem Geschwindigkeitspotential © = ®, + D, gehörige Wellenoberfläche wird dann 
in einiger Entfernung hinter dem Störkörper dargestellt durch 


= KrIRt—e. Toll) 
p(& m) ds leeres sin? O*] h - sec? O* See? [K* - h] | 


EHlT, y) — 


x sin [K* (@ — £&) cos O*] cos [K*(y— n) sin O*] dO*, 
Kr [Kt —e- sin? O*]se O*%4 [K*h], 


(5.42), 


0 rel — e..heı 
0 = IE Ne 


woraus die Änderung des Wellenbildes bei Fahrt gegen oder mit Strömung abzulesen ist. 
Zu bestimmen ist nun noch der Wellenwiderstand. Er berechnet sich nach der Formel 


R* = |re: m) DSH er oe (5.43), 


VE fd D wir op Er 
ee, (Ma FE Nr er WE: (3.44) 


—h 


zu setzen ist. Da an den Grenzen des Bereichs S der Druck p(z, y) gleich dem Atmosphärendruck 
ist, wird 


ie &,m) 
[re a) has dn’ -| [pe gear - are» \av=o. 
- 2 Sm) 
* 


so daß der letzte Term bei = keinen Beitrag zum Wellenwiderstand liefert. Läßt man nun 


noch die Bestandteile fort, die nach zweifacher Integration über die Druckfläche den Beitrag 
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Null ergeben, so erhält man 


m...) [p@ ma | [pe.n)aar 


a2 


K*cosO*{K*—e-Tg[K*h 
x [1-jer- eantorThnserOr Sec |KFH] z 
Z % IOABE, 
x cos [K* (&' — £&) cos O*] cos [K* (n’ — n) sin O*] dO*, 


K* = [K* — e- sin? O*] sec? 9* Tg [K*h], 


0 für FF —e-h<1 
9 = 1 
AaIt COy ne fü Fr 1 ut h 1 
ee 
oder wegen 
K*— 2. Tg[K*-h] = se 0*2g [K* - h] 
und 
* in2 * 2 * 2 %* 2K*h 
[X —£.-sın [0] ] A « SEC [0] Sec [K h] = Sin [2 K*h] 
mit 
ie U* 1 en ER . S z ? 
P=| DE) dEdn; — | (pen) IR* (cos O* + 7sin O*)]d£dn (5.46), 
V* P —sın 
schließlich 
ee EEE 2(3,) sec 9* K* Sin? [K* h] 
er Kü ven pRrn —aKen dO*, 
O%* 


K* = [K* — ec sin? O*] sec? O* Tg [K*h], BAM 
0 für FF —e-h<1 


1 | ze 
are OR für PT ee 
[5 el für Fefl- sıE> 


9 = 


c) Hognersches Drucksystem. Ebenes Problem 


Hier entfällt die Abhängigkeit von y, die Geschwindigkeitskomponente » ist identisch 
Null, und die Gleichungen (2.17) und (3.6) des früheren Textes reduzieren sich auf: 


oD od 
u=— N D=— -—, 
X 02 | 
AR ı1X9, 48) 
0oD de f» g 
N = u, che: 
pi, ) 0 l | er 14 “ 0x Fi | 
wo 
2) N od _o 
Or 
ist und 
l dp „| D LTE ; 
[& dr n = ” ( % ee) a NT 
in  UrToD | far an Kae RN; 
Adn) ? B 2 Er dz Perle ER rare (5.50), 


—ı 


wo p(x) das auf die Strecke S der Wasseroberfläche wirkende Drucksystem ist. Der Zusammen- 
hang zwischen Bodenform und äquivalenter Druckverteilung bleibt hier noch unbekannt. 
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Wendet man zur Darstellung des vorgegebenen Druckes das Fouriersche Integraltheorem 


für eine Funktion einer Veränderlichen an,.d. h. geht man aus von der Darstellung des Druckes 
in der Form 


Po). | D(&) de j SED IERF 0 Dan ee (5.51), 


so wird man für das Potential ®, auf den Ansatz 
DE) = 2. |P@ ae | 7) exp [[ K* (ec —£)]Co] [K*(z+h)] dK* . . (5.52) 


N R A ö ® 
geführt. ®, genügt der Randbedingung =( sowie der Potentialgleichung. Geht 


= 


i —h 
man mit (5.51) und (5.52) in die Randbedingung (5.49) hinein, so folgt 


DEREN 2 ER A 
FRI. D8 Sf [KFHIIER RE SIR NT) 


gesetzt ist. Somit gilt also 


Er = f Cof [K* (@+h)] 
D, <: 0: „50 a | [K* Rh] {K* _R: . Tg [K* h]} sın [K* (x — &)] dK* (5.53). 


Bestimmt man mit diesem ®, das Wellenbild weit vor bzw. hinter dem Drucksystem nach (5.50) 
[hier ist p(x) = 0], so wird 


Re nern 
*Ge) = Eu u. | PO de | ge RE. &g [R# n]°° [K* @ — | aK*. 
5 0 


Die einzige reelle Polstelle des Integranden liegt bei 
K* = Ki%g [K* h] 
für den Fall, daß Kf-h>1 ist. Es folgt dann, daß das Wellenprofil in hinreichender Ent- 
fernung vom Drucksystem dargestellt wird durch 
Kr Bere Kr Ah]: 
| est A NL * (m — 
OF | POL mens ren REDE, 

wo K* die reelle Lösung der Gleichung 

ei SEE RR Va a ee 2 a Ze A (3.54) 


ist und das obere Vorzeichen für die Wellen vor dem Störkörper, das untere für die Wellen hinter 
dem Störkörper gilt. Mithin lautet das Potential ®,, das die Wellen vor dem Störkörper aufhebt 


Ile Co] [K* (c+ N] 
Fi ee Ken — Kir Re end, (5.55). 


a 
1 


8 
Ken KM. 20TK* eh] 
Die zu unserem Geschwindigkeitspotential D=@, + ®, gehörige Wellenoberfläche wird dann 
in einiger Entfernung hinter dem Störkörper dargestellt durch 


BE Te K hl, en 
= ,.,|P@ I—K*:h&c® [K*h] en: ee #0:58). 


Kr aa lK*-n] 
Zu bestimmen ist nun noch der Wellenwiderstand. Er berechnet sich nach der Formel 


N (5.57), 
| dE’ 


s 
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wo 
de*  v* 


ea 5 od Ps 1 | 
dx g BF 


Pr 77, o v* "de 
—h 
zu setzen ist. Da an den beiden Grenzen des Bereichs $ der Druck (x) gleich dem Atmosphären- 


druck ist, wird 
1 e = 1/2 


% dp ET. 
eo gez), 


8 


En 
so daß der letzte Term bei nn keinen Beitrag zum Wellenwiderstand liefert. Läßt man nun noch 
die Bestandteile fort, die nach zweifacher Integration über die Drucklinie den Beitrag Null er- 


geben, so erhält man 


i> K*.1K*—e-Tg [K*-h 
Re. | Bio ae | PIE) at cos [ne ea N ee ae 


K*=Kt:-Tg [K*:h]. 
Mit 


P= jo de, U Hi V+r = » [P® BEE (5.59) 


ergibt sich schließlich: 
N 
1—K*h-©e@[K*h] 


R* > ULR®2 V *2 
= — 4 7 r. 
29 4 
K’= Kg KH. I]: 
6. Vergleich des Wellenwiderstandes und des Wellenfeldes 


für strömendes und stehendes Wasser 


Zur Untersuchung der Formeln für strömendes Wasser und stehendes Wasser wollen wir die 
sich auf strömendes Wasser beziehenden Größen mit einem *, die sich auf stehendes Wasser 
beziehenden Größen ohne * bezeichnen. Die Geschwindigkeit in stehendem Wasser ist also v, 
die Wellenzahl X usw. 


A) Zunächst betrachten wir den Wellenwiderstand. Diesem lassen sich nun verschiedene 
Formen geben [10], je nachdem, ob wie bisher O* oder K* die Integrationsvariable ist. Der 
Zusammenhang zwischen diesen beiden Größen ist nun gegeben durch 


K* = [Kf —e- sin? O*]se® O*Tg[Kthl, Ki=L+te, 
D 
oder 
r G 9 m — 
K* — & sec? O* + | To [Re Hl 20 (6.1), 


und durch partielle Differentiation folgt hieraus: 


dK* 2{K* —e.Tg[K*h]} Sin [2 K* h]tg @* 


AO Sin [2 K*h]— 2 K*h 2. See (6.2) 
oder 
dK* _4g sec? O*tg O*- Gin? [K*h] 
dO* vr Gn[2K*hl—2K*h (6.3). 
Ferner gelten die Relationen: 
*2 
08 O* — _glv r 
cos Ver IK* h] Tg [K*h], 
| ar Be 
sin O* = 1 nn u en gjv* TalK*h 
| K* —e.%g [Ken 5 8 ] 
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a) Michellsches Schiff 


BANIE beschränken uns hier auf ein Elementarschiff oder ein zum Hauptspant symmetrisches 
Schiff. Dann haben wir also nach Formel (5.26) für den Wellenwiderstand 


n/2 3 
8 4 
R+=28 I VO LEO IR* (R HD] az ag 
on 5 


x (ber {U(E) K* Cof [K* (© + N] —b- Sin [K* (EC + A)]} de’ ae’ 


2 K* .cos O* 
©in [2 K* h|—2K*h 


cos [K* (£&' — £E) cos O*] dO*, (6.5), 
K* = Pa sec? O* 1 | Tg > h] 3 
0 für Pen ee 


= 1 2 
0 ÄnCIcoSee ee 3: I Rn: 
e | für Bde yI, h>1 


in strömendem Wasser, wogegen in stehendem Wasser 


* 


r/2 
8 6 0) 
R = Ft Ie@ Co [IK (h + 9] de ae | ozgssi [K(h+ ON] dr ar 
(cH 8 E 
K? cos © s 
x EnPKHI—2K 7, c08 [& (€ — &) cos 9] dO, (6.6). 
e 0 Fur =’ 
z a rd sg arc cos {m fur Pl 
h 


Führt man in (6.5) bzw. (6.6) anstelle von ©* bzw. © die Variable K* bzw. K als Integrations- 
variable ein, so erhält man in strömendem Wasser ; 


ee I U Cof [K* (h + O]dEat 
T 0& 
Kon 


4 | | “ {U(E) K* Cof [K* (€ + h]I—b- Sin [K* (€ + h)]} de’ u 


(6.7) 
K* . cos O* 
DE ee: - = a x (BE, * * 
* Sin [2 K* HJ {K* — e- Tg [K* - h]} tg ©* 608 [K* (& —£) cos O#]/dK*, 
K* = [Kf* — e- sin? O*] se 9* Tg [K*-h]; Kr, = Rd Ta [Kı. hl; KR >= mn AB: 
und in stehendem Wasser: 
40 Son : .[L[.9 7 
er v— Co [K(h+H]dEdt | | v-,Cof [K (€ + hy] dE’ dd 
TT 0E dE 
Kon =, s 
K.cos © (6.8). 
Fe A ;O]dK, 
er ger ae 3 
KR: seo [K-hH:; K= KB - TolKirib Ki- na 


Für Kt-h<1 bzw. K,:h<1list Kf,=0 bzw. Kyn—=0, d.h. im überkritischen Bereich 
werden bei den Integralen (6.7) und (6.8) die unteren Integrationsgrenzen Null. Dies legt nahe, 
zwecks Vergleichung des Wellenwiderstandes die Formeln (6.7) und (6.8) im unterkritischen 
Bereich zu benutzen, wogegen man im überkritischen Bereich zweckmäßig auf die Formeln (6.5) 
und (6.6) zurückgreift. 


20* 
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Setzen wir jetzt in ähnlicher Weise wie Weinblum [7] voraus, daß die Integranden in 
(6.5) und (6.6) bzw. in (6.7) und (6.8) ungefähr gleich sind, so unterscheiden sich die Integrale 
(6.5) und (6.6) bzw. (6.7) und (6.8) im wesentlichen durch die untere Integrationsgrenze. Bei glei- 
cher Integrationsgrenze werden wir also in strömendem und stehendem Wasser in erster Näherung 
gleiche Wellenwiderstände erhalten. Das bedeutet: Der Wellenwiderstand für eine bestimmte 
Geschwindigkeit v in stehendem Wasser entspricht dem bei strömendem Wasser bei der Ge- 
schwindigkeit v*. Letztere ist aus der Beziehung 


ei (6.9) 
An € BR er a Pe 
im unterkritischen Bereich d.h. für F, < 1 bzw. 
Fr 1 
FM 1 N (6.10) 


im überkritischen Bereich d.h. für F, > 1 zu bestimmen. Für F, = 1 ergibt sich aus den beiden 
Relationen dieselbe Geschwindigkeit v* aus 


SErFano ohre 2 
VE, ze 


b) Hognersches Drucksystem 
Führt man mittels (6.3) X* als Integrationsvariable in (5.47) ein, so erhält man 


es 9 -12 
R* > i U*2 y*r K*? (( v* TS [X * h 
N sr K* — e-Zg[K*-h] 8 | 
Kon 
I 1 | 
a I ee er... (6.11) 
1 y 2 «© * s 
Kran [K#.1]3 [KX* h] er: 
wo Ky, die Wurzel der Gleichung 
7 G 
Kh=Kt:Tg[Kt,-h], Ki — er a (6.12) 
ist. Ferner ist | 
U* 1 ge Ensat = e 
„ =, | DE, Eu [K* (& cos O&* + sin O*)] d& dn. 
Die entsprechende Formel in stehendem Wasser ist 
+ g — 1/2 g 1/2 
R- —_ [wıwrelı sem) —lı Teen! ler 
Irog. RK Tree 
Kon e 
6.13). 
Kon Ro 29 [Kahl], Kom. Be 
U a ar C08..; N | 
” =5 I P&n) ri [K* & cos 9 +7: sinO)]| d&dn 
Die früheren Formeln für den Wellenwiderstand waren 
ei 9 749 u 
: 9 pı za" Sec O* K*. Sin? [K* h] 
2% — U* ı yx2 D* 
70.g \ | ) Sin 2 K*h] — 2 K*h io 
% 
ak — ER: O* ı .|7 73% 
# v*2 * SEC a N [X h] E | a (6.14) 
| 0 für Fer ne 
AR jarc cos | für Fry1 
"rien ul Be 
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in strömendem Wasser bzw. 


rı/2 


0, p2 - . sec O- K?. Sin? [K h] 
fe + va)? dO 


neo: ©inPKhl—2Kh 


& el) 
0 fürn, 1 


I saca c = 
e pa 9° 22, = arc cos Im für, >1 
Fi, 
in stehendem Wasser. 

Setzen wir jetzt wieder voraus, daß die Integranden in (6.11) und (6.13) bzw. in (6.14) und 
(6.15) ungefähr gleich sind, dann unterscheiden sich die Integrale (6.11) und (6.13) bzw. (6.14) 
und (6.15) im wesentlichen durch die untere Integrationsgrenze. Bei gleicher unterer Integra- 
tionsgrenze werden wir also in strömendem und stehendem Wasser in erster Näherung gleiche 
Wellenwiderstände erhalten. Das bedeutet: 

Der Wellenwiderstand für eine bestimmte Geschwindigkeit v in stehendem Wasser ent- 
spricht dem in strömendem Wasser bei der Geschwindigkeit v*. Letztere ist aus der Beziehung 


PR 


un 
F} +e-hF 
im unterkritischen Bereich, d.h. für F, < 1 bzw. 
FR 1 
Fer gr 


im überkritischen Bereich, d.h. für F,> 1 zu bestimmen. Für F,=1 ergibt sich in beiden 
Fällen v* aus der Relation 


Ben 


Fi = ——— zu = ®3t4)o- + 
re 


Der physikalische Sinn der gemachten Annahme über die Gleichheit der Integranden von R* 
und R ist, daß bei Fahrt auf strömendem Wasser mit der Geschwindigkeit v* Wellenhöhe, Länge 
und Form dieselben sind wie auf stehendem Wasser bei der entsprechenden Geschwindigkeit v. 


Die von uns aufgestellten Zuordnungen der Geschwindigkeiten v und v* bedeuten, daß man 
im unterkritischen Bereich gleiche minimale Wellenzahlen in Schiffslängsrichtung verlangt, 
wogegen man im überkritischen Bereich den minimalen Richtungswinkel der Wellen als charak- 
teristische Größe auszeichnet. Da zwischen minimalem Richtungswinkel 0% bzw. ©, und dem 


Öffnungswinkel des Wellenfeldes ö* bzw. ö der Zusammenhang 6* = + & — 0%) bzw. 
Ü= +5 — 6) besteht, heißt das letztere, daß man im überkritischen Bereich den Öffnungs- 


winkel des Wellenfeldes als charakteristische Größe ansieht. 


c) Hognersches Drucksystem. Ebenes Problem 
Hier haben wir für den Wellenwiderstand in strömendem Wasser die Formel: 


g 
a8 K*.K# Tg [K* h] 
* 0, Fo Er ea 2 2 . . 3 

i > = e Bere ker] 1 
E 
r p’*2 
2 wo 
ä K* = K* Tg [K*-H], Kate, 
; Ur +iV*=,, "D(e) ers Be. ; =... (6.12% 


Ss 
wogegen in stehendem Wasser 


| PB K-Ry-%g IK h] 
R-, 04m, Ey “= nl ee  AENERE eic (6.18) 
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mit 


K=K,-2g[Kh]l, R=) 


De = | den ee (6.19) 


u T D(e) de 


ist. 


hat in beiden Fällen denselben Wert. Hieraus entnimmt man: 


Ist R der zu einer gewissen unterkritischen Geschwindigkeit vin stehendem Wasser gehörige 
Wellenwiderstand, so bestimme man die Geschwindigkeit v* aus der Beziehung , 


Fi? 
— —=1+e:-h- Fi KEEA) 
Fi 
Dann ist R* ae: der zu v* in strömendem Wasser gehörige Wellenwiderstand. 


B) Wir IR nun über zur Betrachtung des Wellenbildes weit hinter dem Schiff. Hier 
haben wir dann für die Wellenoberfläche die folgenden Formeln. 


a) Michellsches Schiff 
In strömendem Wasser 


* I K* sec O* Co) [K* (h +9] Sin [K* h] 
re N j ee an 0 En Kth—2Kth 


x cos [K* (x — £) cos * cos [K* ysin O*] dO*, (6.20). 
0 für Ff-Yl—e-h<1 
= 1 
arc cos es air Fi-yl—e:h>1 
oder 
4vk (f  K*. Co] [K* (h +9] 
*(2,y) = — (ET) dEdE 
ODE 1 ET 2 2 
% 6.21 
cos O* Kia a cos O* { en ak 
ig O* cos [K* ( — £&) cos O*] cos [K* y - sin O*] dK 
Kin = K$:-Tg [Kt,- hl 
wo 
K* = [Ki —e- sin? O*] see O*Tg [K*hl; Ki=-,+e; 
p*a 
60 = U0) 
in stehendem Wasser 
a/2 [ a 
uff 3 : K sec O Co) [K (h + &)] Sin [K - h] 
a) = o(&, &) de d£ | 
a: Ein? KhjJ—2Kh 
8 ” 
x cos [X (& — &) cos 9] cos [K - y - sin 9] 40, + MD. 


U A | 
1 | N 
arc cos für A, 1 


F)f 


= 


er ie 


ah u A te v 


a 
> 
& 
7 
a 
ner, 
$ 
F 
= 
F 
Rn 
+ 
F 
gr 
: 
A 
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oder 
ab K-Co [IK (k+9] 
a) 7 [ee 018, &) de de | a 
B . (6.23), 
x Par) cos [K (ec — £) cos 9] cos [K -y- sin 9] dK, 
Kon = Kg: Tg [Kon hl 
wo 
K=-Ksa0.2[Khl, K=-, d&H=v.l. 
b) Hognersches Drucksystem 
In strömendem Wasser 
zı/2 
6*(z, y) 4 P(&,n) d& d ga sec? O* . Sin? [K* h] 
T, — ——— ’ v 
9 TEEN ER De Sin P K*hl—2K*h 
s CH 
x sin [K* (e — £) cos O*] cos [K* (y — n) sin O*] dO* “2... (6.24) 
0 für Ffyl—e-h<1 
* = 1 
arc nen, für Fayl—e-h>1 
oder 
a, y)= ll n)d£ las Try sin [K* (e—£) cos O*] cos [K* (y—n)sin O*]dK*, 
(6.25), 
Kön= Ko: [Kon hl 
wo 
* = [Kf —e- sin? O*] se O*:Ig [K*h]; Kf = te, 
in stehendem Wasser 
rı/2 
A 7, K-see0- Ein2[K-h] 
8 E71 
x sin [K (e — £) cos ©] cos [K (y — n) sin 9] d®, “2... (6.26) 
0 für, <1 
0. 1 2 
arc cos Im für F,>1 
oder 
au Zi] P(&, n) de dn er EA sin [K (x — £&) cos 9] cos [K (y — n) sin O0] dK (27, 
Koh 
Kon = Ko: Tg [Kon h] 
wo 
K= K,se® 9% [Kh]; K=-5. 
c) Hognersches Drucksystem. Ebenes Problem 
In strömendem Wasser 
: s I ee, R* d 6.28 
.,, PO) Ken. Seen | @—sld, RE)" 
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wo 
K*= Kt: Tg [K*h], Kr=te, 


in stehendem Wasser 


ME LLCER an 
=, ) PO Eee 
wo 
K=K,-Tg[Kh], R=.. 


Hier ließe sich nun zunächst ein Vergleich genau so durchführen, wie beim Wellenwiderstand. 
Man würde als Ergebnis erhalten, daß den nach (6.9) bzw. (6.10) zugeordneten Geschwindig- 
keiten v* und v das gleiche Wellenbild entspricht. 


Um hier noch weitergehende Resultate zu erhalten, wollen wir das Wellenbild in strömendem 
und stehendem Wasser bei gleicher Relativgeschwindigkeit (v* = v) untersuchen. Und zwar 
interessieren wir uns bei einer Elementarwelle mit dem Richtungswinkel 0* = © für die Wellen- 
länge und die Wellenamplitude. 


I. Am einfachsten ist die Aussage über die Wellenlänge zu erhalten. Mit v* = v gilt nämlich 


9 
Kent: bzw. Ku= 


Da bei Bergfahrt & > 0, bei Talfahrt e < 0 ist, gilt also für 0* = © 
K+ —e-si0* = K,+e-co®20>K, bei Bergfahrt, 
K+ —e-s0*=K,+e.-co?0<K, bei Talfahrt 
Die Wellenzahlen K* bzw. K bestimmen sich nun aus den Gleichungen 
K*.h K-h 
* BER ce Ban 
(K,+e:co82O)hset 9 TEN AIR K,:h-sec© 


Die Lösungen K* : hbzw. K - h ermitteln sich in bekannter Weise dadurch, daß man die Gerade _ 


1 1 
durch den Nullpunkt mi a Er BER 
en Nullpunkt mit dem Steigungsmaß (K,T 00) bzw KEG mit 
der Zg-Kurve zum Schnitt bringt. Haben für das betrachtete © beide Gleichungen Lösungen, so 
ersieht man hieraus unter Berücksichtigung von (6.30), daß 


K*> K bei Bergfahrt, | 
K*<K bei Talfahrt | ' 


Wegen des Zusammenhanges zwischen Wellenlänge und Wellenzahl 


Ian E30: 


= Tg [Kh] 


Be 2.4: 


an 2n 
Pe K* bzw. =7 


folgt ohne weiteres 

A* <A bei Bergfahrt, | 

4*>1 bei Talfahrt 7 m a en (6.32), 
d. h. bei gleicher Relativgeschwindigkeit gegenüber Wasser sind die Wellenlängen bei Bergfahrt 
kleiner, bei Talfahrt größer als in stehendem Wasser. 

II. Wir gehen jetzt zur Untersuchung der Wellenamplituden über. 
. a) Beim Michellschen Schiff wollen wir die Ausdrücke (6.20) und (6.22) vergleichen für 
ein Elementarschiff, dessen Oberfläche also durch 
y = I) » F2(2) 


gegeben ist. Hierin sind f(x) bzw. f,(z) Polynome in x bzw. z. Beachtet man 


0*&,0 = Gr u) +5-2-1O)], 


506,88), 
(6,0) =v e 6) 


; 
2 
Pr 
i 
4 
3 
3 
E 
. 
4 
E 
4 
4 
3 
3 
L 


Zei dr 
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wo &-f,(&) SO für diein Frage kommenden {, so sieht man, daß für v* = v und O9* = 9 
bereits die Berücksichtigung der Tiefenverteilung des Deplacements bei Bergfahrt eine kleinere, 
bei Talfahrt eine größere Amplitude der Elementarwelle ergeben würde als in stehendem Wasser. 
Wenn wir also in Formel (6.20) o* = o annehmen und hierbei die Änderung der Wellenamplitude 
noch in derselben Richtung erfolgt, so würde dieser Effekt bei Berücksichtigung der Tiefenver- 
teilung des Deplacements nach (6.33) nur verstärkt. Bei den weiteren Überlegungen beschränken 
wir uns jetzt der Einfachheit halber auf den Fall eines prismatischen Elementarschiffes, wo also 
f.(&) = List. Ferner denken wir uns die aus dem Additionstheorem folgende Umformung 


cos [k (e — £) cos 9] cos [k ysin 9] = 5 [cos [k{@ — &) cos# + ysin 9}] 


+ cos [k { (x — &) cos 9 — y sin 9}]} 
in die Integranden von (6.20) und (6.22) eingeführt. Werden jetzt in den Formeln (6.20) und 


(6.22) die Integrationen nach & und £ ausgeführt, so sehen wir, daß sich das Wellensystem in 
beiden Fällen aufbaut aus elementaren Kosinus- und Sinuswellen mit den Richtungswinkeln 


©* bzw. © und — ©* bzw. — 9, 05 < O* < 5 O0 = welche vom Bug oder Heck 


des Schiffes ausgehen. Diese einzelnen elementaren Kosinus- und Sinuswellen haben für v* = v 


und 9* — © bestimmte von K* hbzw. Kh abhängige Amplituden. Wir berücksichtigen jetzt, 
daß die Formeln (6.20) bzw. (6.22) aus den Formeln (6.21) bzw. (6.23) durch Einführung der 
Integrationsveränderlichen ©* bzw. © anstelle von K* bzw. K hervorgehen. Bei dieser Sub- 
stitution werden die Integranden von (6.21) bzw. (6.23) mit einer Funktion multipliziert, die im 
Falle (6.21) mit der auf der rechten Seite von (6.3) stehenden Funktion identisch ist und die im 
Falle (6.23) aus der in (6.3) rechts stehenden Funktion durch Ersetzen der Größen v*, O*, K* 
durch v, ©, K hervorgeht. Diese beiden Funktionen sind nun bei festen v*, O* bzw. festen v, © 
monoton fallende Funktionen von K* hbzw. Kh. Die Untersuchung der Amplituden der von 
Bug und Heck des Schiffes ausgehenden elementaren Kosinus- und Sinuswellen mit dem Rich- 
tungswinkel O* — © läuft damit bei v* = v bis auf Faktoren, welche monoton fallende Funk- 
tionen von K* hbzw. K h sind, hinaus auf den Vergleich von: 


Sin [K* h] — Sin Ix* h (1 — .) Sin [K hl — Sin IK h ( = 3) 


K* h ©in [K* h] = Kh- ©in [K-h] ae 
Wegen (6.31) ist also die Funktion 
Sin x — Sin [a - x] T 


x- Sinz a re yet 


auf Monotonie zu untersuchen. Nun ist 


a x— Ein [a = Sinz- Sin [x-2] +xCojz- nu [ 2] — Sin?2—oax-Co| [a x2]- Sinz Er, 
dx . 


z-GSinz { . Sin? x 


da für jedes x > 0 
Sin x 

Sina x] 
Damit ist also auch die in (6.34) auftretende Funktion eine monoton fallende Funktion, d. h. also 
wegen (6.31), daß bei gleicher Relativgeschwindigkeit gegenüber Wasser die Wellenamplitude der 
Elementarwellen bei Bergfahrt kleiner, bei Talfahrt größer ist als in stehendem Wasser. 

b) Beim Hognerschen Drucksystem wollen wir die Ausdrücke (6.24) und (6.26) vergleichen. 
Der Einfachheit halber nehmen wir an, daß das Drucksystem auf einen rechteckigen zur x-Achse 
symmetrischen Bereich wirkt. Ferner setzen wir 


1+xbtgr< + o&x6otg [ax] RE 


Die,n)= Const, 


voraus, was nach Weinblum [9] den Widerstand auch für eine komplizierte Druckverteilung in 
der Größenordnung richtig wiedergibt. Führen wir jetzt in (6.24) und (6.26) die Integrationen 
nach & und n aus, so sehen wir, daß sich das Wellensystem in beiden Fällen aufbaut aus elemen- 
taren Sinuswellen mit den Richtungswinkeln ©* bzw. © und — O* bzw. — 0, O4 <= O* s z : 


ee “ von denen je eine von den Eckpunkten des Rechteckbereiches ausgeht. Beach- 


ten wir nun noch, daß auch hier die Formeln (6.24) bzw. (6.26) aus den Formeln (6.25) bzw. (6.27) 


; 1; A 
N nr 
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durch Einführung der Integrationsveränderlichen @* bzw. © anstelle von K* bzw. K er 
gehen, und daß die dabei in den Integranden auftretenden Faktoren wie vorhin bei festen vr, 0 
bzw. festen v, © monoton fallende Funktionen von K* h bzw. K'h sind, so haben wir zur Unter- 
suchung der Amplituden der von den Eckpunkten des Rechteckbereiches ausgehenden Elementar 
wellen mit dem Richtungswinkel O* = © bei v* — v zu vergleichen: ; 
1 1 
Kr und K' } 

Hieraus ist wegen (6.31) ersichtlich, daß auch hier bei gleicher Relativgeschwindigkeit gegenüber 
Wasser die Wellenamplitude der Elementarwellen bei Bergfahrt kleiner, bei Talfahrt größer ist 
als in stehendem Wasser. - 

c) Im ebenen Fall haben wir beim Hognerschen Drucksystem die Ausdrücke (6.28) und 
(6.29) zu vergleichen. Für v* = v kommt dies auf den Vergleich von, 


Tg [K* h] zg IX A] 


ag a NETTE 11769 1 ae TTS Ee 


> 2. (6.35) 


hinaus, wo 
Kt. n= Ric hz6 [REF] bzw. Kh=K,htg [Kh] 


ist. Eliminiert man mit Hilfe der beiden letzten Ausdrücke K# h bzw. K,h und nimmt noch 
einige Umformungen vor, so sind die Ausdrücke 
1 1 


ee a, — und —p—: 2. : (636) 


Khl— or 
ee Tee] 
zu vergleichen. Hierzu untersuchen wir die im Nenner vorkommende Funktion 


E 

ofa& — = — 
Eotg Sin?r 
auf Monotonie. Nun ist 


d 2.212 Co) x 
dx Rx IT Ein S Sin? e Sin ) 2 


z.-Cotgr — 1=botgr- ce — Tgr) >0 


für x > 0. Damit ist also die in (6.36) auftretende Funktion eine monoton fallende Funktion von 
K*h bzw.Kh, d.h. also wegen (6.31), daß bei gleicher Relativgeschwindigkeit gegenüber 
Wasser die Wellenamplitude der Elementarwellen bei Bergfahrt kleiner, bei Talfahrt größer 
ist als in stehendem Wasser. 


für jedes x > 0, da 


7. Zusammenfassung 


Bei Linearisierung der zugrunde liegenden Differentialgleichungen und Randbedingungen 
gelingt es, die Auswirkung eines der Tiefe nach veränderlichen Geschwindigkeitsprofils, welches 
das in fließenden Gewässern wirklich vorhandene gut annähert, auf Wellenfeld und Wellenwider- 
stand bei Fahrt eines Schiffes zu erfassen. Für die durch das Schiff hervorgerufene Störungs- 
bewegung werden die Geschwindigkeitskomponenten wie auch der Druck mittels einer harmo- 
nischen Funktion dargestellt. In unserer Theorie wird das Schiff wie auch sonst üblich einerseits 
ersetzt durch eine Quell-Senken-Belegung auf der Mittschiffsebene (Michellsches Schiff), zum 
anderen durch ein auf die Wasseroberfläche wirkendes Drucksystem (Hognersches Druck- 
system). Die für Wellenwiderstand und Wellenfeld angegebenen Formeln gestatten bei vor- 
gegebener Schiffsform die direkte Berechnung dieser Größen. 

Auf Grund einer Hypothese über die Integranden der Wellenwiderstandsintegrale, die in 
ähnlicher Weise von Weinblum [7] benutzt wurde zur Ermittlung des Flachwassereinflusses 
auf den Wellenwiderstand von Schiffen, werden die durch die scherende Grundströmung bedingten 
Unterschiede im Wellenwiderstand gegenüber stehendem Wasser diskutiert. Im unterkritischen 
Bereich F,< 1 erweist sich die maximale Wellenlänge in Schiffslängsrichtung, im überkritischen 
Bereich F,> 1 der Öffnungswinkel des Wellenfeldes als maßgebende Größe für die Änderung 
von Wellenwiderstand und Wellenfeld durch die scherende Grundströmung. Schließlich wird 
noch nachgewiesen, daß, bezogen auf die Relativgeschwindigkeit in der Wasseroberfläche, die 
scherende Grundströmung die Amplituden des relativ zum Schiff stationären Wellenfeldes bei 
Bergfahrt (Fahrt gegen Strom) verkleinert und bei Talfahrt vergrößert, verglichen mit den ent- 
sprechenden Amplituden in stehendem Wasser. 


\ 
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Diese sich aus der Theorie ergebenden Änderungen des Wellenwiderstandes und des Wellen- 
feldes durch die Grundströmung stehen in guter Übereinstimmung mit den Ergebnissen von Ver- 
suchen, die mit einem Modell mathematischer Form durchgeführt wurden. Hierüber, wie auch 


gebenen Formeln, wird in einer späteren Arbeit berichtet werden. 
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Manuskripteingang: 26. 3. 1958 
Anschrift: Dr. F. Kolberg, Würselen bei Aachen, Kaiserstr. 69 


Über eine Rand- und Mittelwertaufgabe der 
Poissonschen Differentialgleichung”) 
Von Friedhelm Eicker**) 


In einem zylindrischen, durch zwei Querschnitte begrenzten Gebiet beliebiger Grundfläche ist eine Pois- 
sonsche Differentialgleichung Aw(x, y, 2) = f(x, y, 2) gegeben. Die Lösung soll am Rande nur von der 
Zylinderhöhe abhängen und auf den Deckflächen verschwinden. Außerdem soll ihr Mittelwert auf jedem 
Querschnitt des Zylinders verschwinden. Es wird gezeigt, daß das Problem eindeutig lösbar ist, wenn f(x, y, 2) 
stetig und zweimal nach z differenzierbar ist. Für das betrachtete und auch für allgemeinere Gebiete wird 
die Aquivalenz des Problems mit einer Integralgleichung 1. Art nachgewiesen. 


The Poisson equation Aw(x, y, 2) = f(x, y, 2) is to be solved in the domain between two cross-sections 
of a right cylinder with an arbitrary base. On the surface, the solution is required to depend only on the 
height z. Furthermore, the mean taken over any cross-section is required to be zero. It is shown that the 
problem has a unique solution if the function f(x, y, 2) is continuous and if the partial derivative fz, exists. 
For the domain under consideration (and more general ones) it is shown that the problem is equivalent to 
an integral equation of the first kind. 


Entre deux sections transversales d’un ceylindre quelcongue une equation differentielle de Poisson 


> 
Ao(z, y, 2) = f(x, y, 2) est donnee. A la surface du cylindre la solution ne doit dependre que de la hauteur 
du corps, sur les sections limitantes elle doit disparaitre. En outre, sa valeur moyenne sur chaque section 
transversale du cylindre doit disparaitre. Il est demontr& que le probleme a une solution unique au cas d’une 
fonction continue et doublement differentiable vers z. Pour le domaine considere et des plus generaux V’equi- 
valence du probleme a une Equation integrale de premiere espece est demontree. 


B ImAHHApnyecKof 06NACTH HPON3BOJIBHOTO TOIIEPEUHOTO CEYCHHA, 06a OCHOBAHHA KOTOPOoH 
TePIIeHANUKYAAPHEI K 06Ppa3ylollmM IMMHNPpa, 3amamo 1mddhepenmmaspbHoe ypaBHenne.1Iyac- 
coHa Aw(z, y,2)= f(x, y,2). IlpennosoraetcHa, uUToO pelteHme Ha Kpalo 06NACTU 3ABHCHT TOJIBKO 
OT BbICOTbI IMJIUHAPA MH HA OCHOBAHNUAX OÖ6palmmaeTca B HyAb. KpoMe Toro TpeöyerTca, YTOOBILHA 
KAKIOM IOMEPEeYHOM CeyeHHN IIMIUHAPAa 06Pamaloch B HYJIb TaR;ke u CpenHee 3HayeHne. 1Io- 
KAasbIBAaeTcA, YUTO MIA CYIUCECTBOBAHHS MU OJIHOBHAYHOCTH PeleHNA, NOCTATOYUHO HeipepbIBHOCTH 
u mmuhdhepeHnmmpyeMoctu 10 2 A0 BTOPOTO HOPANKA BRIINOYHTENBHO YIIEHA f(x, Y, 2), HApyIuaroımero 
ONHOPONHOCTB. „ln 60JIee 0o6MMX oÖJLacTel NOKABbIBAaeTca, UToO 9Ta 3amaya paBHoCHJIbBHA HEeRO- 
TOPOMYy HUHTeTPalbHOMY YPABHEeHWIO IEePBOTO POJa. 
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1. Problemstellung !) 


Gegeben ist ein zylindrisches, einfach zusammenhängendes Gebiet T beliebiger Grund- 
fläche zwischen zwei parallelen, auf den Mantellinien senkrechten Ebenen. Es werde ein kar- 
tesisches Koordinatensystem zugrundegelegt, dessen z-Achse zu den Mantellinien parallel sei. 


*) Über diese Arbeit hat der Verf. auf dem International Congress of Mathematicians in Edinburgh 
1958 berichtet. 

**) Mathematisches Institut der Universität Mainz. 

1) Die Aufgabe ist der numerischen Meteorologie entnommen; die Anregung dazu gab Herr Dr. 
K.Hinkelmann, Deutscher Wetterdienst, Frankfurt/M. 


über explizite Berechnungen des Wellenwiderstandes und des Wellenfeldes nach den hier ange- 


. 
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T werde in z-Richtung durch 0 und z, > 0 begrenzt. In T ist die Poissonsche Differential- 
gleichung 
Aofz, y, 2) = I(@, y, 2) 
zu vorgegebenem f(x, y, 2) unter folgenden Bedingungen zu lösen: 
a) der Mittelwert von ® soll auf jedem Querschnitt D des Zylinders verschwinden: 


= Jona: SE SZ: 


b) ® soll auf der Oberfläche M des Zylinders nur von z abhängen und auf der oberen und 
unteren Deckfläche verschwinden: 


o| = g() , g0) = ga) = 0. 
M 


g9(z) ist jedoch, abgesehen von den Randbedingungen, nicht vorgegeben und soll erst durch die 
Rechnung bestimmt werden. 


Es ist insbesondere zu untersuchen, ob und wann eine Lösung vorhanden ist, ob sie ein- 
deutig ist und ob bestimmte Forderungen für f(z, y, z) erfüllt sein müssen. 


2. Existenz und Eindeutigkeit der Lösung 
2,1. Zunächst ist in 


@*lz, y, 2) = SS! G&, y, 2,8, °) I& 1, d) ds dn dd 
zu jeder stetigen Funktion f(x, y, z) (auf schwächere Bedingungen für f wird verzichtet) durch 


Vermittlung der Greenschen Funktion G(, y,2,&,2,{) = — 1 + harmonische Funktion 
eine Lösung der Gleichung ar 
Aare} 


mit der Randbedingung &* = 0 vorhanden. Durch Mittelung ergibt sich daraus die mindestens 
zweimal differenzierbare und für z = 0 und z, verschwindende Funktion 


ve) = [] w*ta, y, 2) dx dy ee De ee 


Eine Lösung ® des Ausgangsproblems versuchen wir nun durch Addition einer in P harmo- 
nischen, auf dem Rande stetigen Funktion u(z, y, z) zu erhalten, für die gelten muß: 


a) = —ML) 2 200: m ee ae 
und 
b) um=Mz) WEM 2 a ee 


Gibt es eine solche Funktion u, so erfüllt 


o=o*+u 
offenbar alle gestellten Forderungen. 


2,2. Es seien zunächst einige für das folgende wichtige potentialtheoretische Sätze zu- 
sammengestellt. 

1. Ist u@x, y, z) eine in @ regulär harmonische Funktion, so läßt sie sich in jedem inneren 
Punkt P von G in eine gleichmäßig und absolut konvergente, beliebig oft gliedweise differen- 
zierbare Potenzreihe entwickeln. Diese ist also eindeutig und insbesondere identisch mit der 
Taylorentwicklung um P. Sie konvergiert mindestens in der größten in G gelegenen Kugel 
um P. 

Analog zum Satz über die analytische Fortsetzung in der Funktionentheorie gilt: 

2. Stimmen zwei in dem (einfach zusammenhängenden) Gebiet G harmonische Funktionen 
in einem Teilgebiet X überein, so sind sie in ganz G identisch. 

' Hierfür mag folgender Beweis gegeben sein. Angenommen, K sei ein echtes Teilgebiet, 
K das Komplement in G; dann gibt es einen Randpunkt P von K und K, der innerer Punkt 
von G ist, und eine ganze in G@ gelegene Kugel um P. In dieser gibt es eine zweite Kugel Q, die 
noch Punkte aus K enthält, deren Mittelpunkt P, jedoch innerer Punkt von K ist. Nun sind 
aber die Taylorentwicklungen der beiden Funktionen um ?, identisch und stellen ihre Funktio- 
nen nach 1. in der größten in G gelegenen Kugel um P, und damit in Punkten von K dar. 
K muß also leer sein. 
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Es sei angemerkt, daß sich eine harmonische Funktion mindestens bis an die Oberfläche 


der Konvergenzkugel der Reihenentwicklung um einen beliebigen inneren Punkt des Regulari- 


tätsgebiets eindeutig fortsetzen läßt und daß mindestens ein Divergenzpunkt auf dieser Kugel 
auch eine irreguläre Stelle der dargestellten Funktion ist. 

3. Ist nunmehr u(z, y, z) eine in dem in 1. bezeichneten Zylinder T harmonische Funktion, 
die auf den Deckflächen verschwindet, so läßt sie sich über diese hinaus in ein gewisses Gebiet L 
nach dem Spiegelungssatz der Potentialtheorie harmonisch fortsetzen. Spiegelt man das ganze 
Gebiet T an der unteren Deckfläche und ordnet Spiegelpunkten entgegengesetzt gleiche Funk- 
tionswerte zu, so erhält man in dem gespiegelten Zylinder T wieder eine dort harmonische Funk- 


tion ü&, y,2) = — u(x,y, —2). In dem Gebiet LCT stimmt ü mit der aus T fortgesetzten 
harmonischen Funktion u überein und stellt somit nach 2. die eindeutige Fortsetzung von u 
in ganz T dar. Offenbar läßt sich dieses Verfahren auf einen beiderseitig unendlichen Zylinder 
anwenden, in dem dann eine bezüglich z ungerade und mit 2 z, periodische, in z also beschränkte 
harmonische Funktion u entsteht. 


2,3. Die in 1. und 2,1. aufgestellten Forderungen sind den folgenden äquivalent, die sich 
unter der Voraussetzung der Harmonizität von u(x, y,z) ergeben. Ist u ferner auf M stetig — 
nach 2,2. gilt dies dann auch in dem unendlichen Zylinder über D —, so existieren für u und 
jede beliebige Ableitung bezüglich z gleichmäßig konvergente und differenzierbare Fourier- 
entwicklungen. Man hat 


u, y,2) = 3 u,@, y)sin”? . Be ern nn re EL) 
mit Be ; 
= | und sin Bee TAEN er (5) 
und z.B. 
pr Da u 


Die letzte Gleichung ist richtig, weil man in (5) die Differentiation nach x vor das Integral 
ziehen darf. (4) darf also gliedweise nach jeder Veränderlichen beliebig differenziert werden; 
insbesondere folgt die Vertauschbarkeit des A-Operators mit dem Summenzeichen. Für u, folgt 
daraus das Bestehen der Differentialgleichung - 


0° 0° ke ge \ 
= uU) | ud, LT HEHE 
(dat am a (6) 
Ferner konvergiert die Reihe 
=; ._ WIU2 
ge) = 3 g,- sin —, 
v—i 20 
wenn g(z) als stetig und stückweise glatt vorausgesetzt wird, gleichmäßig, woraus nach (3) ge- 
schlossen werden kann 
Ds constiauf F- 
unter /' die Projektion des Zylindermantels auf die x, y-Ebene verstanden. Schließlich lassen 
sich beide Seiten von (2) in gleichmäßig konvergente Fourierreihen entwickeln, denn auch 
die Mittelung kann man mit dem Integral in (5) vertauschen unter Erhaltung der zweimaligen 
Differenzierbarkeit nach z, welche die gleichmäßige Konvergenz garantiert, ebenso wie in (1). 
Somit folgt als dritte Bedingung für u, 


N a oe ER rend re 
y, aus 


En OT 
vl) = 2,7, : sin — 
v—1l “0 


Es läßt sich hier schon erkennen, daß eine Lösung des in 1. gestellten Problems falls vor- 
handen nur eindeutig sein kann. Denn die Differenz zweier verschiedener Lösungen wäre eine 
harmonische Funktion u(z, y, z) mit verschwindendem Mittelwert und nur von z abhängigen 
Randwerten, die sich wie oben entwickeln läßt. Sind die u,(x, y) ihre Fourierkoeffizienten, so 
gelten die obigen Bedingungen, jedoch ist statt (8) u, = 0. Nun ist 

vlg) ssin = = 
“0 


Mn vu 
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eine harmonische Funktion, die in dem Zylinderraum zwischen beispielsweise 0 und z - 


Randwerte und infolgedessen nach dem Mittelwertsatz auch Innenwerte von nur einem Vor- 
zeichen annimmt. Das ist aber mit u, = 0 nur verträglich, wenn u, = 0. Dies gilt für alle », 
infolgedessen ist auch u = 0, was mit der Eindeutigkeit gleichbedeutend ist. 


2,4. Zu den in 2,3. umgewandelten, neuen Bedingungen kann man nunmehr eine Lösung 
angeben. nen 
Für v, = u, — 9, lautet die Differentialgleichung (6) 


at) 0 +9: 


v, verschwindet am Rande, g, ist ein zunächst noch freier Parameter. Sind 9; 4 die zu 1 nor- 
mierten Eigenfunktionen bzw. Eigenwerte zu (6), so existiert die absolut und gleichmäßig kon- 
vergente Entwicklung für v, 


ra. St _Apde y) je 
v,(%, y) vg, _ 4— ( n]z)® 2) x 


somit 


2 u 
u, y) =9q,-\1— u . u: 
a | 20 & vrlzo)? —4;, 
Mittelung ergibt mit [dedy=F 
D 


EEE) =... (r- (3). in, )- Re 
u = al) —ı = ala) —A; > 
letzteres wegen der Vollständigkeitsrelation für die konstante Funktion 1: 
ey nz 
Da sämtliche Eigenwerte negativ sind, ist die Summe rechts in (9) stets kleiner als 0 und 
größer als — F; sie geht mit wachsendem »? gegen null. Gemäß (8) bestimmt sich g, nach (9) zu 


& Y 
9, A; ’ 


en (v fa a Pf 


D 


und u, wird 
“ 

u (X y) 2 ) PA (v 2) — FRANE 
re — 4 war Ki 
a 4; ; do 


Damit sind in eindeutiger Weise die Funktionen u,(z, y) gefunden, die den notwendigen Be- 
dingungen (6/78) genügen. Es ist noch zu zeigen, daß die mit diesen Funktionen gebildete 
Reihe (4) harmonisch ist; sie erfüllt dann nach Konstruktion die Bedingungen a) und b) aus 32,1. 

Nach dem Weierstrassschen Konvergenzsatz der Potentialtheorie konvergiert eine Folge 
von in 7 harmonischen, auf dem Rande stetigen Funktionen gleichmäßig gegen eine im abge- 
schlossenen Gebiet 7 harmonische Funktion, wenn die Folge am Rande gleichmäßig konvergiert 

Mit den Funktionen (10) ist für jedes » 

u,(z, y) - sin 
a! 

eine in 7’ harmonische Funktion, ebenso alle Partialsummen von (4). Auf der Oberfläche M 
nehmen die letzteren die Werte 


N 
Yr EN 
a | ‚sin 1 Ve ar he 


an. 
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9° 
worin die Konstante X = >> A ‚| ist. Die Summenfolge (11) ist also dann sicher gleich- 
i 


mäßig konvergent, wenn y(z) viermal stetig differenzierbar ist. Hierfür ist aber hinreichend, 
daß f(x, y,z) zweimal nach z differenzierbar ist. Damit sind die Voraussetzungen des obigen 


Satzes erfüllt. (4) ist somit in T gleichmäßig konvergent und stellt dort tatsächlich eine har- 
monische Funktion dar. 


Damit läßt sich als abschließendes Ergebnis formulieren: 
Die in 1. beschriebene Aufgabe hat für jede stetige, nach z zudem zweimal stetig differenzierbare 
Funktion f(«, y, z) eine eindeutige Lösung ® = ®* + u, worin u durch (4) und (10) und ®* in 
2,1. angegeben ist. Die anfangs unbestimmte Funktion g(z) wird dargestellt durch die gleich- 
mäßig konvergente Reihe 
= Y . 92 
= — b - sin 
9%) 2 A; Z. 


PIE en 


und ist stetig. 


3. Äquivalenz mit einer Integralgleichung 


Für die Anwendung ist es von Interesse, ob auch für allgemeinere als für zylindrische 
Gebiete 7 eine Lösung existiert, die entsprechenden Bedingungen wie den 1. genannten genügt. 
Diese Frage kann hier nicht beantwortet werden, jedoch 
wird die Äquivalenz des Problems mit der Lösung einer 
Integralgleichung 1. Art für die Funktion g der Randwerte 
aufgezeigt. 

Gegeben sei eine Schar glatter Flächen D(f) mit dem 
Scharparameter f, welche einen gewissen Raumteil stetig 
überstreicht und in der zwei verschiedene Flächen keinen 
Punkt gemeinsam haben. Zwischen den zwei Flächen D(t,) 
und D(t,) liege das einfach zusammenhängende Gebiet T. Bild 1 
Sein Rand R soll auf jeder Fläche D(f) mitt, St St, eine 
geschlossene Randkurve R’(f) mit nichtleerem Innengebiet erzeugen. Zwischen D(t,) und D(t,). 
soll R glatt sein. — Das praktisch interessierende Beispiel für T ist ein durch einen Kegel- 
mantel begrenzter Ausschnitt einer Kugelschale (Bild 1). 


Die Randwertaufgabe 
20 12.0, 2m 2, o= gt), auf R, 
in der die Randwerte nur vom Parameter { abhängen, wird gelöst durch 


26», 
An o(P) = — | «», 0) : [0 dg -! | GP, 0) 


en ING 


Ode u (19), 


P bedeutet den Punktx, y,z, Q den Punkt &,n,5, G(P,0Q) die Greensche Funktion des Ge- 
biets T. Die Bedingung 1. a) lautet nunmehr 


o’—= ([vod-=0, nee 
Dt) 


d.h. der Mittelwert von & verschwindet auf jeder Fläche D(f). Das liefert aus (12) 
CP, 
ONg 


|." 10a = | OO nn Haar. ar (13); 
iv R 
Die Vertauschung der Integrationsreihenfolge in dieser Formel ist erlaubt, wie man etwa 
auf Grund des Fubinischen Satzes sieht, wenn die Integrale über die uneigentlichen Integranden 
im Lebesgueschen Sinne als Integrale über den kartesischen Produkten T X D(t) bzw. RX Dh) 
aufgefaßt werden. — Die linke Seite ist eine Funktion von { allein und werde mit F(f) bezeichnet. 
In dem Integral rechts läßt sich das Flächenelement do folgendermaßen schreiben: 


I(s, t 
dor dsndl= Dan, 
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unter s die Bogenlänge auf den Randkurven R’(f) verstanden. I soll die Bogenlänge auf den 
orthogonalen Trajektorien zur Schar R’(t) auf R sein, deren Existenz wir voraussetzen wollen; 


2 B soll eine analytische, beschränkte Funktion auf R sein. Dann läßt sich (13) schreiben: 


RER Sei 
FÜ) = | I ask) dr. 


4 Re) H 


Dies ist eine Integralgleichung 1. Art für g(r) mit dem unsymmetrischen und nicht be- 
schränkten, aber integrablen Kern 


"OG(P, 0)” Ol(s, T) 
= — _— — er BE LET 7 N 
Kit, ı) J no BR 
Ir) 
Falls diese Integralgleichung 


& 
Fl) = [ Kit, r) g(r) dr 
F 
eindeutig lösbar ist, existiert auch ein » der verlangten Art. 


Der Verf. ist Herrn Prof. Grunsky, Mainz, für mehrere hilfreiche Bemerkungen sehr zu 
Dank verpflichtet. ö 
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Existenz periodischer Bewegungen eines n-fachen Pendels im Falle, 
daß einige Wurzeln seiner charakteristischen Gleichung 
ein Vielfaches einer anderen sind 
Von @. Bradistilov und G. Boyadjiev 


In einer seiner früheren Arbeiten beweist @. Bradistilov das Bestehen von n Scharen periodischer Be- 
wegungen eines Systems aus n hintereinander gekoppelten, in einer Schwingungsebene liegenden Pendeln, 
unter der Voraussetzung, daß keine einzige Wurzel seiner charakteristischen Gleichung ein Vielfaches einer 
anderen ist. In der vorliegenden Abhandlung wird dieselbe Frage betrachtet, aber unter der Voraussetzung, 
daß p— 1 Wurzeln der charakteristischen Gleichung ein Vielfaches einer anderen darstellen. Es wird 
bewiesen, daß für diese Wurzel überhaupt so viele Scharen periodischer Bewegungen bestehen, wie die An- 
zahl der reellen Lösungen eines die Bedingung der Periodizität ausdrückenden algebraischen Systems vierten 
Grades aus p— 1 Gleichungen mit p — 1 Unbekannten angibt. 


In one of the previous papers by one of the authors (G. Bradistilov) it has been proved that for a 
system of n coplanar pendulums coupled in series there exist n families of periodic motion provided that 
none of the roots of the characteristic equation is an integer multiple of any of the other roots. In the present 
paper Ihe same problem is considered under the assumption that p — 1 roots of the characteristic equation 
are multiples of another root. It is proved that for this root the number of families of periodic motion is 
equal to the number of real solutions of a system of p — 1 algebraic equations of 4th degree containing p — 1 
unknowns, a system that represents the conditions for periodieity. 


Dans un de ses articles precedants G. Bradistilov prouve lexistence de n systömes de mouvements 
periodiques d’un systeme de n pendules couples Vun apres l’autre et se trouvant sur un seul plan des oscil- 
lations, tout en supposant que pas une seule racine de son equation characteristigue ne soit un multiple 
d’une autre. Dans Varticle present la möme question est mise en consideration, en supposant cependant 
que p — I racines de l’equation characteristique presentent un multiple d’une autre. Il est prouve que pour 
cette racine en general autant de systömes de mouvements periodiques existent que s’elöve le nombre de 
solutions reelles d’un systeme algebrique du quatrieme degre de p — 1 equations avec P— 1 inconnues, 
syslöeme qui exprime la condition pour la periodicite. 


B o1HoN u3 CBOUX ÖoJlee panHuux paoor I‘. Bpammeru1oB HORAasbIBaeT CYINeECTBOBAHHE 
N CEMEHCTB IEPUONMYECKUX BUSREHMÜ CHCTEMBI, COCTOAIMEH U3 N HOCHENOBATEIIBHO CBAZAHHBIX 
MAATHUKROB, JIEFKAIUMX B OMNHOH INIOCKOCTU, B UPENIOJNOKEHHH, YTO HM OJMUH KOPeHb Xapak- 
TEPHCTUMECKOFO YPAaBHEHUA CHCTEMbI He ABJIAETCH KPATHbIM APyToro Kopnun. B Hacrosmei 
pa6oTe paccMaTpuBaeTcaH Tore BOIPOC, IPMMEM OAHAKO JNONYCKAeTCAH, YTO pP — 1 KopHeü 
HBJBHOTCH KPATHBIMU APYTOTO. JloRasbIBaertca, UTO MIA HTOTO KOPHA BOOÖMEe CYINECTBYET 
CTOJIBRO CEMEHCTB HEPHOJUYECKUX IBUSKEHUM, RAKOBO YUCHO NEeÜCTBUTENBHBIX pelueHnHunüä HeKROo- 
Topoä auıreöpanyeckoä cHCTeMbI 4-oü cTemeHun, coctronmei us p —1 ypasueunä c p — 1 
HEN3BECTHBIMH, BbIPA7RAlmImeh YCJIOBHE UePHOAHYHOCTH. 
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5 Aus der Arbeit [1] ist bekannt, daß das n-fache, in einer senkrechten Schwingungsebene 
liegende physische Pendel n Scharen periodischer Bewegungen um seine stabile Gleichgewichts- 
lage besitzt, wenn man voraussetzt, daß die charakteristische Gleichung des Systems der Be- 


wegungs-Differentialgleichungen des n-fachen Pendels Wurzeln aufzeigt, die einfach sind und 


keine unter ihnen als Vielfaches einer anderen erscheint. In [2] ist dieselbe Frage von dem 
zweifachen physischen Pendel behandelt, wenn aber die eine Wurzel der charakteristischen 
Gleichung ein Vielfaches der anderen ist. Es hat sich erwiesen, daß in diesem Falle überhaupt 
eine Anzahl von drei Scharen periodischer Lösungen besteht. 

Im Vorliegenden wird die allgemeinere Aufgabe gestellt, die Existenz periodischer Be- 
wegungen eines n-fachen Pendels nachzuweisen, vorausgesetzt, daß eine oder mehrere Wurzeln 
der charakteristischen Gleichung Vielfache einer der übrigen sind. 


$ 1. Bezeichnen wir mit o, den Winkel, den das »-te Pendel mit seiner Gleichgewichtslage 
schließt, mit /, das Trägheitsmoment desselben Pendels in bezug auf seinen Schwerpunkt; mit 
a, den Abstand zwischen der Drehachse und dem Schwerpunkt des v-ten Pendels, mit b, den 
Abstand zwischen den Drehachsen des » — 1-ten und »-ten Pendels. Setzen wir 


n N n 
2 1 1 
Adam, =D" >. M,, B, =4,m. DD, m, Vl, 2er) 
k=v-+1 k=v+1 


und 9, =/y, wo 4 ein kleiner Parameter ist, so erhalten wir aus [1, (8)] für die Bewegungs- 
Differentialgleichungen des n-fachen physischen Pendels 


u—I n 
Er leo, Zu) FA y, Ze) et AT Beles 
f 2n—3 ya 


i £ ® 1 
+4 sind (p, — y) vl = —9B,y, +AgB,% ( 1)% One! W=|1,...;,n) ). 


n=2 


Speziell für A = 0 ist (1) ein Erregungssystem 
v„—1 


NE re ee EEE ION 
= k=v-+ 


dessen charakteristische Gleichung nach [1] lediglich rein imaginäre Wurzeln + g'(k=1,...,n) 
aufzeigt. 
Nehmen wir an, daß p— 1 dieser Wurzeln Vielfache von o, sind, beispielsweise 
0 = k;0 FEED EHE,D): 


Wo ky,ky...,k, ganze positive Zahlen sind. 
Das gesamte Integral des Erregungssystems (2) ist 


w=-_L%L, (C,c0so,t-+D, sin 0,0, 
ve N ee) 


wo = Zu en G, Qu sin Qu e. D,; 0. €08 0, l) 


H= 


wo L,, (u ey n) die Gleichungen 
v—1 N 
B, On 2) b; Dis + (A, + J,) Qu a B,] 0 4 b, 0 = ‚Bi Lea =) (v = L SCH! n) (4) 
1 k=v-I- 


erfüllen und die dadurch gebildete Determinante A = |L,,| verschieden von Null ist. 
In [1] stellt man fest, daß für jede Wurzel 0, ($s=2,...,n), wenn 0,/o; kein Ganzes ist, 
bei kleinen Werten von A je eine Schar periodischer Bewegungen des n-fachen Pendels besteht, 


welche Bewegungen von einem Parameter, dessen Periode bei A — 0 entsprechend zu n neigt, 

Ss 
abhängig sind. Bei vielfachen Wurzeln, beispielsweise wenn 0; —= %;0,, ist die Funktional- 
determinante — sin "2... sin Tr 4 gleich Null und der Beweis von der Existenz perio- 
discher Lösungen für de le kommt außer Kraft. 


$2. Unsere Aufgabe besteht darin, den Beweis anzutreten, daß auch bei der Wurzel o, 
periodische Lösungen bestehen. Bei Aufstellung des Beweises wollen wir von einer Anmerkung 
aus |1] Gebrauch machen, nach der die Bewegung periodisch ist, wenn das n-fache Pendel 


bei seiner Schwingung zweimal durch die Gleichgewichtslage geht. 
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Wollen wir zu diesem Zwecke das Partikulärintegral des Erregungssystems näher betrachten: 


2 ” 
vO=L,sngat+3%L,N,sink,at; 
Pr A EN en‘ 


p 
vol) =04L,00sat+ 0% Luk, N, cos k, 0, t 
= 


dessen Periode = ist und die Anfangswerte 
1 


4 
Y,0(0) =0, %0(0) == L,; f 0 2, Ir k, N, R 


sind, wo N, Na ..., N, als unbekannte Konstanten auftreten. Die letzteren sind so zu er- 
mitteln, daß die entsprechende Lösung des Systems (1) periodisch wird. 
Das Integral des Systems (1) mit den modifizierten Anfangswerten 


p n 
0)=0, WO)=-aLı t+aLL,kN, ta) + & In 3 
u=! u= 


kann für hinreichend kleine Werte von A und a, (n=23,...,N) mit Potenzreihen von / und 
&, in folgender Gestalt 


v0 =) + EP + RIO + RN +] @=h en): 


dargestellt werden. 
Hier sind P,,(), Q,(), R,,() usw. unbekannte Funktionen von £, welche für t = 0 auf 
Grund von (6) folgende Anfangswerte annehmen: 
P, (0) = 0, P, (0) = Ou L,,» Q,(0) = Q,(0) ve 0, R, „(0) > R, (0) => 0. 


Wenn man das Integral y,(t) aus (7) ins System (1) einsetzt und die Koeffizienten vor 
den entsprechenden Potenzen von «, und / gleichsetzt, so erhält man für die unbekannten Funk- 
tionen P,,(d, Q,(b), R,,() usw. Linearsysteme von Differentialgleichungen 


v—1 n 
B, 2 br Yk m (A, tr J,) 2. 2 In = Ken, Yk Tr g B, A F,() nV (8), 
— =r+ 
wo F,(t) für die Funktionen P 


vu 


(D, Q,(t) und R, (f) entsprechend gleich sind: 
F+O=0, (Beh 


® v—1 1 ir E 
EA. ==4B, _ +B, 2 de 5 (Wo — Yeo)* Pro — Bro — Yo) vi 


N 1 j k h 
+4 BR Bi 5 (dr — Pro) Yo — (Pro — Yko) Yko : , » (10), 
R OFKt) 
I P. l = u = 2 (00 5 v= er na 177° = . . . . . . . . 
(tl) IN, (u PP b=1 ‚n) (11) 


Das homogene System der Differentialgleichungen von (8) ist gerade das System (2), sein 
Gesamtintegral ist also (3). In Anbetracht der Anfangsbedingungen der Funktionen P, „(t) so- 
wie dessen, daß diese Funktionen Partikulärintegrale von (2) sind, findet man 


P,,() =L,,sin o,t Vals) Amer ee 


Die Funktionen Q,(l) und R,,() werden sich ebenso wie die Partikulärintegrale des Sy- 
stems (8) mit Anfangsbedingungen y,(0) = y,(0) = 0 ergeben, deren rechte Seiten durch (10) » 
und (11) gegeben werden. Um das Partikulärintegral von (8) finden zu können, werden wir 
das Verfahren der Variation der Konstanten anwenden. Dieses Integral wird von der Art (3) 
sein, wo die Ableitungen von C,(l) und D,() (u=1,...,n) das Gleichungssystem erfüllen 

F,() 


Bi; ti > L, N ».: ‚ 
2 ur (C, cos 0,E + D,sin 0,0 ei 2 e (— C, sin 0, + D,cos 0, = r 7 (13). 
4 u=l1l [2 


v 


Die aus den Koeffizienten der Unbekannten gebildete Determinante wird am raschesten 
berechnet, indem man {= 0 annimmt, weil aus dem System der Differentialgleichungen (2) zu 


. . . — b 
ersehen ist, daß sie eine Konstante darstellt. Sie beträgt D=— Kan — 


en 


BIZUrT irn x 


Hr f u BESTEN TÜREN 
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Nachdem gewisse Umwandlungen nach Spalten vorgenommen werden und die Laplacesche 
Regel angewendet wird, nehmen die Determinantenwerte in den Zählern der Lösungen des 
Systems (13) folgende Gestalt an: 


— singt. A ar < FF, 

ö En N SE EN EN I = 10.B (el n) 
u cosg;t- A 20 a 
‚= 


++. Qi-1' Bi4l On gB, 
wo 4,; der zum Element L,; gehörige Minor in der Determinante 4 ist. 
Dann sind die Lösungen des Systems (13) 


v d; 4; ie 2 d; 4, ; . 
en singt, D=,- oosgit (Del) 
FR, 
wobei A, = ——_ Au: 
f , 2 gB, # 


Hieraus ermitteln wir C,; und D, und erhalten für das Integral des Systems (8) mit An- 
fangsbedingungen y;(0) = y;()) =0 (i=1,...,n) nach entsprechender Umwandlung 


t 
1 .. i 
v.= 4 | |&, 0022.41 sine, e—w du a N: 
0 


Die gesuchten Funktionen Q, und R,, ergeben sich aus (14), indem man F‘, entsprechend 
durch die Ausdrücke (10) und (11) ersetzt. 

Für genügend kleine Werte ng A| es [&.| erscheint das er y,(t) als periodische 

2n +2Ö 


Funktion mit modifizierter Periode er Integralperiode = (5), wenn 
1 


d.h. wenn 


p p 
—L,,sinöfF &L,,N, sin k,6 = 2,2, sink,d.&, 


u=2 


Fr 5 L, „sin ea +6) ou ‚ale ) ; 3 Rn ar wi 
Ak=Pp-+1 ur 2 


Die Zeichen F entsprechen dem geraden oder ungeraden Ak, 
Nachdem wir nach den Potenzen von ö weiterentwickeln, erhalten wir das System 


3 tm +24 94-0, 08) 
! 1 


ee) 


1 n 
< TOQu 


n p 
€, ar: > L, u I, nn Une > L, u k, & E En = Di L, u 93: —— Aus 9, =. > I, Au \ 
u=2 u=2 u=p+l 01 u=2 
Die eingeklammerten Glieder enthalten ö wenigstens in erster Potenz, die Glieder außerhalb der 
Klammern in höherer Potenz. 
Das System (15) ist gegenüber einer Anzahl von n — p + 2 Unbekannten a1, + - » Om 
6, 42 homogen. Um eine Lösung zu erhalten, soll die Matrix 


Ines: allen Det ga) = 9, | 


| 
| Inp+ı ua Inn Imı min (en er Nn) ln er 7 9, 

von niedrigerem Range als n—p-+ 2 sein, d.h. sämtliche n — p + 2-reihigen Determinanten 

dieser Matrix sollen gleich Null sein. Wenn aber p — 1 dieser Determinanten Null sind, soistes klar, 

daß auch die übrigen Determinanten gleich 0 sind. Folglich, damit das System (15) eine Lösung 

hat, genügt es, daßp— 1 (n—p+ 3)- reihige Deter minanten gleich Null sind. 
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Dies seien folgende Determinanten: 


Lip+41 ...Lın L.3&G+m)—& Q+09ı | 
x Into ... In»+1,n In-+1, 1 + (in p+1 7 No) nr On—+1 + Oa-r41 
Lsp+1 .. Lion Lu &* (+ Ns ) —ı Q, + 9, 
(=n—p+2,...,n) . (16). 
Nehmen wir nun an, daß die Determinante 
| In p+ı | Lin Li Tr & 


La +1,pP+1°.* In—»+1,1 Inp+1.1 E I p+1 | 


verschieden von Null ist. Dann können wir aus den ersten n— p + 1 Gleichungen des Systems 
(15) &p41» - - - » &m und ö finden, die folgende Gestalt haben werden: 

ER (op: 00); — A 1 00); (=p+Ll..n.. 
Diese Werte setzen wir in die Determinanten T, ein. Es ist klar, daß die Gleichungen (16) für 
genügend kleine Werte durch 5,(A), .. . , %,(A) erfüllt werden können, wenn ihre freien Glieder Nullen 
sind und die durch die Koeffizienten von a, ..., a, gebildete Determinante verschieden von 


Null ist. 
Die freien Glieder sind 


Lip+1 ... Lin Lu # &ı Qı 
ER ae ee an 2. RA ee Re ee 
R Un oe ... In p+L,n In s+H1 + En On—+1 ( ) 
Lsp+1 Au Lyn L;ı + & Q, 
(s=n—p+23...,n) 
und die Koeffizienten von a, ..., &, Sind 
OT, 
0,, = s=n—Pp + 23..,0D, WER. .25BDL 
su ON, ( Air, P) 
Die aus obigen Koeffizienten gebildete Determinante soll verschieden von Null sein: 
|or | 
ENTE 
ul = 5m, #0 0 en Arne x ee 
Die Bedingungen (17) stellen eine Anzahl von p— 1 Gleichungen gegenüber den unbe- 


kannten Konstanten N,...,N, dar. Um diese Gleichungen finden zu können, sollen wir Q, (= 
[4 
welche Funktionen von N,, ..., N, sind, durch die Ausdrücke (14) ersetzen. Nach den entspre- 


N 
chenden Berechnungen ergibt es sich, daß (2) Polynome dritten Grades in bezug auf N,, 
Ne end, Eı 


In (17) setzt man Q,, ..., Qu-p41, Q, ein und so wird ein System aus p — 1 Gleichungen 
von p — 1 Unbekannten vierter Potenz gegenüber N,, ..., N, erhalten. 
4 
SEE „Om hehe NE Na... Nr=0 (++. + „s4 (m=1,...,p—1) (19. 
Pan. ER 
p 


Dieses System hat höchstens 4?! reelle Lösungen für N,..., N,. Eine weitergehende Unter- 
suchung dieser Wurzeln ist fast undenkbar, weil die Koeffizienten Ay isig...ip Komplizierte Ausdrücke 
sind und die Anzahl der Gleichungen p — 1 ist. Trotzdem könnte man Folgendes schließen: 

Wenn das System (19) eine Anzahl von s reellen Lösungen Nas: +, Ny, zuläßt, wo s S 4-1 
und wenn für diese Lösungen die Ungleichheit (18) erfüllt ist, so bestehen für das n-fache physische 
Pendel in der Umgebung der u s Scharen periodischer Bewegungen für die 

a 


Wurzel o,, mit modifizierter Periode ‚ welche durch die Annäherungsgleichungen 


01 


: 2 . 
L,singt+ 3 L,N,sink, gt 


u=2 


Pl) = A YA 


’ 


NENEWEN.*,. 
9, d —4 Yrs(t) wA 01 


ausgedrückt werden. 


» 
L,,csg2+3L,,N,,k, cos k. 0 ı 


u=2 
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. Insbesondere können etliche oder sämtliche der Zahlen k,(u=2,..., p) gleich 1 sein. Be- 
zeichnet man mit r die Anzahl dieser Zahlen, so besitzt die charakteristische Gleichung des Sy- 
stems (2) eine r-fache Wurzel. Ist für diese Wurzel die Ungleichheit (18) erfüllt, so bestehen für 
das n-fache Pendel s Scharen periodischer Bewegungen für die Wurzel 0- Wennr<p—1 
werden diese Scharen durch die Annäherungsgleichungen ausgedrückt: 


9,:(d = 4 Y,(d ND A [(Z, ı ur > Ir Nu) sin al Eile 2 Lou, Nags sin RK. 01 l, 1 
9) = Ay) Bra la + I L,N,)osgt+IL,N,.sok,adl J. (21), 


wo u, die Nummern der Wurzeln, betragend o,, und u, die Nummern deren, die nicht gleich 0, 
sind, bezeichnen. 


Wenn r=p—-1, so sind die Annäherungsgleichungen 


p 
0) = A Y,s(d) > 7 IL, ı = 2 Lu N,sl sin 01 t ’ 
u= 


. ’ DE 
9% =Amd) RrAalLı+ 3% L,.N,sl cos ot 
u=2 


$3. Wollen wir nun den Sonderfall betrachten, wo die eine Wurzel der charakteristischen 
Gleichung ein Vielfaches einer anderen ist, beispielsweise 0, = k o.. 

Dann wird das System (19) zu einer Gleichung dritten Grades in bezug auf den einzigen Para- 
meter N reduziert: | 


PN NEFDNFBEN Fps Vene er (23), 
und die Bedingung (18) nimmt die Gestalt an: 
dP(N) 
an re ee ee (24) 


Nach Durchführung der entsprechenden Berechnungen ergibt sich, daß der Koeffizient p; 
beikZ# 1,3 gleich Null, und bei k = 1,3 verschieden von Null ist. Man kann also folgende 
Schlußfolgerungen ziehen: 

1. Beik # 1,3 besitzt die kubische Gleichung (23) eine Wurzel Null, weil p, = 0. Durch 
direkte Berechnung wird für diese Wurzel nachgeprüft, daß die Ungleichheit (24) erfüllt ist. Dar- 
aus können wir schließen, daß bei k # 1,3 eine Schar periodischer Bewegungen des n-fachen 


Pendels in der Umgebung seiner Gleichgewichtslage mit modifizierter Periode en besteht, 
welche Bewegungen durch die Annäherungsgleichungen 9ı 


ein) SAL, sinot, BANDS Au Le, Cosa. ee Teer nun) 
ausgedrückt werden. 
Sind die restlichen Wurzeln der Gleichung (23) reell und ist für sie die Ungleichheit (24) er- . 


füllt, so bestehen für die Wurzel 0, noch zwei Scharen periodischer Bewegungen für das n-fache 
Pendel mit den Annäherungsgleichungen 


9) = Ay) & All, sing t+ L,,Nsin kg), 
od) =Ay) ZA (L,ıcosgt+kL,N coskod) 


2. Beik =3 ist das freie Glied der kubischen Gleichung überhaupt verschieden von Null. 
Dies besagt, daß in diesem Falle die kubische Gleichung (13) wenigstens eine von Null ver- 
schiedene Wurzel hat. Folglich, wenn für diese Wurzel die Ungleichheit (24) erfüllt ist, so be- 
steht eine Schar periodischer Bewegungen für das n-fache Pendel um die stabile Gleichgewichts- 


lage mit modifizierter Periode ae ‚ deren Annäherungsgleichungen von der Art (26) sind. 


Ger) 5 ar 


1 
Sind auch die übrigen Wurzeln der kubischen Gleichung reell und außerdem einfach, so hat man 
noch zwei Scharen periodischer Bewegungen des n-fachen Pendels, deren Annäherungsgleichungen 
ebenfalls von der Art (26) sind. 

3. Beik = 1 besitzt die charakteristische Gleichung eine doppelte Wurzel 0, = o,. Da das 
freie Glied der kubischen Gleichung (23) überhaupt verschieden von Null ist, so besitzt diese 
Gleichung wenigstens eine reelle Wurzel, die verschieden von Null ist. Ist also die Ungleichheit 
(24) für diese Wurzel erfüllt, so besteht eine Schar periodischer Bewegungen des n-fachen Pendels 


2 ER 
um dessen Gleichgewichtslage von modifizierter Periode erh mit den Annäherungs- 
gleichungen gı 
‚D=iyl) <A(lLı +L.N)singt, 
N) ri FE SR 
AG) a Ay) N 40 (L,ı 2; L. N) C0S 0] t 
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Sind auch die übrigen Wurzeln der kubischen Gleichung reell und einfach, so werden noch 
zwei Scharen periodischer Bewegungen vorhanden sein, deren Annäherungsgleichungen ebenso 
von der Art (27) sind. 


Anmerkung 


Das Verfahren kann auch auf ein System aus physischen, in einer Schwingungsebene 
liegenden Pendels angewendet werden, welches mit gleichbleibender Winkelgeschwindigkeit (0) 
um eine außerhalb der Ebene liegende Achse rotiert. Dabei bleiben die Ergebnisse dieselben. 
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Torsion of a large thick plate with rigidity varying linearly 
with depth and rigidly fixed at one face 
By S. B. Dutt 


In der vorliegenden Arbeit wird das Torsionsproblem für eine ausgedehnte dicke Platte behandelt, deren 
Steifigkeit in Richtung der Dicken-Koordinate veränderlich und deren eine Deckfläche festgehalten ist. 
Die Untersuchung wird auf den Fall ausgedehnt, daß in beiden Deckflächen vorgegebene, in radialer Rich- 
tung veränderliche Schubkräfte wirken. 


The present paper deals with the torsional problem relating to a large thick disk the rigidity of which 
varies in the direction from one face to the other, one face being assumed fixed. T'he investigation is extended 
to the case where both faces are acted upon by shearing forces depending on the radial coordinate. 


Dans Varticle present on traite le probleme de torsion pour une vaste plaque Epaisse dont la raideur 
dans la direction de la coordonnee de l’Epaisseur est variable, et dont un des deux recouvrements est fixe. 
L’examen est elargi sur le cas que des forces axiales donnees et variables dans la direction radiale agissent 
aux deux recouvrements. 


B Hacroamei pa6orte Hccsenyerca 3anaya Kpy4eHun Mı1A 60NbMOÄ TOJICTOÄ ILTACTHHKH, 
FRECTKOCTb KOTOPOH H3MEHNHETCH B HAIIPABJIeHUN KOOPAHHATEI TOJINMHBI H BEPXHME OCHOBAHHE 
KOTOPoüÜ 3aKpelsIeHH0. Mlccse1oBaHne pacıpocTpanserTch Ha cıIyyali, KOT1A B 000HX OCHOBA- 
HUAX NeÜCTBYIOT CHBUTAIOIIME CHJIbI, H3MEHSIOIMMeCH B PANHANbBHOM HaupaBJIeHHn. 


1. Introduetion 


Investigations on torsion of a semi-infinite homogeneous elastic medium by a shear stress 
were carried out by Reissner and Sagoci [1] on considering a mixed boundary value problem 
in the theory of elastic equilibrium and they obtained an explieit solution by introducing in a 
suitable manner a system of oblate spheroidal coordinates. The same problem was reduced by 
Sneddon [2] to a pair of dual integral equations by the use of Hankel Transform and a solution 
was obtained for the static case. Yi-Yuan Yu [3] investigated the problems of torsion of a 
semi-infinite body and a large thick plate of finite size on expressing the applied shearing 
force in terms of aFourier-Bessel Integralinstead of a series, the latter being used in the case 
of a bar of finite size. 

In the present paper we want to discuss the similar problem of torsion of a large thick 
plate with rigidity varying linearly with depth and rigidly fixed at one face. The analysis 
has also been extended to the torsion of a plate by two shearing forces prescribed at the bound- 
ing surfaces. 


2. Formulation of the problem and boundary eonditions 


We are to determine the components of the stress and displacement vector in the interior 
of an elastic solid with varying rigidity bounded by planes z = Oandz = hwith positive direction 
pointing into the body. At the boundary of the body z = 0, torsional shearing force is prescribed 
over a circular area with its axis coincident with z-axis, while the face z = h is held fixed. 

It has been shown by Michell that an isotropie bar of circular cross section of varying 
diameter when twisted by terminal couples, only the tangential component v of the displacement 
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vector is different from zero and all other components of stress vanish identically except two 
which on using cylindrical coordinates (r, ©, z) are given by the relations 


OD a) 
u Tem ü Se z EL Re ee ARE, [09) 


u being the shear modulus. 
In our problem the medium being one of varying rigidity, u is expressed in the form, 


Pe TR a nr (2). 
Eq (1) then assumes the form 


= 0 
To; — (lo + Mm 2) = To = (lo + 1 2) | a a 8). 


Inserting the expressions of Eq (1) in three equations of equilibrium shows that two of them 
are identically satisfied and the remaining one 


_ er Tg; 
I 


gives on using the Eq (3), the following differential equation for v, 


O®v 1 v 0°Vv u er 
a a ne ng (5). 


To obtain the boundary conditions, we assume the shearing stress at the boundary z = 0 to 
be acting within a circular area of radius a and is given by the function, 


2 T,o 
02 ar T 


N EEE (4), 


Zi, 20, T, = Fr) =kKr), 
ZN, R>.@, ee er (6) 
and ZN DE) 


k being a constant. 


3. Solution of the problems 
To solve Eq (5) we assume, 


Re en u re Re. 
The Eq (5) then assumes the form 
Br ZUR ERZ az u dz 
- ———- —— — N) Ta ea 8); 
dr? r dr ee ® 


whence we get, 


a ee Ba 1 5 
— + — — en Re Oi A are: 9% 
er k m R S 
and 
d2Z m dZ : 
——, _- BE a 10), 
tm Ft = 
where 42 is the separation constant. 
The solution of the Eq (9) is 
Rerılir) (11), 
where J, is the Bessel function of the first kind and first order. Putting 
EIERN (12) 
Yı 


we have, dö = dz and the Eq (10) becomes 
d2Z 1 dz 


AN ee Pe 13), 

de 4 Ede RZ=0 (13) 
whose solution is given by 

va Fa A ae Was. LE Va ee Er Be Er (14) 


where /,, K, are modified Bessel functions of the first kind and zero order and of second kind 
and zero order respectively. 
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Thus the final solution of Eq (5) is obtained on using Eqs (11) and (14) in Eq (7) 


Jar) [Ay I(A2) + By Kol L)] dA 


i E „Gay 
= j Jar) IA, Ih ee) +B, K,(a® u) Alla 
0 


From Egs (6) and (15), we have atz=h 
(„#0 ta2) B,K er) ee En (16) 
A, lo (A HE +D,%0 m 
and using Eq (15) in Egs (3) we get 
5 4m 
ee | Jar) f AT, (a Ko a) 738 u) d....(dm, 
0 


Ne 1) Jar) h As (a a +1B,K, (a ee) 41. (U 
0 


The shearing stresses and displacement at z = O0 are given by 


u | 1 Jar) IA, %; ( Fa B,K, ( ze) ENT: (18), 
een, | 1. Jr) IA, 1? ( ke) RR, (a =) di (18.1), 
(e EX n) IA, fl ( ke) EOOR, (a ze) PN (18.2). 


The coefficients A, and B, are, in general, functions of the parameter and may be determined 
from the boundary conditions (6), i. e., (7e,)z-o = F(r), F(r) being a function of r. Expressing 
this in the form a Fourier-Bessel Integral [5] as 


FO-SIALANAyEWIRDM seen (19) 


and comparing the same with the expression in Eq (18) and using boundary conditions in Eq (6), 
we obtain 


ie I ka 
Ana) BR, (a) I | Ryıaa ee 
atı ih #1 ih Bu y® J(Ay) dy ag 2(A a) (20) 


From Eqs (16) and (20), A; and B; can be obtained as 
SIR Jy(A a) Ky(A hy) 


A, = R 7 (21), 
AR (ak) Kan) + Idah) RK, (akt) 
u kı 
2 
Bu = JA 9) Jul iu 2 SE (22), 
MEN Eu) KA) + Ih) K, (a = 
Kı Mk’ 
where 
hı = Ko af Kı h 


Hı 


with the use of these coefficients in Eqs (15) and (17), 


we obtain the displacement component v 
and the stress components 7g,, 


Zub 
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4. Torsion with shearing forces on the two bounding planes 


The foregoing analysis can be extended to the case of a large thick plate undergoing a 
torsion with prescribed torsional shearing forces acting on the bounding planes, z—= + h, where 
the half-thickness h of the plate can have any finite value. 


The boundary conditions are 


Ge) = Fr), rsa, (ae nee ee (23) 
in which F and G are given functions of r. 
From Eq (17) stresses at z= + hand z= —.h, can be written as 


Gola t+ mh j AIQ@r) [A LA) — B,K,(@h)] dA, 


R ESP BE (24) 
Treo) = — (wo tm) SAJsAr) [A, IA h,) + Bi; Ky(A h,)] dA 
and z 2. 
(Toz)z= —n = (u — u h) [A Jı(A r) [A, IA h,) — B, Ky(A h3)] dA, s 
Si (25) 
Ho)-—n = — (WM) {4 Jar) [A, Io(A he) + B, Ko(A ha)] dA 
where h, = We } 


U 
Using the Fourier-Bessel Integral Eq (19) and remembering the boundary conditions 
Eq (23), we obtain 


(ode = En) =/3 Jar) aa BE HONd (26), 
ae ii ALarnd J YOU een: (27) 
whence we have, 
Yı hı [A, IL@h)— B,K,(Ah)] = N EDIT UdU ee (28), 
ke, A; hab) Bah)] = Sy AAMay 22: (29). 


From Egs (28) and (29) coefficients A, and B, can be obtained and substitution of the same in 
Egs (24) and (25) gives us the expressions for stresses acting on the bounding planes of the thick 
plate. 


5. Discussion 


(i) To get the solution of torsion problem of an isotropic large thick plate rigidly fixed 
at one face subject to a torsional shearing force prescribed over a circular area at the free sur- 
face with its axis coincident with z-axis, we have to make in Eq (17) after substituting for A, 
and B, from (21) and (22) noting that for large value of x 


1 


I(a) = ae ae Pe: 0) 
Van 
we then have, making u, > 0 
Le" [ ern h) 
o,=ku | 72a) Ar) osh (A) I RN (31), 
Be | IA a) Isar et GR 5 U ER 1.1), 


sinh} (2 —h 


>) 
ee An (31.2 


_— er (Aa) A(Rr) 


which agree with the results of Yi-Yuan Yu [3] as given in Eq (16), if we remember that he 
took the origin at the face rigidly fixed, 
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(ii) To get the solution of torsion problem of an isotropie semi-infinite body subject to 
the torsional shearing force over a circular area at the face z = 0 with the positive z direction 
pointing to the body, we should note that only k functions can occur in the expression for v 
as given in Eq (15) and the second boundary condition (20) need be satisfield after substituting 
for B, from Eq (20) in Eq (17) and making u, > 0 and using Eq (30), we get 


1, ka aD Tuner N . (82), 
10] “ 
ka] La In ee (32.1), 
j 0 
ka Ä 
ae | Lad Lüne MM... 2.822) 
es 


which agree with the results of Yi-Yuan’ Yu [3]. 


6. Numerieal Caleulations 


To calculate the shearing stresses and displacement we take for example the case of the 
earth. Seismological obsvervations (Bullen [6]) reveal that within the earth the value of u 
varies from 0.63% 1012 dynes/cm? to 3.03x 10? dynes/cm? from the surface up to the depth of 
2898 Km., below which is assumed the existence of fluid layer, mainly from the fact that S-waves 
are not propagated through it. 


With the above data we proceed to calculate the stresses and displacement on the surface 
orsthesearthi gr eat zi =. 


at z= 0, rigidity u, = 0.63% 107 gynes/cuPr, 
at z=h=289 km, rigidity = u, = 3.03x 10%? dynes/cm?., 


il a — 0.08283 x 10% dynes/cm®. 


whence, we get, 


& _7.605x1®cm and = al — 36.58x10°cm. ... . . (83). 


Kı 


These large values forming the arguments of K and I functions, we obtain the expressions for 
stresses and displacement from Eqs (18), (18.1) and (18.2) on using the relations (30) in the fol- 
lowing forms 


Go) o=-ka/A,AadIl@ntanhAhh ....2.2.2.22:.0684M). 
0 
Hero =ka/IAadJlar)tanınAhadA........%. Al), 
0 
ka Rs 1 
(VD), = — F ;: Ja) Jı(Ar)tannAh)dä ...... . 842). 


Remembering that for large values of h (2898 km), tanh (A h) tends to unity and we thus find 
the above expressions all agreeing exactly with those given by Yi-Yuan Yu [3] for torsion 
of a semi-infinite body. So the numerical calculations and results are exactly the same as given 
by Yi Yuan Yu [3]. Hence we conclude that the variation of rigidity within the earth is 
practically without any effect on the stresses and displacement on its surface, 
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Poiseuille and Couette flow of non-Newtonian fluids 
By James Serrin 


Vor einigen Jahren fand Rivlin einige exakte Lösungen der Bewegungsgleichungen einer nicht- 


Newtonschen Füssigkeit, die den wohlbekannten Lösungen der Navier-Stokesschen Gleichungen ent- 
sprechen, welche von Poiseuille und Couette angegeben wurden. In der vorliegenden Arbeit wird die 
physikalische Bedeutung der Rivlinschen Lösungen geklärt und eine Methode zur Bestimmung des Koeffi- 
zienten der Querviskosität vorgeschlagen. Ferner wird gezeigt, wie sich gewisse Experimente von Merrington 
und Weissenberg auf Grund der Resultate dieser Arbeit erklären lassen. 


Some years ago Rivlin discovered several exact solutions of the equations of motion of a non-N ewtonian 
fluid, corresponding to the well-known Poisewuille and Cowette solutions ofthe Navier-Stokes equations. 
In this paper the physical nature of Rivlin’s solutions is clarified, and a method is proposed for determining 
the value of the cross-viscosity coefficient. In addition, it is shown how certain experiments of Merrington 
and Weissenberg can be explained by means of the results of the paper. 


Il y a quelques annees que Rivlin a trouve quelques solutions exactes des equations de mouvement pour 
un liquide non- Newton, solutions qui correspondent aux solutions bien connues des equations de Navier- 
Stokes indiquees par Poiseuille et Couette. Dans l’article present la signification physique de la solution 
de Rivlin est Eclairee, et une methode pour la determination du coefficient de la viscosite transversale est 
proposee. En outre, il est demontre comment certains excperiments de Merrington et Weissenberg 
peuvent etre expligues sur la base des resultats de cet article. 


HeckKoAbBRoO AeT TOMy Haaan PHBJIMH HAalMeJI HeKOTOPLIe TOYHbIe PelleHnA ypaBHeHnü 
ABEFREHMA HEHIOTOHOBCKOH FKUJIKOCTU, COOTBETCBYFOIIHE XOPOIMO U3BECTHBIM PeIUeHuAMIIyaceına 
u KyarTa aıa ypasHenna Hapbe-Crorca. B HacTonamei pa6oTe BhIACHAeTcH Pusmyeckuü 
xapakrep pemenHnüi PusJmna m TIpennoraetcaH MeToN AA ONPAelIeHNuA 3HayeHuAd Koabeu- 
IMeHTa MOoMepeyHoü BA3KOCTU. KpoMe TOTO IOKA3bIBAeTCA, KAK HEKOTOPBIE ONEITBI MeppHuHr- 
ToHa u Beliccen6epra MOTyT ÖbITb OÖACHEHBI IIPH IIOMOIIM Pe3yJIbTAaTOB HAaCTOAIMeH Pa6oTkI. 


The theory of non-Newtonian fluids is based on a stress-deformation relation of the form 
Ta lEERDH DE egetrreee tl), 


where T is the stress tensor, D is the rate-of-deformation tensor, and a, ß, y are certain scalar 
coefficients!). We shall be concerned here with incompressible fluids, so that «& is identified as 


the negative of the dynamic pressure, x = — p. The coefficients ß and y are assumed to be known 
functions of the scalar invariants of D and of the thermodynamic state of the fluid. That is 
za, 11:7), ERSTE 


where II, III denote respectively the sum of the principal minors of D, and the determinant of D, 
and T is the fluid temperature. 
The fluid motion is governed by the familiar Cauchy law 


en =of+aivT a ee (2); 
in which v is the velocity vector and f is the gravitational force, f=—-yk. Some years ago 


Rivlin ([2], [3]) discovered several exact solutions of equations (1), (2), corresponding to cer- 
tain well-known solutions of the Navier-Stokes equations. We are interested here in two of 
these exact solutions, namely those corresponding to Poiseuille and Couette flow: our pur- 
pose is to clarify the physical nature of these solutions, and especially to point out their similar- 
ities and differences with the classical solutions. Some steps in this direction have already been 
taken by C. Truesdell?), but the present analysis seems to be more complete. 


1) The reader unfamiliar with the subject of non-Newtonian fluids will find adequate references and 
background material in the extensive monograph of C. Truesdell [5], or in [4], $$ 58—60, 65. 
2) [5], $72. This work is, in fact, the inspiration for the present paper. 
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Finally, we observe that equation (12) offers an excellent opportunity to deter- 
mine an experimental value for the coefficient y in a second order theory. This 
statement will be expanded in the work below. 


1. Poiseuille flow 


In eylindrical polar coordinates (r, ©, z) we seek, following Rivlin, a solution of equations 
(1) and (2) of the form 5 
„-%=0, „=w=wm), OSTSa...2.» 02... 


This corresponds to steady laminar flow down a straight pipe of diameter 2 a. For such a motion 
the deformation tensor D takes the form 


1 0) 0 w > 
D-z|0 0-0], 
m) 0 
and therefore 
0 0 1 1 0 0 
f $ 1 2 
dc 0 0 rd 0 u18 
1 0 0 0 0 i! 


1 
where $ = ß > w'2, 0; r) = ß(r, T), y=y(r, T). Since the motion (3) has zero accele- 
ration, Cauchys law reduces simply to of + divT = 0. Moreover, using the techniques of 
tensor analysis?), one finds that the covariant component of div T is given by 


ce 1720 or 
2 +7 k 2 | 
(div T), = € = Fre) 73 — it ae EEE EI 


where x, = r, x, = 9, 2, = z. The equations of motion thus become 


opt = *) ai | 

or r Fra hl Ho 5 Bw) (5) 

op en ) o/1 e ur 
er Te ee + (gro) ER 


in which the ©-dependence of all quantities is taken zero. Equations (5) can be integrated pro- 
vided that $ and y do not depend on the temperature, or if the temperature is assumed to be a 
function of r alone. Indeed in this case one finds the following unique solution of (5), 


p=(lCc-eNer iQ, Pwi=lr, Ze Se 
where C is an arbitrary constant, and 


cm A | yr 
Ir) = ap 4 2 dr 
The second equation of (6) is a differential equation for the velocity w(r), to be solved subject 
to the boundary condition w(a) = 0, (i. e., adherence at the pipe wall). Due to the fact that ß 
will in general involve w’, the solution of this differential equation may be difficult to obtain 
explicitly. 


To complete the solution, we observe that the stress components are given by 


TE 1 ok \ 2 Pr — 1 
IT = 7% p- 1 ywa=lodg ORET Te z=-Cr, 
and that the total mass flux M is 
M=—/[fowrddO= —no[wdr)=no c/rp- er 
1) (0) ; 


where in the last equality we have used the boundary condition w(a) = 0. Notice, in particular 
that the normal stress rr on the pipe wall varies linearly along the length of the pipe exactly as 
in the classical solution. 


°) The reader is assumed to have an e i i ; ; 
Ah ) £ ave an elementary knowledge of this subject in order to derive 


J 
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P 


Quadratic dependence. In the special case where the components of T depend at 
most quadratically on the components of D, it can be seen that ß and y must depend, at most, 


on the fluid temperature. Supposing in fact that ß and y are constant, integration of (6) yields 
the familiar formula 


C 
= eo (N U N 
while from (7) follows 
20€ 
8 u 
Since y does not appear in these formulas, it is clear that standard viscocity or velocity measure- 
ments will not indicate the presence or absence of non-linear viscosity in a fluid. 

The effect of y is felt, however, in anon-uniform distribution of pressure across 
asection of the pipe. To investigate this effect in detail, we may suppose the fluid to issue 
from the pipe into an atmosphere at pressure p,, the latter exerting a force on the output cross- 
section in the amount z a? p,. Since ß and y are constant, the stress components are given by 


fr =22=(eg— ()z— 2yD?r? — Const. ER, (10), 
where D = (/8 u, and the following force balance holds at the exit section (taken to be z = 0), 
N 


M = 


RN SR I Sl RO 


a 
rap, =—- 3 [22 dr) =nyDeat + nat Const. 
N) 


Using this equation to eliminate the Const. from (10), we obtain 


ge 2=eohelr Di —2N pe ee in 

Let P be the force per unit area which the fluid exerts on the pipe walls; that is, P= —ır 
evaluated at r= a. In virtue of (9) and (11) we have then 
are+) 

—p = |\—) |— ee ende ee 

P—p (<) pt’ (12), 


where pü = p, — 0 g zand I" = W/r ao = volume of flow per second per unit cross section area. 
Formula (12) may serve as the theoretical explanation of an interesting experimental 


phenomenon discovered by Merrington [1], namely the tendency for a fluid stream to swell at 
the exit section of a viscometer. Indeed, according to (12) we have 


EEE Dar 


at the exit section, indicating that excess pressure in the fluid can be caused simply by a posi- 
tive coefficient y. Formula (13) implies, moreover, that the swelling should be emphasized by 
high flux and small pipe radius, a fact actually observed by Merrington. 

An experimental verification of (12) could easily be secured by measuring the variation of 
P—-p,along the length of the pipe, for various values of /'. By plotting (a/I') (P — p$) 
against the pipe length z, one should then obtain a series of straight lines, each intercepting the 
axis z = 0 at the fixed height y; moreover, from the slope of these lines one can read off the 
coefficient of viscosity u. The internal consistency of the theory is immediately decidable from 
the range of values of u and y which are experimentally obtained for various values of /. To 
the author’s knowledge, the required experiments have not yet been performed. 


2. Continuation 


The conclusion that the fluid stream should tend to swell at the exit section of a pipe can 
also be derived without the assumption that 8 and y are constant. Indeed, in the general case 
let us write 


tT=22=(eg—(C) (CP [y 8° r dr — Const, een. 0 LOy® 
Then the following force balance holds at the exit section, 
rap = — nf dr) = —n (Ca [yPß?r’dr +na:Const. 
Ö {) 


where a straightforward integration by parts was used in deriving the final equality. Using this 
equation to eliminate the Const. in (10) leads to the formula 
{2 d 


Pen rt 44 a a | De. 


#7 = 232>= (ag-'C)2 
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and finally 


a 


GY r? f 
P-n=)| Ga... r022 22. 


0 
at the exit section. Thus whether or not ß and y are constants, the fluid stream 
will tend to swell ify>0. 
Now let /’ denote as before the volume of flow per second per unit cross section area. Then 


da 
TEN terre, 
0 


and so by the Cauchy inequality (assuming y > 0) 


[73 


a 4 
2 [ p3 3 ”2 5 3 
es =) u a (5 dr zur | Fr. 
a Miny J PB 
Combining this inequality with (13)’ yields the interesting result 
P—p, 2 (Tl? Miny. ns Ders Dee 


Inequality (13) not only exhibits a remarkable similarity to formula (13), but also shows as 
before that the swelling should be emphasized by high flux and small pipe radius. 


0 


3. Couette flow 
Again using cylindrical polar coordinates, we seek a solution of (1), (2) of the form 


Dr m 0, ta =IT@lr). - . wu a 
In this case 


1 0 ro 0 1 ro? 0 0 
D-- To 0 0 3 2, 0 To? 0 s 
= 0 0 0 0 0 0 


s 1 
(using physical components), and $ = Bl A 72 073,06; r) = ff(r, T), y=y(r, T). With ihe 
help of (4) the equations of motion are easily reduced to the following system, 


6 AR re ) 


Tr ör ' Or i 
“u Op Le arts = 
re eh 


VE een 12 


Assuming that and y do not involve T explieitly, or else that T is a function of r alone, these 
equations have the integral 
P=-09z2+ffr), Po =CG. nu. 
where C is an arbitrary constant and 
Ir) = 39] ri ro ro®:dr. 
The only stress component whose value is affected by y is zz, as the reader can easily check. 
Quadratic dependence. Setting ö — constant, and y = constant, as before, integration 
of (15) yields the velocity distribution 
v=ro=Art+Br, (A = G/aP) sera er 
‚exactly as in the classical Couette flow. A difference arises in the stress component 72, which 
now takes the form 
22 =—p=o092—Ayr—of(A+ BRW%r dr. 
To investigate the effect of y on the motion, we suppose the flow to take place between two 
concentrie rotating cylinders, with the upper fluid surface being left open to the atmo- 
sphere. We shall be interested in the shape of this free surface. Now Couette flow involves a 


3 
r 
1 
- 

; 
E 
| 

| 


En Eu 
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non-vanishing stress component T,o, so that (16) can only approximate the actual flow under 
consideration; nevertheless, as long as the free surface remains relatively horizontal (i. e. at rela- 
tively low speeds of rotation) the discrepancy should not be serious. Assuming, then, that the free 


_ surface remains nearly horizontal, its locus is given by 22 = constant, and its slope satisfies 


de _ Ar 
dr gr 


4 | B 
1 a — —- =—?f. 
ja+ 9m 2, 0=4r 
We shall consider two cases in detail: () The inner cylinder (r = r,) remains at rest, and (II) the 
outer cylinder (r = r,) remains at rest. The constants A and B for these cases are easily deter- 
mined, and in particular we find that 


= = (z/7,)° in Case I, 
— (r/r,)? in Case II. 


When y = 0, the slope in Case I varies from 0 at the inner cylinder to a positive value at the 
outer cylinder, while in Case II it varies from a positive value at the inner cylinder to 0 at the 
outer cylinder. When y > 0 these slopes are obviously decreased somewhat. This leads to the 


<b 


following general conclusions: 


Y 


Ne 


1a 


Fig.1. Typical fluid cross seetions in Couetteflow of a non-Newtonian fluid 
1a. Theinner eylinder atrest. 1b. The outer eylinder at rest, and r, > 2 VrIe. le. The outer eylinder at rest, and r, <2 VyIe- 


Case‘l. When »>0O the fluid will.tend to climb the inner eylinder, and 
will generally also climb the outer cylinder, unless y is very large (Fig. 1a). 

Case II. The fluid will tend to fall slightly at the outer cylinder. At the 
inner cylinder the fluid will rise or fall according as the quantity 


4 ee INH. la) 


is greater or less than zero (Fig. 1b, c). In the important case (r, < r,) this leads to the 
eriterion: the fluid will tend to climb the inner cylinder if 


Te 2 lo) we el 


and to fall away from the inner cylinder otherwise. 

The above conclusions are derived under several approximating assumptions, the most 
important being that the free surface should not greatly vary from the horizontal. This qualifi- 
cation should alway be borne in mind, though in all probability the general conclusions remain 
valid in a considerably wider range of circumstances. ‘The conclusions above are in general 
agreement with certain experiments of Weissenberg [6], and very probably furnish the correct 
theoretical explanation of his results. So far there has been no experimental check of the in- 
equality (17), though this would certainly supply a strong test of the preceding considerations. 
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The effect of square and triangular notches with fillets 
on the transverse flexure of semi-infinite plates 
By S. Shioya!) 


Die Aufgabe der Abschätzung des Binflusses einer quadratischen oder dreieckigen Kerbe auf die Span- 
nung in un dünnen, ee Halbebene wird auf Grund der Kirchhoffschen Biegetheorie 
und der komplexen Darstellungsmethode von Muskhelishvili behandelt. Zur Bestimmung der Koeffi- 
zienten, die in der Durchbiegungsfunktion der Platte enthalten sind, wird die Störungsrechnung angewandt. 
Numerische Ergebnisse für einige spezielle Fälle werden mitgeteilt. 

The problem of estimating the effect of a square or triangular notch on the stresses in a semi-infinite 
thin Bi eh to bed nu is treated on the basis of Kirchhoff’s theory of bending and by 
using the method of complex representation due to Muskhelishvili. T’he perturbation method is used 
to determine the coefficients contained in the bending function of the plate: For some special cases numerical 
results are given. 

Base sur la theorie de flexion selon Kirchhoff et sur la methode complexe de representation selon 
Muskhelishvili on traite l’&valuation de Vinfluence d’un cran quadratique sur la tension dans un demi- 
plan fin avec effort.de flexion. Afin de powvoir determiner les coefficients contenus dans la fonction de 
flexion de la plaque, on applique la calculation des erreurs. Des resultats mumeriques pour quelques cas 
speciaux sont communiques. 


3anaya OUeHKH BJIIMAHUA KBANpaATmyecKof HaceyYkKu Ha HAaIPAsKEeHUA B TOHKOU, TONBEPTae- 
Mo H3THÖe IIOJLYILIOCKOCTH, pellaeTcH Ha OCHOBAHUH Teopum uarudanun Kupxroda Hu KoM- 
ILIEKCHOTO MeTona IpencTapieuun Mycxesnmmsusm. Isa onpenerenun KoehhHLMeHTOB, 
comep;kamuxcaH B PYHRUHM W3INÖA INIACTUHKU, IPHMEHHETCH BUYHCHIEHUe BO3MYINeHHÜ. 
Ina HeKOTOPbIX YACTHbIX CAIyYaeB COOÖIMAIOTCH YUCJIEHHBIE PEe3yJIbTAaTbI. 


1. Introduetion 


The problem of the stress distributions that arise when a notched isotropie plate is bent 
by transverse couples has already received some theoretical attention in certain special cases. 
Lee [1]?) has investigated the influence of two symmetrically disposed hyperbolic notches on 
the transverse flexure of thin plates. Neuber [2] also has discussed similar problems. Tamate 
and Shioya have given the theoretical solution for a semi-infinite plate with a semi-eircular 
notch under plain bending [3] and for an infinite strip containing two symmetrically placed 
semi-cireular notches and subjected to plain bending [4]. 

The present paper gives a theoretical solution for a semi-infinite plate notched by a single 
square or triangular notch with fillets and transmitting a constant bending moment. The analysis 
is developed on the basis of the Poisson-Kirchhoff theory of thin plates [5] and by means of 
Muskhelishvili’s complex variable method [6]. A method of perturbation is adopted for 
the determination of parametric coefficients involved in the solution. Numerical examples also 
are worked out in some detail, in order to clarify the effect of the notch. 


2. Preliminary Results 
When the middle plane of an isotropie homogeneous thin plate free from lateral load lies 
originally in the complex plane of z = x + iy and the plate is assumed subject to the limitation 
of the Poisson-Kirchhoff theory of thin plates, its deflection w normal to the z-plane is 
governed by the biharmonic equation. As is well known, the solution of the biharmonic equation 
may be expressed in terms of two analytic functions 9,(2) and y,(z) of the complex variable z 


w=Re [292 +u@)]), HF Ertl 
Overbars are used to indicate conjugate complex quantities, and ‚Re‘ denotes „the real part 
of, 
With the aid of mapping function 


z='ulß), en ee 
curvilinear co-ordinates (o, ©) are introduced into the z-plane and we now have the deflection 
of the plate given by 
w— Re |[wo(£) Ar) Fr] re ee (3), 
where 
| rI)= ne Ad, 2A =SYO A 
anc 


ve) = yılold)] = yıl) . 


) Postgraduate Course of Mechanical Engineering, Tohoku University. 
) Numbers in square brackets refer to the Bibliography at the end ofthe paper. 
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The bending moments M,, Mo, twisting moment.M,., and shearing forces Q,, Oo, all 
per unit of arc length with reference to the new co-ordinates (o, ©) are expressible by the formulas: 
M, + Mo=—4D(1 +») Re [90], 
M—M,—2iMas =2DU-„2lEHA) +FO ONE. - - 
9. — 10 = —4 DL PHKe lv 
in which 
PH) = PAWd, FO =YAIW) 
and D = Eh?/12 (1 —»?) is the flexural rigidity of the plate; E, v and A being, respectively, Young’s 
modulus, Poisson’s ratio and the thickness of the plate. 
For a free boundary of a plate, the boundary condition is given by 


— (JE) +) TO SI d+TE::.. 2... (8), 
e=(1—»)/8-+») 


and C, is a real constant and C, a complex one. When the region occupied by the middle plane 
of the plate is simply connected, these constants may be taken to be zero. 

With these formulas we shall seek the solution for the case where the semi-infinite plate, 
having a square or triangular notch with fillets, is subjected to plain bending in the direction 
of the straight boundary. In this case, the straight boundary and the rim of the notch are assumed 
to be stress-free. 


where 


3. Solution of Problem 


Let the middle plane of the plate occupy the right half of the z-plane, with the exception 
of the interior of a notch whose center is at origin (Fig. 1, Fig. 2). 


Fig. 1. Shape of Triangular Notch Fig.2. Shape of Square Notch 
Bv means of the mapping function 
* met 
et) ET. er <Lth) 


and by varying the value of the real constant R and m, the straight boundary of the plate and 
the rim of a square or triangular notch with fillets will be mapped onto the line Oz de n/2 and 
o= 1, respectively. Namely, when 0O<m< 1lor —1l1<m< 0, the corresponding notch 


resembles in shape a square or triangle with rounded corners, respectively. 
99 


302 8. Shioya, The effect of square and triangular notches with fillets .... 


m 


Then the boundary conditions of the present problem are expressed as follows: 


() 9=+nß 


dee rer 2), 


od) 
— (1/e) KE) + he sd re ee ra 


The load applied at the infinite part of the plate is M, = constant = M,, with M;.;=M„=0, 
which may be easily transformed into 


() e=1 


ü) > © s 
M,+Me=M,, | 


Me — M,— 2i M,e = M, e'°] 
To satisfy these conditions, let us consider the functions 


le) = 9) + "pC) al 0 ee (10), 
Ye) = pl) + "y) 


where !o(£) and !y(£) correspond to the deflection of the plate, without the notch, subjected 
to plain bending and can be written as 


In) = IM RIED ACH Z 5) 


Yye&) = IM, Re Da—mlc+7 .) a 


Uo(d) and !!y(£) are the respective parts of the functions which are analytic in the entire region 
of the semi-infinite plate surrounding the notch and therefore include only terms of negative 
powers of £. Taking into account the requirement of the symmetry of the problem, we may 
write Igp(2) and !y(£) in the form of series: 


I yes ; den —(@n+2 —(2n+ 
pL) = ab Tre ES (en+2 4 ( 


er M,R c(£) en ae, > j 
ud) = El PAR RE de De = in wre 
> _ an} r—(2n+2 ei 22 2n 
a 


where qg4s Dans Con and d,, are real constants. 
Substituting Equations (11) and (12) into Equation (7) and equating to zero the coefficients 
of various powers of go, we obtain 


On = — An; dyn au: Dan 
€ € 


In a similar way, Equation (8) gives 


1 2 s 
abe ner fiene Bela Be end 


4 6 


x a 
- 2 5 {{e2Ri® 4 me-2nie) _. (e-En+9i0 m ent 9io)) 
GE = n-+2) 


oo 


un 2 Ä Er 5 en ji ((e@n- 1,3822; me BR STE) + (ern +8 m e2n+3)ie)) 


n=(0 


2. a = — e(2n-+2)io + z (dan- 2)i® = a) 


n=0 


a Sn bi [erde + e@n+1)i0 er 5 (e® en Bio a en a . (14) 
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from which it follows that 


ml +3)eo © + (1—m)cos3 © + m (1m) 0os5 ol 


+ I |a en 3 0411- ana) RO +men+20 


— e(2n-+2) 
m 1+m m 
Hat omg en + 4) a +2 yen—3 [0] 


er nen) O+mEn+N9 


en ja ai ve, cos (2n +3) 0 I) ” 
(15). 


al + m) 1-3) o+a + m) sin 3 9+74 + m) sin 5 e) 


ee ang RO + in @n+ 20 


3 nf (2n+4) Ola +3 ]3% (2n—3) © 


en) m @n-Do+m@n+nE 


202: +1) 
m 1+m 5 a 
= un sin(2n + 3) 0 = 


Besides, sines and cosines of even multiples of © can be expanded in Fourier series of sines 
or cosines, respectively, of odd multiples of © valid between O = — rz/2 and 9 = n/2: 


cos2n® ee c03.(2:5-- 1):@,, | 
L er Zzn 
sinn En Baer | 


Substituting Equations (16) into Equations (15) and then equating to zero the coefficients of 


4 cos (2s +1) ©, sin (2s + 1) ©, we get 
; 
1 m 1 m 1 m 
„1m (+ 3) ja m) E an li +- 3)+ 30m} + Zr 
| u Mm er Be > m, Im)  M'@m 
- Al — e@n+DV@n—I)@n+l) aE]3 e(2n+2)|(2n+3)(2n+5) 
| — 1)" agu 4 A u ne 
Br on een 
m 1 1 m m l—n 
bad + m + m), 4 Bi ee "ee ti a 


1 AN i+m —Y"r&s+ 1) an 
+ ( Ur@s re) Dash Dana — len e@n-+: a On—2s—1)@n+2s+i) 


=0 


4 m 1l+m (— rt? 2 s +1) an 
Abt cenray 


ee / e(2n +2) et | 
NH Qs+ 1a DEAD ug 
Er en 0 ann 
61) 


22% 
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ee; "(1 3)+ "(4m 13) +2,14 m) ee 


4. & 1m (— 1)" 2n 
Deo Da 

4 (m in Je pen den = (— 1)" (An + 2) Gym 
ee ee Nee) 

4 S&m(— 1" 2@n—2)a _ 
a2 Leone 


m l+m (18), 


1 m 
ut tm + tm A { Fr) en 


im 4X im MH ana 
er u (st jr" tz I "elent a s—l)@n+2s+1) 


a N LE (— 1rt3 (2n + 4) Qyn 

en re 

ER (— Drts(2n + 2) Qgn a IPA au 
m @n-2sH)@n+2sH3) n,3(2n—2s—3)(2n+2s—1) 


9) 
where ö, , is Kronecker’s delta. 


To solve these equations, we may apply a method of perturbation, in which m and e are 
the perturbation parameter. Let us assume the expansions 


( — > SI Na er m, => 2, DENE ME. ie er (19), 


r=0 9=2 


where Na and "5% are independent of m and e. Substituting these into Equations (17) and 
(18), and then equating the coefficients of the same powers of m and e on both sides, we get 


= 014 in uEN 
2 U UN 76 rl I) 2.E 1)* DIT Ma9 


3 
l 
[=) 


= B 2 
+, ap (— 1)% Br Br (a2) B 
(20), 


n 2 N, Re E 
2 (— Ijrtetl a2 (al) = 16 I, Al — d2,,) 2 ( — JJkte+1 2 2: Van 


+) Z- Ver @Kk— 25 +30 Ga) 


[e,0] 


3 n r p 7 hu 1 
3 1, map) = gg Our — dan) dan + (1—8,,) | 5 Ir 


1 
(r —1)h(r—2) er —1 —1 
—Vptp +[ en 


wo 


n=0 


1 | 1 R 1 1 
AB : DE) + = (r— a) AR (= 5 DH + — aka) 
13 1)® I oo NaPp—D + (1 —6, ;) iS 1)% 1 0 (Naw—2) 
ge Pak "2% Op,3, zn ) ) Pr, Ya, 


+2 Pot, R— 1) (ak + — 80 Val) 


s (9,0 1 g— S ‚1 2 
eb Bas 1)* (2 Ö0, — u, ) ( ale re (1 — dp) 2. 1)* 5 “2, 2] x 


# 
. 
r 
: 
& f 
F 
Br 
= 
} 
4 
f 
| 
S 
; 
e 
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= NS 8 7 7 
2 1jr+s+1 023 (Nato) z rs et tönt dr (21), 
P1A 

+ U dere + Barn) + 15 daran e=1) 
7 ur 1 5 a 

— 48 2r,s SS 16 (2r ne Ög,p > (1 0) 8 15 at Ir 2) Ö,s 
ee 885 a 2s+1 | 

_— EHEN (1 --6,,) — 303) 3). DPD A DR —D + DH 6, , 

1 A 1 = @ 
£r Ne a more il 2. 2, Unna: 2% Sure 


+(—6,,,) 2 (— 1je+>+1 (2k— 25 + 3) 2: Qaß—» 
+ = (— 1+s+1f@n +7) (An — 3) — 4 (s +2) (s — 1)} 02642 «262 1a) 
ar 2, HH MSOCKLF2SE HN 2er N) <-S@k2sNnN a 


we (di —b;, DPA (— N {(2 k 4 1) (2 k 48 5) en (52 u 1)} Rn N) | 


sZ1) 
eat L 3 Ns 4 s 122 
He a 4%0s +7 9: )0&,p +(1—6,,9 63,2) =6 s+ D = en 
A [6.0] 
ET IE pn ur ach Pa (22), 
1 1 3 35 7 2s+1 
‘ Ip) = gs Ö1,r Ö1,» Ar 120 + dns Fr 1202 Zr une + on = (der —4),8 Fr Ö2r—3),5) 
2s+1 7@Qs+1 2s-+1 
= 3 Öar—1,s + - L 12 si Öor,s + 1 Berne)! Op da %,) 


1 2s+1 4 @ 
x E a ner Opa + 375 en + = 2s+ u jjerenl a nal 


4 » . 4 @ 
—— @s+ 1) yet ara nad —— (25 + 1) > (— 1jrtsHt 92649 0 -Daß) 


+ 2 s+1) 5 a Ve | I) EN er A (23), 
n=0 

where 
28 1 28 1 
%, = ent Non > 
2" (2n—2s+1)(A2n+2s+3) 2» (2n+2s—1l)2n+2s-+|) 
628 ee a EL 3 us Ben B 28 — - = er pi e 
= (an —2s—3)(2?n+2s+1)(2n+2s+7)’ ?n (2n+2s+1)2n+2s+7)’ 4) 
Bu- iron Bam oe 5 02, 
2n (2n+3)(2n-+5) ?n (2n+2s+1)(2n+2s+3) @n+2s+5) 

== 1 28, — = E = 1 nn 
Mn öön-1Oonfb)’ Man @n—2s—1) @n—2s+5) @n+2s+5) e=1) 


With these relations, the coefficients Pa) and ®b%) in the Series (19) can be determined suc- 
cessively frrr=0,1,2,..., andp= 1,2,3,..., respectively. In particular, the coefficients 
al) and Mb), corresponding to the solution for the semi-infinite plate, notched by a single 
semi-circular notch, was already obtained by the authors [3]. 
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Fig. 4. Bending Couple M, on the Minimum Section. (vr = 0.3) 


10 20 30 40 
ro — 


Fig. 3. Bending Couple M, on the Minimum Section. (» = 0.3) 


Be 


Fig. 5. Bending Couple Mg at the Rim of the Notch,. (» = 0.3) 


Fig. 6. Bending Couple M,, on the Straight Edge. (» = 0.3) 


Fig. 7. Stress Concentration Factor ky versus e 
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It may be noted that the coefficients of a@) are invariant for all values of r and p in the 
infinite sets of Equations (20) and (21). Hence, if we solve these equations once, the solutions 
in successive cycles are obtained only by the calculations about the constant terms. 

The foregoing solution is valid in so far as the Series (19) is convergent. For brevity, no 
proof of convergence will be pursued here, but it will be illustrated by numerical examples which 
follow. 


4. Numerical Examples 
The foregoing solution is worked out in some detail. The calculations of Nam, nd are 


carried out until r=7andp=5, but their values are not shown here for brevity. By using 
these values, it is now rather straight forward to calculate the stress resultants in the neigh- 
bourhood of the notch. 


D/Xo 


Fig. 8. Stress Concentration Factor ky, versus r/Xo 


(1) Triangular Notch With Fillets: For example, we may take the values of m as 1/3, 
1/5 and 1/10. Then we have the shapes of the notch as shown in Fig. 1. On working out the 
values of M, and Q, — (0M,o/0s) on the rim of the notch for m = 1/3 and v» = 0.3, where s 
denotes the arc length of the notch, the greatest value found is less than 0.04 M, (M, for © = 0). 


Therefore, the values of parametric coefficients obtained are thought to be sufficient to get 


the stress resultants with satisfactory approximation. When m is smaller, the approximation 
is more satisfactory. 

The values of (i) M, and M, along the line of symmetry, Ox, (ii) M« on the edge of the 
notch, (iii) M, on the straight edge are calculated for these values of m and for v = 0.3 and 
shown graphically in Figs. 3, 4,5 and 6, respectively, where x, is the depth of the notch, shown 
in Fig. 1, and o is the distance from the origin, O, to given point in the plate. 

The maximum bending couple occurs at the apex of the notch. Since M, = 0 on the same 
edge, the maximum couple is given by 


Maa/M, = 1-21 — g/l — m} 3% MY Bee er 2) 
p=1r=0 
where 


[0,0] 
— )' Na) 1 
U de 


By use of this equation, the stress concentration factor k, = Mynas/M,1s calculated and plotted 
versus e in Fig. 7. It is interesting that the values of k, changes almost linearly with e. For 
the limiting case of m = 0, the results agree with those for the semi-eircular notch [3]. The 
stress concentration factor k, is also plotted against r/x, in Fig. 8, where x, is the depth of the 
notch and r is the radius of curvature of the notch at its apex. 

(2) Square Notch With Fillets: Let us take the values of m as —1/3, —1/5 and —1/10. 
The corresponding shapes of the notch are shown in Fig. 2. The distribution of stress resultants 
near the notch for these values of m and for » = 0.3 are shown in Fig. 9. When |m| is small 
in this case, the maximum value of the bending couple, occurring at the apex of the notch, 
decreases gradually with increasing |m|. This may be due to the decrease of the curvature at 
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x. However, when |m| increases further, the point of the maximum bending couple 
ae a shift towards © . 2 and the value of the maximum couple also begins to nn 
The stress concentration factor k, for this case is plotted against e ın Fig. 10. Fig. 
shows the relation between k, and r,/x,, where x, is the depth of the notch and r, is the radius 
of curvature of the notch at © = 7/4. 


. EZ = 


1,0 


0,15 0,20 025 030 
E—— 174 % 
Fig. 10. Stress Concentration Factor k, Versus e Fig. 11. Stress Concentration Factor k, Versus r./Xo 
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Centre of flexure of a beam of orthotropic material having 
' a section bounded by an ellipse and its major axis 
; By A. Chakravorti 


Diskutiert wird das Problem, wie man das Biegungszentrum eines aus orthotropem Material bestehenden 
Stabes bestimmt, dessen Querschnitt von einem Ellipsenbogen und der größeren Ellipsenachse begrenzt 
wird. Es wird angenommen, daß die Last parallel zur geraden Kante wirkt. _ f 


In this paper the problem of finding the centre of flexure of a beam of orthotropic material having a 
section bounded by an ellipse and its major axis is discussed. Here the load is supposed to act parallel to 
the straight edge. 


On discute le probleme comment il est possible de determiner le centre de flexion d’un baton consistant 
de materiel orthotrope dont la section transversale est limitee par une courbe elliptique et par Vaxe elliptique 
la.plus longue. On suppose que la charge agisse parallelement ü Varete droite. 


PaccyztnaetTca BOUPoc ONpeneJIeHusn IIeHTPa u3Tu6a CTeP>KHA 3 OPTOTPOHHOTO MaTepuara, 
ceyYeHHE KOTOPOTO OTPaHnyeHo Ayroü u 6oNBIUefi OCbIo HJnsıHrca. Ilpennosnaraetca, YTO IuHuA 
MelCcTBUA HATPY3KH MapaslıejIbHa IIPAMOMy P&öpy. 


Introduetion 


Though the position of centre of flexure for such simple cases as that of an isotropic semi- 
circular section with a load parallel to the straight edge has already been determined (Sokol- 
nikoff [1]) the corresponding cases for non-isotropic material have not yet been discussed by 
any previous investigator. By appropriate modification of the method suggested by Love [2] 
for solving flexure problems of non-isotropic material, the problem stated in the title has been 
solved. 


Nomenelature 
24,2 — Rectangular coordinates. 
l,m, n — direction cosines of outward drawn normal. 
0%, 0, 0; = normal components of stress parallel to x, y and z axes. 
T — shearing stress. 
T(&, y) — stress function in rectangular co-ordinates. 
Toy Tyz Tzz — Shearing stress components in rectangular co-ordinates. 
u, D, w = components of displacement. 
Yyz Yan Yay — Shearing strain components in rectangular co-ordinates. 
O2» Oy» ©; — elongation along rectangular axes. 


Let a beam of uniform cross section with a straight edge have one end z = 0 fixed and 
the other end z = !loaded by some distribution of forces that is statically equivalent to a single 
force W acting parallel to the straight boundary. At the end z = 0 axis of x is taken along 
straight boundary, axis of y perpendicular to the straight edge and z axis along central line of 
the beam. The origin is at the mid-point of straight edge. 


Here x and y axes are parallel to principal axes of inertia at the centroid. 
Following the semi-inverse method of St. Venant, we assume as in case of isotropic beam 


W 
%=y=%y=0 and re DEE DE EL RN 


where I = ff x dx dy integration being taken over the entire cross section. 
R 


From equations of equilibrium we get 


u I Eee, 
and 
ZERNNE AT ARE N EUR IE VRR a anen: 


From (1.2) we see that r,,, ,, are functions of x, y only. It will be seen that (1.3) is satisfied 
if we assume 


We, a nel, 


oT 
Tr, = ay + S(Y) 9] 


But for non-isotropic bodies having 3 orthogonal planes of symmetry 


Tyz = Lyy:» Tyx = My; Toy =Nyy 
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and 


In this particular case we get 
= 610, ,=F9%, ,.=G%> = Ly:> %s=MY.- 
From compatibility equations we have 


Jene E 
Mare ar L v: ray\m ”- u DE 1 a ee: ( % 
u u E ls \ 16 
2A le ge) DE (1.6), 
ll + au) 
Uzeiz Er L ’r: + M ee ee (1.7) 
From (1.6) we get 
| WE = ee en | 
End ine > 
a a a 
er tm TMmi: 
Relation (1.7) is identically satisfied. 
From (1.5) and (1.6) we get 
102702 1@7 10S Ww 
= eo — Y—y)—2%9 - -- ++. 1.8), 
L 9 "M a? M 2y 2F, 1 (y — Yo) %o (1.8) 


where y, is y-coordinate of centre of inertia of cross section. Local twist at point (x, y) of cross 
section is 


io 1 a en. nr ep ER Braun Eau: 
a: 2. a) 2a te yt L%& M% 
1 ars 1er a w 
her oe "M ap ma TU-W+n. 


We define centre of flexure as the position of load point for which the mean local twist is zero. 
Then we have &, = 0 in (1.8). 

Now let (@,y, I) be the coordinates of centre of flexure. Since in the present case W is 
directed parallel to a principal axis (x axis), x coordinate of centre of flexure is zero. y coordinate 
of centre of flexure is given by 


—3W=/ a. Yu) dedy, 


when € is the boundary of the cross section and R is the region within boundary C. 
® “7: Tr oT 
TE: 
Ze 


x 
T% Bulisagsr Wash 
=) PB T&y)+n| Ta) —zyS@) + ya a 3, 3 T(&, y)}| dx dy 


R 


x 


4 SW) - dx dy 


72 e) 


Ye R > Mi 
=2 | T(x, y) dx dy + y Tl, y) de + I: T(&, y) — xy S(y) + er © dy 
; 6 Ö 


R 


(by Green’s Theorem) . . . 2... (19). 
Again we have the boundary condition 


I, +mr,=0 


1. 
dy|oT W Be oT 
St Se REN ed 
ds | oy sy) 27a | =” = 
or 
S | dy 
$ 
W— eg = K 


The function S(y) may be Fe re 


: 
j 
| 
j 
? 
| 
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We choose S(y) = z x? on C so that the above condition becomes Be = 00n.GC- We 


ds 
can take T=0onC, as ne constant term does not affect the stresses. 
Then from (1.8) we see that stress function satisfies the equation 
Par (ar Ww as | 
TA —2F, TW7 Wu ay’ L=0,8n.G2 211.30) 
= dy 
S(y) = ER on C except that SW) may take any value where ge 0: 
Rn from (1.9) and (1.10) y is simplified and we find that 
ee | Te» ara + ! ayd 1.11 
= wJ U y Te ER er (1.12). 
Now we put 
YLı=z, ET N (1.12), 
then by simple calculation we can show from (1.10) 
@teeoh 2F,W — 1 ds 
os H+ a = — ——i...... 1.15 
9 0 IyM n—/M YM dy, re) 


Let the equation of the boundary C of cross section be transformed to C’ by transformation 
(1.12) and region R within C be transformed to region R’ within C’, then we get following 
boundary conditions 


B=0-0me 
Ww 2 ’ dy, n - 
and S(y,) = 3," C' except that S(y,) may take any value where le 0 (ds, being an 
1 
element of arc of C’). 
By transformation (1.12), we can show that 
1 

= — —— T(z,,y,) dx, dı Hana: de rel): 
Tal] (&, Yı) dx, dy, 3 IyrsMm3 2 Yı dy, (1.14) 


Section bounded be an ellipse and its major axis 


2 2 
Let”, = ae — 1 be equation of the ellipse and’ y = its major axis. We make use of 
transformation e +iy=ccosh (£ +in) 
dl Ccosh cos, Be e= mb, 
ı 
3 
„= 1 &) - = 5 (cosh 2 2E—cos2n). 5 
IF E=a be equation of boundary of ellipse, then a = ccoshx, 
b=csinha; n=0 andn=n give the equation of the major axis. = > 
In this case ° 7 
4b ch sinh a 
RETTET 
As in (1.12) we put s 
& =yLxz=yYLecosh£&cosn, y=YMy=yYMecsinh&sinn. Fig. 1 
Let 
u +iy, = ccosh (& in) 


then 
Er 6, C0Sh £&, 6087, = 6sonh&,sinn, 
Equation of the ellipse is now transformed to another ellipse given by the equation 


yr 


an 


a 
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If & = o, give the boundary of this ellipse then 
YLa=ccohs, yYMb=uasohu und q=L@—MP®. 


In this case also „= 0 and n = give equation of major axis. Now 


1 en) Eu a 91 
ee 1) = -L(coosh2&, —co527) -» - »En.,H4 (2.1). 
From (1.13) we see that the stress function T satisfies the equation 
et @T wi 7, We 2.2 
Ge Sen n/MW+ Tan ee 
a, Mi a 
and 1° 0.000 (1. 8.100 Tatruma N, 
and on eo ken Vader). Sr Regener (2.3). 
(2.2) is satisfied if we assume a particular solution of T 
Anl + R)+ By . (2.4). 
‚-Aylrat us N 
Comparison with (2.2) gives 
2 6 W 2 F, a? W 
rede: map): ae a ea 
> 
Here = en - 


(2.4) gives T,=0 on y =0. But T=0 on &, = a, also. We assume 


T=-Tst4 DA, snh @n KIA sn Au De (2.6). 
n=0 


Since T=0 on .& =, hence we get 
0=AyMbsinn + BM B?sin®n +2 Asnyısinh@ rn +1)u,sin@n-+1)n. 
n=0 


Multiplying both sides by sin (2n + 1)» and integrating between limits (0, z) we get 


8b BM 
2n ) == — = RR er 
a ne u aen—-N@en+N1@n+t3) n=12,3,...) 
and (2.7). 
ee N 2 ] 
AyMi+ Asichy = — a a. 


Thus we get 
cosh? &, cos?n | sinh? & sin? 


Ir — A cz sinh & sin nl I aa M DR 


+ Bei sinh? & sin? 
ie > Anyısinin@n+1)&sin@n+1)n (2.8), 
n=(Ü 


where A, B, Aa. are given by (2.5) and (2.7). In this case the centre of flexure lies on y axis 
and by using (1.14) we can show that y i.e. y coordinate of centre of flexure is given by 


2 - 4AabYLM: „Bab YLM3 


= — 
LMW 15 8 
y (2.9). 
16? BM £ 1 2abyLM 165 
2 / » ee Se) SEE = 
+ Antyün+iye pt a a 7, 
In case of isotropie semi-elliptic section we put 
. o i 
LM, —2F,M= ag. ? u [o = Poisson’s ratio] 


ö- 


nn S 


aaa da au u 
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then 
I y en EIER en, 
x [Boot @n + Da 2] a NIE PENERIGEN 2107: 
where 
er rip Br 3 perl 


In case of isotropic semicircular section of radiusr, we put a=b=randa— + » so that 
cothma— 1 when mis a positive integer, then centre of flexure y of this beam is given by 


te) ee RER N 


In obtaining this value we have noted that 


> 1 BE 
—(@&n—1)(@en+l%@n+s3% 8 128° 


Putting o = 0.3, we find that (2.11) reduces to 
y=.9 117° (appfoximately)” „2 000. 00 RE (Zl2): 
This agrees with the result as given by Timoshenko and Goodier [3]. 
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Visco-elastic steady flow 
By S. K. Sharma 


Durch Abänderung der Spannungs-Dehnungsbeziehungen für zäh-elastische Stoffe werden einige Pro- 
bleme der stationären Strömung in geschlossener Form gelöst. Dabei wird ein außergewöhnlicher Effekt 
gefunden, der weder beim Newtonschen noch beim nicht-Newtonschen Lösungssatz erfaßt wird. Die 
Überlegung läßt 'sich auf das Studium von Gelen anwenden, die von lyophilen Lösungen gebildet werden 
(im Grenzfall z. B. Tafelgelee). Die Ergebnisse erweisen sich als mit den Experimenten als übereinstimmend. 


Using a modified form of the stress-strain relations for visco-elastic materials, some problems on steady 
flow have been solved in a closed form. An extra normal effect lacking in both the Newtonian and the non- 
Newtonian approaches has been found. T'he analysis is shown to have application in the study of gels formed 
by Iyophillic solutions, an extreme illustration of which is table jelly. T'he results are found to be in good 
agreement with experiments. 


En appliquant des relations modifides de tension et de deformation pour des materiaux visqueux et Elastiques 
on a resolu quelques problemes du courant continu d’une facon definitive. On a trouve un effet extra-normal 
qui ne correspond ni aux solutions approximatives Newtoniennes ni non-Newtoniennes. On peut appliquer 
l’analyse aux etudes de gelees se formant aux solutions Iyophylliques dont l’exemple extreme est la gelee de 
table. Il s’avere que les resultats sont bien conformes aux ewpervences. 


IIpnmeHsa BUNON3MeHEeHHbIE OTHOLIEHUA HAIPAsKeHUA NM NePopMalımm 1A BABKOYIACTUYHBIX 
MaTepHaloB, |yHAJIOCh PeIIHTb HEKOTOPBIe 3anaym HempepbIBHOTO TeyeHnst B KOHeYHOM BHJIe. 
Bei 06HapyskeH CBepXHOPMalbAblä IPPeRT, HecoBMecTHMbIa HM C HIOTOHOBCKUMH HN C 
HEeHIOTOHOBCKUMH IOAXOMaMM K peireHMIo. AHaJıma OYeBUJIHO MOSKET ÖbITb HePHMHeH KR 
U3y4YeHNIO 30J1eH, O0PAsyIoINHXCH y JIMOPHIBHBIX PACTBOPpoB, KpaliiuuMm IIPUMePOM KOTOPBIX 
ABJIACTCH OÖBbIUHOE CTOAIOBOE Kelle. ÜOkasbIBAeTcH, YUTO Pe3ybTAaTbI XOPOMO COTJIACYIOTCH C 


OIIBITAMH. 
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1. Introduetion 


The experiments (Weissenberg — 1947) on highly viscous liquids show that the ordinary 
stress-strain relations used in classical hydrodynamics are not sufficient to explain all the col- 
loidal suspensions in the steady laminar flow. For example, in simple cases like shear between 
two parallel plates, some normal stresses are also required in addition to the shear forces. In 
recent years, this problem has been engaging attention of many authors, like Reiner, Rivlin 
(1948), Oldroyd (1950) and Truesdell (1952). The improvement proposed by them over 
classical results consists in, what is called, the second order theory. The second order theory 
approach seems to be complete so long as the fluid does not exhibit appreciable elastic effects. 
A further approach is made by Yoh-Han Pao (1957), who has taken elasticity of fluids into 
account but his approach becomes restricted due to the assumption of zero hydrostatic pressure 
and linear dependence of stress on displacement gradient. _ 

We have taken a more general form of non-linear stress-straih law and attempted to give 
some qualitative results for fluid exhibiting appreciable elastic behaviour. The results obtained 
are more general since those for newtonian and non-newtonian fluids may be deduced as parti- 
cular cases. We have also sketched the velocity profile for some visco-elastic cases which 
affords a comparative study of the visco-elastic and viscous non-elastic phenomena. Shear rate — 
apparent viscosity curve is also drawn in case of an ideal type of starch solution. The results are 
in good agreement with experiments. 


2. Basic equation 


We take the constitutive equation of an incompressible visco-elastic material, also known 
as Maxwellian fluid as follows: 
en rn: 
where $/ is the extra-stress tensor, given by 
G=s—pß, 
Si En Erlen ae u 21.2 Per 2 er u ze 
d; is the rate of deformation tensor, 
t isthe relaxation time constant, 


u is the coefficient of viscosity, 
and ) is cross viscosity-coefficient of the material. 


The equations of motion, in the absence of extraneous forces are: 
ov' i Ri Ü 
0 ag en Km Dawn t  aen . (2.2), 
where o is the density and p is the arbitrary hydrostatic pressure. 
The equation of continuity is 
Ve=ld, N. va me 


3. Reetilinear shearing flow 
The velocity field in this case is 
vl) real ee Rn 
where / is to be determined. 
Ihe rate of deformation tensor assumes the form: 


| Eid 
2d=|f.. Dub] ne 
Kr © 


where dash denotes differentiation with respect to the argument. 
We assume the extra-stress tensor as 


% sy 0 
sel: 5 I rt 
0 d 0 


when we neglect inertia terms. 
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Equations (2.1) and (2.2) reduce to: 


> os} Hi 
stell, 2 A | (3.4), 
OS ; 
et ee ee 0), 
os}, 
Y Y — 1f2 
le Bi (3.6), 
OP 052.054 
a a a een (3.7), 
op Os, Os, 
Een Dil Be a er ee (3.8), 
op 
and a VER ee 8.9): 
All other equations are identically satisfied. 
From equation (3.6), we get: 
SI Al] a a A EE (3.10) 
From (3.5) and (3.10), we get 
ST RE ee ey (3.11) 
Equations (3.4) and (3.11) give: 
A IT ED RAD Se a (3.12). 
Using results (3.10) to (3.12), equations (3.7) and (3.8) reduce to the form: 
op an [20 19% [24 
ae (3.413), 
op € Li 100 
7 ZA WS Se Meer: (3.14). 
Applying condition of integrability to (3.13) and (3.14), we get 
0 „ 1 „ 0 ' 1 
7, el Bauen Ar 0 RER (3.15), 
which gives 
BEREITETE EN BE (3.16), 
where —c is the arbitrary constant pressure gradient along axis of x. 
Equations (3.13) and (3.16) give 
0 
=—e I ee MER (3.17) 
Integration of equation (3.17) gives: 
DD Er en NER (3.18), 


where D is a function to be determined. 
From equation (3.14) and (3.18) we get 


Dy)=2Aff”, 
De N a EL ER I (3.19), 


where k is an absolute constant to be determined. 
From equations (3.18) and (3.19), we have 


jr Eu RL SE Ge EEE (3.20). 


Wr Sry Sr) ER EEE (32T) 


or 


Integrating (3.16), we get 


where FE is an absolute constant. 
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In case of non-newtionian flow, equation (3.21) takes the form | 
BI-=ETTY a (3.22). 


Equation (3.21) can be rewritten in the form 


= =) fs + (&) f=F-y 5) ur Mae ie a 


where 
F = E(ejp), 
or we have the equation 
a a a ee END a (3.24), 


where 7* is a dimensionless parameter and y = (c/u). 


The elastic effect is then governed by 7*, r* = 0 corresponding to the viscous non-elastic 
case, and any non zero value of 7* giving the visco-elastic effect. 


Differentiating equation (3.24) ®- r to f and rearranging we get the equation 


B:*f?+7%f] u EEE NN (3.24.1). 
Integrating equation (3.24.1) with respect to f, we get the equation 
I I 
ei ee ER a = (3.24.2). 


Finally, therefore, we get the equations giving velocity profile in the form 
PP AD= FETT N (3.25), 
3 1 
2 En ala rn re ie. o . (3.26), 
where F and G are constants to be determined from the boundary conditions, and on applying 
the boundary conditions, viz. 


Var, 1 a (3.27), 
we get 


Bes, 1 
M=N(, G= 4 T* Pa + y* Pa De Sch saure ur (3.28), 
where p, is the valueofp= faty= a. 
Equations (3.25) and (3.26) become 
ep pay mee ee DE 


3 Inne. 3 Tan 
WR al ee | *"Btz pi) u ea 


\ 


The velocity profile for all values of r* is nearly parabolic. Consideri >» parti 
an ya r all values of 7* is nearly parabolic. Considering the particular case 


we get 
Pa = —3.J3. 
The equations (3.29) and (3.30) give 
—y= 02) HP 77 De (3.31), 
4 a MS 
= (p* — 118.6) + 9 (PAD 


For viscous non-elastic case, velocity profile is the parabola: 


= —2(f—8), 2 Vie (3.33). 
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ne 


We can now draw out the following table: 


Table A 
tt =. 

p Y f p Yy f 
—3.3 4.0 0.0 1.0 —1.0 6.7 
—3.2 3.8 0.5 1.5 —1.6 6.0 
—3.1 3.7 1.0 2.0 —2.2 5.0 
—3.0 3.5 1.5 2.5 —2.8 3.5 
—2.5 2.8 3.5 3.0 —3.5 1.5 
—2.0 32 5.0 3.1 —3.7 1.0 
—1.5 1.6 6.0 3.2 —3.3 0.5 
—1.0 1.0 ‚6.7 3:3 —4.0 0.0 

0.0 0.0 7 

See Fig. 1. 


Fig.1 


In the case 7* = .04, since p. = —2.9, the equations (3.29) and (3.30) take the form 


Eu > MP IDEE a ea 9), 
1 
— f= .03 [pt — 70.7] + 5 [DERART re se: 
The table will be as below 
Table B 
p Y f p Yy f 
—2.9 4.0 0.0 1.0 —1.0 5.8 
—2.6 3.2 m 1.5 —1.6 5.0 
—2.5 3 2.0 2.0 —2.3 3.8 
—2.0 2.3 3.8 2.5 —3.1 2.0 
—1.5 1.6 5.0 2.6 —3.2 17 
, —1.0 1.0 5.8 2.9 —4.0 0.0 
- 0.0 0.0 6.3 
See Fig. 1. 
E- It is obvious from Tables A and B that elasticity damps the speed of the fluid and this 
2 effect increases with the increase in the value of r*. 
E The stress components are: 
{ wege 2ufen, -| 
i W=cx—k, | 
i Berk if, 6) 
= CU, | 
en) ) 
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Comparing them with the results for non-newtonian fluid flow (Truesdell — 1952), we get the 
extra effects due to elasticity as excesses of our stress components over those in the non-newto- 
nian case. These are respectively 27 f’cy along the normal in the x-direction and the difference 
between the values of A f? from equations (3.21) and (3.22) respectively. 


Speeific resistance Z, is given by 


Beuftirfi.. Sal A 
where f’ is the rate of shear. 
It shows that specific resistance is an odd function of the rate ofshear, which is an experi- 
mental fact. 
From equation (3.39), apparent viscosity app given by (t,/f’) is obtained as 


EB en. 


u} 


where s denotes the shear rate. 

This shows that apparent viscosity increases with shear rate, which has been experimen- 
tally verified by Freudlich and Roder and confirmed by Pryce Jones through his experi- 
ments on certain gels, which exhibit appreciable visco-elastic effects. This has been mentioned 
by Andrade (1951). He has also drawn two curves obtained by Pryce Jones to show a rela- 
tionship between apparent viscosity and shear rate in case the fluid is a Starch solution with 
viscosity at zero shear equal to 30 poises approx. 

Taking 

1 = 30 poises approx., 


and assuming for A 7 the value .1 approx., the equation (3.40) becomes: 
Ya FW LUST rn mn a 


The theoretical curve, therefore, connecting the shear rate and apparent viscosity in case of an 
ideal Starch solution may be drawn as follows: 


90 F 
in / 
Theor / 
60 + % 
S 
S 
.Q 
S 
S 
SQ 
N 
SQ 
N 
ou l L | 
v 70 15 


5 
Rate ofshear (sec, 
Fig. 2 


From Fig. 2 it is clear that there is good agreement between theory and experiment, for 
small rates of shearing. 


4. Poiseuille Flow 
The velocity field in this case may be taken as 
u=V,.2=ß,. Zee Bln), 2 . (4.1) 


the system of coordinates considered being eylindrical. ® is a function to be determined. 
The form of deformation tensor becomes: 


24=|o *)0 <TD 
[0 0 0 


where dash denotes differentiation with respect to the argument. 


S.K. Sharma, Visco-elastie steady flow 319 


We assume the extra-stress tensor as follows: 
SO Se 
si =|0 0 IN ee Een a a IB AR (4.3) 
ser 0 5 
The equations (2.1) and (2.2) on neglecting inertia terms take the form: 


|, 98 ‚ 
&trlo| Ze AN A IE ER (4.4), 
r 05, r G (2 
een ED en E (4.5), 
s+t| as | a r 4.6 
2 3 EEE Braga are Se ae (4.6), 
en 
ar 5 + a2 DR 0 2 TE N ER (4.7), 
op 
no paar 34 A SE Pe a a en (4.8), 
and 
op OSy Eos, 5, 
rad, + ET Eee (4.9) 
From equation (4.4), we get 
SER ER Re EREr (4.10) 
Substituting the value of s; from (4.10) into (4.5), we get 
5 rl B) ae a (4.11). 
From (4.6) and (4.11), we get 
RE TEE A N PR EN JE (4.12). 
Using results (4.10), (4.11) and (4.12), equations (4.7) and (4.9) give 
op a 
ee — Se ar läR.d KASPIUEHIT SENSE BE. (4.13), 
and 
op 77 22 I 1 ’ 3 
le +3/iro o)+— (eo SAT) Bee en, (4.14). 
Applying condition of integrability to (4.13) and (4.14), we get 
Ale ; n I=2| RE 
a s (u.o +TAw)r = 2 ww’ + Er El a a ee: (4.15). 
From equation (4.15), we get: 
d 


Sl 


eo +rrody=—e 
ER SR RE EEE es (4.16), 
D 
uoa+iro:= ne 


where —c is the constant pressure gradient along axis of zand D is a constant to be determined 
from boundary conditions. Equation (4.14) therefore gives 


or on integration " 
a a N A Pa Sr (4.18), 


where D is a function of r, to be determined. 
From equation (4.13) and (4.18), we get: 


Aw’? 


TEIL 
Zu ) 1 1.7 3 2 
= [D(r)] Aw ©" + = 

03% 
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or integrating we get I 
Din =r0t +2 | "ar +4 ET 1 N Ra u 


where A is a constant. 
From equations (4.18) and (4.19), therefore, we get 


2 
p=A—cz+ans4a| © dee ee 
Equation (4.16) takes the form, on applying the condition that ©’ is not infinite when r — 0, 


na +Arot=— gr TEEN 2 


The system of stresses Ü is given as follows: 


‚2 
B=Atez—a|® dr, 
ı2 


2 
B=A+tcz—i dr—iw | 
2) ! 9 
N N dr— Twrc, [ nn. (E22), 
{ I 
E=— zer, 
Belle) 


where w’ is given by the equation (4.21). 

The extra normal effect in this case appears as the excess of our 4, {3 and £ over their 
corresponding values in the second order theory. This effect can be calculated for particular 
values of c from equation (4.21). 


In this case also, as in the previous one the specific resistance is an odd function of the 
rate of shearing, which is experimentally varified. 


The particular case when pressure gradient is taken to be zero corresponds to ® —= 0 as 
can be seen from equation (4.21). In this case stresses = 18 = = a constant, as must be 
the case. 

The shear rate-apparent viscosity curve has its equations as: 


Naon TAT E  R  a (4.23), 


which is the same as in the previous case. Equation (4.21) can be rewritten as 
Er 
ppm au. le 


In the same way as in the previous case the equation of velocity profile will be given by equation 
(4.24) together with 


3 Der 1 
Ä Tr pi + 9 PA=A- 5;® We a EEE 
Applying boundary condition, viz., 


BEN Dur 


equations of velocity profile become: 


1 
mp + pP Meat) nn 
3 * 4 4 1 9% 2 2 1 3 $ 
7% (p!— pi) + 5 (p? — pa) y = Me pe 


a 


The velocity profile in this case can be drawn exactly in the manner it is done in the last case. 


5. Parallel Plate Viseometer 
The velocity field in this case may be taken as follows: 


u=0, v=rol), VE EB ill 


4 
; 
2 
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The deformation tensor takes the form 


I 


0 0 0 
2d=|0 DELTEN 0, ERS A ne et (3.2) 
I 
Extra stress tensor may be assumed as follows 
0 0 0 | 
= 0 RUSS DE Ka ee (33): 
NR 3 
The equations (2.1) assume the form: 
058 | 
sI+r 079 23210 NIEREN sale (54), 
+0 ger ] — z 5.5 
S 30 Ba a 0 a a ln (93% 
5 +:l026| in @ (3.6) 
ä a re 
These equations easily give us 
EU N a N ne (5:7), 
Pa PER 1a a 3 Ka zer a a Me aa mg Po lcd PREIS): 
Se DM Rum Bsaneln Bine fe Re 7 (5.9). 
Equation (2.2) assume the form: 
DE Se 
he I ee (5.10), 
op Ga, 
a a a (5.11), 
ED 0 > 
es (3.12) 
From equations (5.7) to (5.9) and equations (5.10) to (5.12), we get 
e A (5.13), 
Besen BETOE ARE Blut hear A 1 (5.14) 
r00 : 
DB EIN NR RE WERT REN DR Den ER Salt einen as (5.15). 
02 
If we seek out a solution for which ($p/0©) = 0 we get 
EL SR ZEN nl ee Ba Re a A RR re (5.16), 
or 
DE EONSTANE ER WER Volk gie) Enns ha Enge (5.17). 
Equations (5.8) and (5.10) reduce to the form: 
“ NARBEN A nr Ag at ac, I (5.18), 
r 
= EEE ee re N (5.19) 


Equations (5.13) and (5.14) give us 
p=B=,@A+2u Ar ZaMAır RE EEE (5.20): 
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The system of stresses { after applying the condition 5 = 0 at the eylindrical boundary r= a 
take the form: 


HR MR HZ) At, 


B= 104 2U Er Ar + Ar —ay At, 
; a 


1 
E= 2 [BA+uDr—Arn)ajA+ „Art —ay) At, 


= urA+irn As 


If we consider & (= A) to be so small that its cube and higher powers may be neglected, the 
extra normal stresses, given by excess of our 5, 18 and £ over their corresponding values in non- 
newionian fluid are respectively: 


1 
ur A? — a), ur A2Br — dd) and Zar A?(r— a). . . (9.22). 


The velocity profile in this case happens to be the same as in the corresponding viscous non- 
elastic case. 
6. Conelusion 

The present investigation is a contribution to the subjeet of visco-elastie fluids. The 
main conclusions may be stated as follows: 

a) the effect of elasticity in case of steady flow of visco-elastie fluids appears as an extra 
stress in the normal sense, 

b) elastieity of the fluid reduces the speed of flow, 

c) specific resistance is an odd function of the rate of shear, and 

d) the apparent viscosity increases with the rate of shear. 
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Einfache Fehlerschranken beim Eigenwertproblem 
symmetrischer Matrizen 
Von Joachim A. Nitsche 


Das orthonormierte System von Vektoren {y:} approximiere die Eigenvektoren ti einer reellen symme- 
rischen Matrix U. Es werden Abschätzungen für die Genauigkeit von y; und von daraus abgeleiteten Nähe- 
rungen für die Eigenwerte angegeben. Die Schranken sind von einfacher Bauart, zur Berechnung sind im 
wesentlichen nur Matrix- Vektor-Multiplikationen erforderlich. 


Let {hi} be an orthonormal system of vectors approwimating the eigenvectors x; of a real symmetrical matri&‘. 
Estimates are given for the accuracy of the y; as well as of the eigenvalues derived from them. The expressions 
found for the boundaries have a simple structure and may be evaluated by operations consisting essentially 
of matriw-by-vector multiplications only. 2 


On suppose que le systame orthorägl& de vecteurs {bj} soit approxi matif des vecteurs propres x; d’une matrice 
reelle et symmetrique U. Des evaluations pour V’exactitude de v; et des valeurs y derivees pour les valeurs pro- 
pres sont indiquees. Les limites sont simples, pour la caleulation il ne faut avant tout que des multiplications 
malrıce-vecteur, 

Ilycetg CNHCTeMA OPTOHOPMNPOBAHHBIX BERKTOPOB {hy} HAeT TIPHÖNMKeHMEe COÖCTBEHHBIX 
BEKTOPOB Y; JICHCTBHTEJIBHON CHMMETPMYUHOK Marpuusı %. Jlamrest ONeHKU TOYHOCTH y u 
BbIBeEHeHHBIX U3 HMUX DPHÖJNMKEHHBIX COÖCTBEHHLIX 3HAUeHHÜ. CrpyKTypa rpaumm 10BOJIBHO 
IPOCTA; WII BbIUHCJIEHMÄ 110 CYIUHOCTU TPEOYIOTCH TOJLBKO YMHOSKCHNFA MaTpull HA BEKTOPBI, 
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Bei der Lösung von Eigenwertproblemen können i.a. die Eigenvektoren und Eigenwerte 
nur angenähert berechnet werden, so daß eine Abschätzung der Genauigkeit erforderlich wird. 
Dabei sind insbesondere solche Abschätzungen zu bevorzugen, für welche sich die Schranken 
in einfacher Weise berechnen lassen. Mühsame Bestimmung einer Fehlerschranke bedeutet einen 
Aufwand an toter Rechnung, der sinnvoller Weise zur weiteren Verbesserung der Näherungen 


. verwendet werden kann. Bei dem Einsatz von Rechenautomaten kommt weiterhin die Frage 


der Vorbereitungszeit in Betracht. Auch von dieser Seite her ist auf einen möglichst einfachen - 
Aufbau einer Fehlerschranke Wert zu legen. Ur 


Im folgenden sollen Abschätzungen für approximativ bestimmte Eigenwerte und Eigen- 


' vektoren einer symmetrischen Matrix angegeben werden, welche diesen Gesichtspunkten Rech- 


nung tragen. Bei der Berechnung der Schranken sind im wesentlichen nur die Operationen der _ 


Matrix-Algebra erforderlich. Sind diese Operationen für eine Rechenmaschine in Gestalt von 


Unterprogrammen vorrätig, so bedeutet die Programmierung der Fehlerschranke praktisch keinen 
Aufwand. Die gewonnenen Abschätzungen lassen sich in folgender Weise zusammenfassen: 


Bekannt seien n ortho-normierte Vektoreny; (=1,...,n), die näherungsweise mit den 
Eigenvektoren g, der n-reihigen Matrix X übereinstimmen. — Für jedes i wird berechnet 

ii Un yi A d; > & = (yi 2 y; Zr (Yi DJ 9?) 1/2, 

Daraus werden die zwei Zahlen 

II q= Min I — ur) » e = Max ge; 

i+k i 

bestimmt. 

Gilt dann 

V g—2ein.> 0), 
so sind die Eigenwerte A, von X sicher einfach und es bestehen die Abschätzungen 

€ € 
re ee ie ar 
A 172 ie: \R, Y: = _gym i q—eyn 


Ist die Ungleichung V nicht erfüllt, so muß mit dem Auftreten mehrfacher oder benach- 
barter Eigenwerte gerechnet werden. Auch hier ist eine Abschätzung möglich, jedoch ist diese 
nicht in der überaus einfachen Weise wie bei A erreichbar. Auf diesen Fall soll an späterer Stelle 
eingegangen werden. 

Die Ungleichungen A la ssen sich auch sinngemäß anwenden, wenn nur m (< n) Vektoren h; 
bekannt sind, diese aber die Eigenvektoren r; (! —=1,..., m) approximieren, die zu den ersten 
m Eigenwerten /, von X gehören. Diese Eigenwerte müssen jedoch einfach sein. Dann kann I, 
II, V und A übernommen werden, nur laufen jetzt bei I, II, V die Indizes i und k von 1 bis m; 
die Abschätzung A gilt aber nur für i=1,...,m—.1. (Die Numerierung i wird gemäß 
hı ZM > °''" 2 Um angenommen). 


1. Modifikation des Einschließungssatzes von Collatz zur Abschätzung der Eigenwerte 


Wir gehen von einem Einheitsvektor ) aus, der den i-ten Eigenvektor x, von U approxi- 
mieren möge. Den zugehörigen Eigenwert A, approximieren wir durch den mit h gebildeten 
Rayleighschen Quotienten 


ea TE eo I) 
Dann gilt 
AHy—uy=e} be ON) Mr a ne 12), 
wobei 3 orthogonal zu } ist: 
Se FE ARE EEE), 
und e durch 
ESTER ALERT 


festgelegt ist. 
Die orthogonale Matrix T sei so beschaffen, daß die Komponenten des Vektors 


y=Th 


einander gleich werden: 
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Dann geht (2) in R 
| Ad—uj=ej HELSWETT PETE (2a) 


über mit 
und es gilt ebenso wie (3): 
a 1 ET 


Nach dem Einschließungssatz von Collatz!) liegt in dem Intervall 


Min A (d), < A) < Max - (A Hd). 
Ö)x Da 


sicher ein Eigenwert von bzw. von %. Dies sei A,. Nun ist » 


Ad. _ [Or 
(0) TR 


Wegen (5) und (3a) gilt aber |(3).| < 1, so daß wir die Abschätzung 
k—Al<efet .. 7.0 2202 ae . (6) 


=u-+teyn ()%- 


erhalten. 


2. Abschätzung eines zu einem einfachen Eigenwert gehörenden Eigenvektors 


Wie in 1. gehen wir von einem Vektor ) (y’ 9 = 1) und den daraus berechneten Werten u, & 
(vgl. (1) und (4)) aus. Wir wissen, daß in dem Intervall [u — en, «+ e jnj] mindestens ein 
Eigenwert von \, es sei A, enthalten ist. In diesem Abschnitt wollen wir annehmen, daß die 
übrigen Eigenwerte A, (i # k) einen Mindestabstand q von «u haben: 


a a I. SER ERD 
und daß die Ungleichung 
gsi. >. a a 


bestehe. Unter diesen Voraussetzungen können wir unschwer eine Abschätzung für die Dif- 
ferenz y — r; herleiten. Dazu setzen wir an: 


y<-u-rarı (d>0).;. 0... wa a 


Wegen der Normierung von )) und der Vektoren r; gilt 
(1 a + FEN. ER ER FE 
bzw. 
gell (ER, 
woraus die Ungleichung 
DSasS a (a . 0.0. 00r Den 
folgt. Setzen wir die Darstellung (9) für ) in (2) ein, so wird 


kA )A- du +2 u A)at=ez. 
Multiplikation mit %, (j # ü) liefert 
v—A)d=eriz. 
Damit ergibt sich | 


RK (at Diem & > (kr 3)” 
(u—4)* 


und weiter 


& 


Zr <Z. 


=] 


Zusammen mit (11) finden wir so für y — x;| = Y(ai)? + 3 (at)? die Abschätzung 
b-als2lı+2). 


ir aa Einschließungssatz für die charakteristischen Zahlen von Matrizen, Math. Ztschr.48 (1942), 


EL N  \ W 


$ 
\ 
i 
\ 
} 
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3. Beweis von A 
Es seien jetzt n Vektoren y; bekannt, die den Orthogonalitätsrelationen 
9; 9a = Öir (re ler an) 


genügen. Wir bestimmen gemäß I die Größen u, und e; und gemäß II q und e. Nach 1. wissen 
wir, daß in den Intervallen [u — en, u; + e yn] jeweils mindestens ein Eigenwert von X 


‚liegen muß. Gilt nun die Ungleichung V, so sind diese n Intervalle punktfremd, d.h. die Eigen- 


werte von /; sind sicher alle voneinander verschieden und es gilt der erste Teil von A, wie sich 
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unmittelbar aus (6) ergibt. Da nun für ein i die Differenz |u;, —A;| (k # i) bestimmt größer 


als g— e yn ist, 


können wir in den Ausführungen von 2. 


u — A| = q—eyn (Ak), 


Me te 
annehmen, und zwar füry=hY; (i=1,2,...,n). Damit ist auch der zweite Teil von A be- 


wiesen. 


Läuft der Index i nur von 1 bis m (<n), so sind bei Erfülltsein von V m punktfremde 
Intervalle gegeben, in denen nach Voraussetzung die ersten m Eigenwerte A, liegen sollen. Daraus 


folgt weiter 


Um +eyn >% 


bzw. 


u — A| =qg— e yn 


so daß die Ungleichung A für is m—1 gilt. 
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KLEINE MITTEILUNGEN 


Zur Berechnung des Körte-Spiegels 

Hohlspiegellinsen stellen Hohlspiegel aus Glas dar, 
die eine brechende Grenzfläche und eine reflektierende 
Silberfläche besitzen. Sie werden nach ihrem Erfinder 
Mangin-Spiegel oder Körte-Spiegel genannt, wobei 
man sie nach der Gestalt der Grenzflächen und der 
Silberflächen unterscheidet. Ein beliebiger Strahl, der 
von der Lichtquelle ausgeht, wird an der Grenzfläche 
zweimal gebrochen und einmal an der Silberfläche re- 
flektiert. Solche Hohlspiegellinsen werden in der 
Lichttechnik verwendet (in den Kinoprojektoren oder 
in der Signaloptik als Reflektoren). 

Sehr oft werden die Prefilkurven der Flächen im 
meridionalen Schnitt mit trigonometrischen Rech- 
nungen bestimmt. Es gibt auch andere Methoden, 
die das Fermatsche Prinzip in angenäherter Weise, 
teils trigonometrisch benutzen. Hier bedient man sich 
des Fermatschen Prinzips zur Ableitung der For- 
meln, die eine direkte Methode zur Berechnung des 
Lichtweges im Hohlspiegel darstellen. Die Bedingung 
des Stigmatismus auf der Achse genügt zur direkten 
Bestimmung der Gestalt der Grenzfläche bei einer ge- 
gebenen Gleichung der Silberfläche und der Dicke des 
Spiegels. 


Bild1 


Sei 4 ein leuchtender Punkt und B sein Bild (Bild 1). 
Dann liefert das Fermatsche Prinzip folgende Glei- 
chung: 

AC+nDO+nDE+EB=K... (l). 


Die Konstante K bezeichnet dabei die Summe der 
„elementaren Lichtwege‘“‘. Schreibt man nun 


AC+nDC=L, BE+nDE=|, .. 0), 
so lautet die Gleichung (1) 
ae GE 2 Re rr la 


Diese Gleichung stellt eine Ellipse dar, welche ihre 
Brennpunkte in A und B besitzt, und deren große 
Achse 2a =K ist. Zwischen den Punkten dieser 
Hilfsellipse und den Punkten der Profilkurve der 
Spiegellinse kann man Beziehungen aufstellen, die die 
Berechnung der Profilkurven ermöglichen. 

Die Radienvektoren der Hilfsellipse schneiden die 
beiden Profilkurven in Punkten, zwischen denen man 
analytische Beziehungen finden kann. Verlängert 
man die Radien so, daß sie sich im Punkte D’ schnei- 
den, so entstehen zwei Dreiecke ODD’ und EDD’, die 
eine gemeinsame Basis DD’ haben und deren andere 
Seiten in demselben Verhältnis n (n = Brechungsindex 
des Glases) stehen. Also 


ED. ODZECH SED zum): 


Diese Proportion folgt aus den Gleichungen (2), die 
das Prinzip der „Rectification‘‘ der Lichtwege ent- 
halten. Man nimmt weiter an 

SEODED = FI EDEDI ER ERENB): 


So ist die Gerade DD’ eine Normale der Hilfsellipse. 
Daher sind die beiden Dreiecke ODE und CO’ D’E gleich- 
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schenklig, woraus man analytische Formeln für wei- 
tere Rechnungen finden kann. 
Sei der Winkel 
&CDD =y, x{DDC=Bß, 
so gilt 
CD:CD’ = sin f:siny=1:n, 

woraus 
« (6) 
folgt. Der Winkel f ist hier aus den Verhältnissen an 
der Hilfsellipse bekannt. Denn es ist 

x DCD’ = 180° — (P +y)=x m), 


Der Punkt D kann auf der spiegelnden Fläche der 
Linse liegen. Seine Koordinaten lassen sich als die 
Koordinaten des Schnittpunktes der Normale der 


einyw sin B «NM nal 


D' 
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In dem angegebenem Beispiel ist die spiegelnde 
Fläche elliptisch, während die brechende asphärisch 
ist. Die Koordinaten der beiden Flächen wurden nach 
den gegebenen Formeln berechnet. 

Es soll ein Körte-Spiegel berechnet werden, dessen 
hintere (spiegelnde) Fläche eine Ellipse ist (die große 
Halbachse a = 530 mm, Exzentrizität e = 0,679245). 
Das Glas hat den Brechungsindex n = 1,492. Im 
linken Brennpunkte befindet sich der Lichtbogen, im 
rechten das Diapositiv. Die Achsendicke des-Spiegels 
beträgt t = 6 mm. 

Der Koordinatenursprung ist im linken Brenn- 
punkte, die Achse OX in der Richtung der großen 
Achse der Ellipse. Die Koordinaten sind wegen der 
Symmetrie des Systems auf den ersten Quadranten 
bezogen. _ 


ie 
— 
BR a 


Bild 2 


Hilfsellipse DD’ mit der bekannten Profilkurve be- 
stimmen. Also zu jedem Punkte der Hilfsellipse kann 
man immer den entsprechenden Punkt der Profil- 
kurve D finden und daraus die Entfernung dieser 


Punkte DD’ berechnen. Denn 

CD:DD’ = sin P:sin x 
und . (8). 
CD:DD’ = sin ß:sin x 


Aus der Figur (Bild 2) folgt für den Winkel der Ge- 
raden CD mit der Achse AB 


- (9) 


PZN: ! 
wo u, der Winkel des Radiusvektors aus dem Brenn- 
punkte A mit der Achse ist. 

Für die Koordinaten des Punktes C kann man daher 
folgende Gleichungen schreiben, wenn X und Y die 


Koordinaten des Punktes D sind: 
Xc=X—C0D:.co9, Ye=Y—CD:sinp (10). 
In ähnlicher Weise ergibt sich: 


Xg=X+0D:.c089,, Yr= Y’+OD:sing,, 


e (11). 
wo 9% =180°’— u -+X 


So lassen sich die Punkte der asphärischen Fläche be- 
stimmen. Aus der Konstruktion folgt, daß die Be- 
dingung des Stigmatismus auf der Achse erfüllt ist. 

Nach diesem Verfahren zerlegt die Gleichung (2) den 
Strahlengang in zwei Teile, wenn wir annehmen, daß 
wir die Gestalt der spiegelnden Fläche kennen. Denn 
die Beziehung ı = AU’ -+-nCD, bei bekanntem 
Punkte D, läßt die Kenntnis der Lage des Einfalls- 
punktes C noch offen, weil die Bedingung der kon- 
stanten Dicke ohne die Bestimmung der exakten Ge- 
stalt der brechenden Fläche nicht erfüllt werden kann. 
Übrigens ist sie vom formalen Standpunkt vielmehr 
eine Ungleichheit, weil die Dicke an den Rändern des 
Spiegels öfters kleiner ist als an der Achse. 

Die Fertigung der asphärischen Flächen ist schwie- 
rig. In Hohlspiegellinsen trifft man so am häufigsten 
die sphärische Fläche mit einer asphärischen oder eine 
Fläche zweiter Ordnung mit einer anderen asphäri- 
schen. 


Koordinaten der 


Ellipsen-Koordinaten asphärischen Fläche 


— 170,0 0,00 — 164,00 0,00 
—169,61 + 14,84 —163,61 + 14,79 
—168,45 + 29,70 — 162,43 + 29,91 
— 166,30 —+ 44,61 — 160,45 + 44,93 
— 163,74 + 59,60 —157,64 —+ 60,02 
—160,14 —+ 74,68 —154,00 + 75,20 
—155,66 + 89,87 —149,44 + 90,49 
—150,25 —+105,20 —143,97 —+105,92 
— 143,84 —+120,70 —137,46 —+121,50 
— 136,36 —+136,36 — 129,69 —+137,25 
—127,73 +152,22 —121,12 +153,19 
—117,83 168,28 —111,14 +-169,32 
—106,55 —+184,55 — 99,73 -+185,65 
— 93,74 201,02 — 86,80 —+202,18 
— 79,23 +217,68 — 72,16 —+-218,89 
— 62,84 +234,52 — 55,65 +235,76 
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Ein Verminderungsverfahren für den mittleren 
Fehler bei Trägheitsnavigation 


Trägheitsnavigation spielt eine wichtige Rolle in der 
gegenwärtigen Technik. Ihr Grundprinzip besteht 
darin, daß wir unabhängig von irdischen Basen durch 
Messen der Beschleunigung des Fahrzeuges auf seinen 
Ort in einem irdischen Koordinatensystem schließen. 
Bezeichnen wir mit o den Fehler des sich so ergebenden 
Weges, mit da den Beschleunigungsfehler und mit dT 
den Zeitmessungsfehler, so können wir schreiben: 


T+aT 7 
o= f v*t)d— [ vlt) dt, 
0 0 
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wobei 
vlt) = Fat) dt und v*t) = 22 fe + Aue 
ist. a einfacher a bekommen wir: 
— 5 da +2v(T)dT 
und den mittleren Fehler 


2 
7= 7 da+ 2547 


mit 


Diese Resultate sind nicht erfreulich, weil hier die 
Zeit quadratisch auftritt. Unser Ziel ist es, diesen 
Fehler irgendwie zu vermindern. Um dies zu erreichen, 
zerlege ich die totale Fahrzeit 7 in n Teile. Dann gilt 
mitt, +, +. +u=T) 

= 12 a © 2 F 

9, Fe da +2»,dT, o =, da+2nmdT,:, 


, > 
Bar: E 
On= 6; dat 2indaT. 


Der ganze mittlere Fehler im Laufe der Fahrzeit 7 ist 
nach dem Gaußschen Fehlerfortpflanzungsgesetz: 


m V2*< > lo;| < EI) 3420 md. 
i i i 


Das erste Glied auf der rechten Seite erreicht sein Mini- 
mum, wenn, =1,= '-: =1t,. Dann wird mitt = = 
n 


da T? 
= Er — +2nVGmamdT. 
Differenzieren wir den rechtsstehenden 
nach n, so stellen wir fest, daß er bei 


— * 
n DV Dina —n 


seinen kleinsten Wert annimmt. 

Wenn wir also die totale Fahrzeit in n* Teile zer- 
legen, wird die Zeit im obigen Sinne gemessen. Damit 
wird die obere Grenze des mittleren Fehlers 


Ausdruck 


2 x 
= 5 T Yomaz da dT 


V3 


in diesem Falle kleiner als bei der kontinuierlichen 
Zeitmessung, wenn 


I El re 
zo) > (Vomaz + Vömaz — v) 
(bei Imaz = v für alle 7) ist. 
In der Tat, wenn wir im ersten Falle den Zusam- 


menhang von 0, T und v in einem 0, T-Koordinaten- 
system darstellen (% = konst), im zweiten Falle das- 


selbe mit &, T (ümaz = konst) machen, so erhalten wir 
im ersten Falle eine Parabelschar, im zweiten eine 
Geradenschar, die einander in 


N „e2et 
N 2 V3 Vz (Vomaz EZ Vomaz —;) 


und in 


—-2y3 3 (Vömaz + Vomaz — 9) 


schneiden. Wenn also 7 > T', oder % = Ymaz, SO ergibt 
die diskontinuierliche Zeitmessung einen wesentlich 
kleineren Fehler als die kontinuierliche. 
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Zur numerischen harmonischen Analyse 


Wenn eine periodische Funktion durch eine end- 
liche Fourierreihe approximiert wird, so werden bei 
der numerischen Durchführung die Integrale der 
Fourierkoeffizienten durch Summen ersetzt. Die 
Güte der Approximation hängt dabei von der Anzahl 
der benutzten Funktionswerte und dem Verlauf der 
gegebenen Funktion ab. Ist dies eine gerade Funktion, 
was hier (nur der Kürze wegen) zunächst angenom- 


oo 
men wird, so seisiey=g+ 956088. 
s—1 


7 
Wird die Periode 2r in 2n gleiche Teile 7 ge- 
oo 
> r 7 
teilt, so ist 9 = 9% + > Qs COS sg, und man be- 


Sl 
rechnet bei trigonometrischer Interpolation die Fou- 
rierkoeffizienten durch die Summen: 


zn 


l 7T E 
=D yewsro, fürlsr<n); 
e=1 
nn Dummy, Innen 
DE —_ 


Sie werden 


Y=ptgntgntgnt'> 

Ar = gr + gen—r + gen+r + Yan—r Flantr t'*» 

An =mtgIn tGntTnt''- 

Damit ist der Einfluß aufgezeigt, den die Anzahl 2% 
der benutzten Funktionswerte und der durch die Bei- 
werte gs bestimmte Verlauf von y(x) auf die berech- 
neten Fourierkoeffizienten hat.!) Die drei Glei- 
chungen lassen sich durch eine einfache Rechnung 
bestätigen. Dabei ist zu beachten, daß zwar das 


2 
Integral [ cos mx dx — 0 ist für alle von Null ver- 


0 
schiedenen m; nicht aber die entsprechende Summe 
2n 


Sa: 
S "cos mo = Diese bekommt den Wert 2n, wenn m 


PA 


Pe! 
eine der Zahlen 2n; 4n;6n;...ist. 


Setzt man y, in die Gleichung für «a, ein, so er- 
scheinen die Summen 


Sie sind Null für s <r. 

Für s> r kann (s + r) oder (s —r) die Werte 2 n, 
4n,6n,... annehmen, mit denen die Summen in der 
Klammer 2n werden. Also bei s=2n—r; s=4n—-r; 
s=sn—r;,...s=2nHtr;s=4n-+Ftr; br en 
Daraus folgt die Gleichung für ar. Eine gleiche 
Überlegung gibt die Gleichungen für a, und az. 

Bei beliebigen periodischen Funktionen treten in 
y(x) noch Glieder p;  sinsx hinzu. Führt man den 
Gedankengang auch dafür durch, so erhält man die 
Fourierkoeffizienten der Sinusglieder 


br = Pr — Pan—r + Pan+r — Pan—r + Pantr — 


Das für die Abweichung der endlichen Fourierreihe 
von der gegebenen Funktion y(z) maßgebende Integral 


UN= 1 (an + a1 608% + + Mm cosnx + b, sinz 


x: + b„ sinn x — y)’ dx 


2) Die dritte Gleichung findet sich, in einem anderen Zusammen- 
hang, auch bei: D. Jackson. Theory of Approximation. New 
York 1930. 
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läßt sich mit den oben bestimmten Werten von 4r 
und 5, ausrechnen. Dabei wird seine Vergrößerung 
erkennbar, die entsteht, wenn die2n Fourierkoeffizi- 
enten durch Summen statt durch Integrale berechnet 


werden. f , 
Werden sie durch Integrale bestimmt, so ist 


U=naıa tat +FPRıtPfrt')- 


Bei ihrer Berechnung durch Summen vergrößert sich 
U um 


u=2n(gn tan +” +nrlan tn +)? 
n—1 


+aptn & Men—r + gentr + gene Flantr 
+++)? + (Pantr — Pan—r + '' .)?}. 


Bei empirischen Funktionen wird man y = % 
t 
+ 3 (qscossx + ps sin sx) ansetzen, weil der Rest 
s=1 


3 (q5 cos s& + ps sin s x) bei hinreichend großem t 


2; 

unterhalb der Fehlergrenze liegt, mit der y meßbar 
ist. Dann enthalten die Klammern in « die Beiwerte 
dn-+ı bis g und Pn-+ı bis pı. Damit wird ersichtlich, 
wie u mit wachsendem n abnimmt. 


Verfasser: Prof. Dr. H.v.Sanden (Hannover), 
(13a) Behringersmühle, Haus 15 


Unbestimmte Integrale, die Potenzen von Zy- 
linderfunktionen enthalten 


Die Integrale des hier behandelten Typs sind bei 
einer elektronenoptischen Arbeit aufgetreten [1]. Nach- 
dem F. Lenz für ein Integral, das aus zwei solchen 
Termen besteht, in einem speziellen Fall eine ge- 
schlossene Lösung angeben konnte [2], wurde versucht, 
möglichst allgemein gültige Integrationsformeln auf- 
zustellen. Die Beziehungen, die sich dabei ergaben, 
sind in den mir bekannten Integraltabellen nicht ent- 
halten. Sie sollen deshalb trotz ihrer einfachen Her- 
leitung hier wiedergegeben werden. 


Es gilt: 


1 
[arzt 2,10 = T?’ı , 


1 
EZ ZI Zup de = — ] 
Dabei bedeutet Z, die Zylinderfunktion der Ordnung p 
vom Argument x; ! und p sind beliebig. Die Bestäti- 
gung von (1) und (2) ergibt sich unmittelbar durch 
Differentiation der rechten Seiten unter Benutzung der 
bekannten Differentiationsformeln [3] 
d, Py ,\y Pay 
da Ze — Pr Zp + Zp-1 = z% Zp — Zpı (3). 


Im Grenzfall 2 = 0 gehen (1) und (2) über in 


za Zu SED) 


u de=InZ, +ple. : (4), 
Zp 


. (). 


Für!=2 und» =1 ergibt (2) einen Sonderfall einer 
Formel von J. Picht [4]. 
Zwei weitere Formeln sind 


Zp,\! zer zZ rn welP (Zp\ 
a Zp-ı Zprı Zy-ı 2ip Zy-ı 
| a 2 ZT A) 
“ \2,) > 122, % —Z;2,) 


Hierin bedeuten Z, und Z,zwei verschiedene Zylinder- 
funktionen von derselben Ordnung p; der Strich kenn- 


Z 
f Hk=—hi,+phx.. 
Zp 


. (6), 
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zeichnet die Ableitung nach dem Argument x. Auch 
(6) und (7) bestätigt man am einfachsten durch Aus- 
differenzieren der rechten Seiten. Dabei ergibt sich (6) 
wiederum mit Hilfe von (3). Beim Nachmehnen von 
(7) ist zu beachten, daß der Ausdruck & (Z, Z2»—Z; Zp) 
nicht von x abhängt. i 

Setzt man Z, in (7) gleich der Besselfunktion Jy 
und Zy = J_,, so entsteht 

1/Ip\ de _ adIy 
Tee JpJ_, 2lsin (pn) 

Mit Z,—= Jp und Z, gleich der Neumannschen Funk- 
tion N, erhält man 


1 (Re __ # fer 

x \N,) Ip N, 21\N,) ' 
Während (9) für beliebige Werte von p gilt, wird (8) 
für ganze p singulär und identisch mit einem der Lom- 
melschen Integrale [5]. Für beliebige Zylinderfunk- 
tionen erhält man die Lommelschen Integrale aus (7) 
mit/!=0 und /=|]. 

Gleichung (6) liefert für den Grenzfall 1=0 die 

Beziehung 


18). 


. 9). 


2 
[5 «- Int sn 
% Zpr1 Zp-ı 2p Zi 
während für p = 0 aus (6) folgt 
Z3 Zo 
a > 
7 =: In x (11) 
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Konstruktion der Doppelpunkte einer Koppel- 
kurve der Schubkurbel 


Für verschiedene getriebesynthetische Aufgaben 
ist die Kenntnis der Doppelpunkte wichtig, die eine 
durch ein gegebenes Gelenkviereck beschriebene 
Koppelkurve besitzt. Besonders in Fällen, wo die 
Doppelpunkte isoliert oder imaginär, also keinen 
reellen Stellungen des erzeugenden Getriebes zu- 
geordnet sind, interessiert eine konstruktive Bestim- 
mung der Doppelpunkte. 

Die Aufgabe, die Doppelpunkte einer Koppelkurve 
eines allgemeinen Gelenkvierecks zu konstruieren, 
wurde erstmalig von L. Eckhardt [2] gelöst, indem 
ein Kegelschnitt mit dem Grundkreis (Fokalkreis) 
zum Schnitt gebracht wird. Dabei werden die für die 
Konstruktion des Kegelschnittes notwendigen Be- 
stimmungsstücke aus den Abmessungen des gegebenen 
Getriebes errechnet. Diese unbefriedigende Lösung 
ist inzwischen von W. Schmid durch ein rein zeich- 
nerisches Verfahren [6] ersetzt worden, dem der 
Gedanke zugrunde liegt, daß einem solchen durch 
die Doppelpunkte laufenden Kegelschnitt in der 
isogonalen Inversion am Doppelpunktdreieck die 
Verbindungsgerade zweier Punkte der gegebenen 
Koppelkurve entspricht. 


Die Übertragung dieses Verfahrens [6] auf die Auf- 
gabe, die Doppelpunkte einer Koppelkurve der 
Schubkurbel zu konstruieren, liefert — wie noch zu 
zeigen ist — keine angemessene Lösung. Deshalb 
wird im folgenden eine andere Möglichkeit der kon- 
struktiven Ermittlung der Doppelpunkte beschrieben, 


ORTE WERE FETTE WERE 
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wonach diese,dem Grad des Problems entsprechend als 
a von Kreis und Gerade gefunden wer- 
en. 

Zur Erklärung des Übergangs von der allgemeinen 
Koppelkurve zur Koppelkurve der Schubkurbel seien 
in Bild la für eine allgemeine Koppelkurve das Drei- 
eck (D) der drei endlichen Doppelpunkte D,,.,s und 


das Dreieck (F) der drei singulären Brennpunkte 


F1,2,s gegeben. Beide Dreiecke befinden sich in iso- 
gonaler Lage [8]. Läßt man nun unter Beibehaltung 
der von den Geraden F\F,, F,F3 mit der festen 
Strecke D,D, gebildeten Winkel «,& + ß und unter 
Beibehaltung der Entfernung D,F, den Brennpunkt 
F, gegen Unendlich gehen, so erhält man diein Bild 1b 


‚ wiedergegebenen und für die Schubkurbel-Koppel- 


kurve zutreffenden Verhältnisse (beim Übergang vom 
Gelenkviereck zur Schubkurbel spaltet sich bekannt- 
lich die allgemeine Koppelkurve in eine zirkulare 
Kurve 4. Ordnung und die doppelt zu zählende un- 
eigentliche Gerade auf): Das Dreieck (F) der Schub- 
kurbelkoppelkurve besitzt zwei uneigentliche Eck- 
punkte F,o,30 in gegebenen Richtungen gegenüber 
der Geraden d, auf der der dritte Brennpunkt F, mit 
den zwei endlichen Doppelpunkten D,,, liegt [4]. 
Außerdem besitzt die Schubkurbel-Koppelkurve einen 
dritten, uneigentlichen Doppelpunkt D;»; die Rich- 
tung nach D;3% gegen d ix:t durch den Winkel« +2 
gegeben. 

Bild 2 zeigt zur Erläuterung weiterer hier wichtiger 
Gesetzmäßigkeiten als Beispiel eine Schubkurbel 
F,ABF,© mit dem beschreibenden Koppelpunkt P 
in einer beliebigen Stellung. Die von diesem Getriebe 
erzeugte Koppelkurve besitzt zwei reelle Doppel- 
punkte, von denen im Bild einer Knotenpunkt und 
der andere singulärer Doppelpunkt ist. Sie liegen auf 
der Geraden d durch F,, die mit dem Gestell FF,» 
den Koppeldreieckswinkel y einschließt [4]. Die 
gleiche Koppelkurve wird durch eine zweite Schub- 
kurbel F, A, Bı F3» erzeugt [3], Bild2. Durch die 
Fußpunkte der von F, auf die Schubgeraden gefällten 


Bild 2 
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Lote läuft eine Gerade w;, die die Winkelhalbierende 
des Doppelpunktdreiecks in D;» darstellt und des- 
halb die Strecke D,D, in M halbiert und beispiels- 
weise mit F} F,. den Winkel $ bzw. gemäß Bild 1b 
mit d den Winkel « + 2 ß einschließt [7]. Bezüglich 
des Dreiecks (D) entspricht sich in der isogonalen 


Bild 3 


Inversion die Koppelkurve der Schubkurbel selbst, 
und in Analogie zum allgemeinen Gelenkviereck ge- 
hören zu zwei konjugierten Kurvenpunkten P, P* 
zwei bezüglich der Stege F, Fy&, Fı F3© symmetrische 
Stellungen FABF,., Fı Aı Bı F3» der beiden er- 
zeugenden Schubkurbeln F, ABF,«%, F} A, Bı F3%. Die 
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Punkte P, P* entsprechen sich in der isogonalen In- 
version am Dreieck D, D, D;&. Dem Punkt P ist der 


‚Schnittpunkt P* der Geraden P*D,, P*D,, P*D,; zu- 


geordnet, die bezüglich der Winkelhalbierenden 
w;,.,3 spiegelbildlich zu den Geraden PD,, PD,, PD; 
liegen, Bild 3. - 

Nunmehr läßt sich ergänzen, weshalb das Ver- 
fahren [6], bei dem eine Gerade invertiert wird, für 
die vorliegende Aufgabe keine angemessene Lösung 
liefert. Jeder Geraden entspricht in der isogonalen 
Inversion am Dreieck (D) eine Hyperbel durch 
D,,2,3®. Dagegen können und sollen die Doppel- 
punkte D,,, als Schnittpunkte eines Kreises mit der 
Geraden d gefunden werden. Jedem Kreis durch 
D,,, sind aber — wie sich beweisen läßt — in der iso- 
gonalen Inversion bezüglich des Dreiecks (D) sein 
zu d spiegelbildlicher Kreis und die Geraden D, D3o, 
D, D;3o zugeordnet. Deshalb soll eine Konstruktion 
beschrieben werden, die die leicht einzusehende 
Gesetzmäßigkeit verwendet, daß von den Doppel- 
punkten D, , aus die Strecken PD;, P*D} bzw. 
PD}, P*D;3» unter entgegengesetzt gleichen Win- 
keln erscheinen. 


Die dargestellten Beziehungen gestatten nunmehr 
folgende einfache Konstruktion: Gegeben ist eine 
Schubkurbel FRAABF,„ mit P als beschreibendem 
Koppelpunkt in einer beliebigen Stellung, Bild 4. Zu- 
nächst ermittelt man die Geraden d, w; mit Hilfe 
der Winkel y, ß im Koppeldreieck ABP gemäß Bild 2. 
Der konjugierte Punkt P* des beschreibenden Koppel- 
punktes P ergibt sich durch die symmetrische Ge- 
triebestellung. Mit ?, P* und den uneigentlichen 
Punkten D;3%», D5» der Geraden d, ı, ist ein Vierseit 
gegeben, und alle Punkte der Ebene, von denen aus 
die Seiten PD;%, PINS oder P*D,o, PD} je- 
weils unter gleichen Winkeln erscheinen, liegen be- 
kanntlich auf einer gleichseitigen Hyperbel [1]. Die 
Schnittpunkte dieser Hyperbel mit d liefern die ge- 
suchten Doppelpunkte D,,. Um diese Hyperbel- 
schnittpunkte zu konstruieren [1], legt man durch 
P, P* Parallele zu d, die w; in M’, M’” schneiden. Die 
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Kreise um M’, M’ durch P, P* bestimmen die 
Büschelpunkte X, Y, und der Kreis k um den 
Schnittpunkt M von d,w; durch X, Y schneidet auf 
d die gesuchten Doppelpunkte D,,, aus. 


Sind die Büschelpunkte X, Y imaginär, Bild 5, so 
haben die Kreise um M’, M’ durch P, P* und der 
Kreis k um M einen gemeinsamen Orthogonalkreis 0, 
en Mittelpunkt O im Schnittpunkt von PP*, ws 
iegt. / 

Abschließend sei auf eine Schubkurbel F, ABF,% 
eingegangen, die eine Koppelkurve mit konjugiert 
komplexen Doppelpunkten D,,, erzeugt, Bild 6. 
Wenn man die oben beschriebenen Konstruktionen 
für die Ermittlung der Doppelpunkte auf ein solches 
Beispiel anwendet, so zeigt sich, daß die Büschel- 
punkte X, Y unabhängig von der Getriebestellung 
immer imaginär sind und der Orthogonalkreis o 
stets den Mittelpunkt M einschließt. Die konjugiert 
komplexen Doppelpunkte lassen sich als Schnitt- 
punkte der Geraden d und eines Kreises / festlegen, 
der zu o orthogonal und dessen Mittelpunkt auf der 
in M auf d errichteten Senkrechten liegt; mit den 
reellen Elementen d, ! kann die isogonale Inversion 
auch im imaginären Dreieck D, D, D;o konstruktiv 
verfolgt werden [5], [7]. 
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Über die Aquivalenz der drei Stabilitäts- 
kriterien 


Es wird im folgenden ein einfacher Beweis ohne die 
Theorie der komplexen Veränderlichen geführt, daß 
das Hurwitzsche [1], Routhsche [3] und das erste 
Schursche [4] Stabilitätskriterium für Polynome 
mit reellen Koeffizienten gleichwertig sind. 


1. Hilfsbetrachtung. Die Stabilität des Polynoms 
22) = 3 aR-k, a, > 0 mit reellen Koeffizienten wird 
k=0 


nicht geändert, wenn man alle seine geraden bzw. 
alle seine ungeraden Koeffizienten mit derselben 
positiven Zahl multipliziert. 

Beweis. Man multipliziere jeden geraden (un’ 
geraden) Koeffizienten des Polynoms F(x) mit einer 
fest gewählten positiven Zahl b. Das so erhaltene 
neue Polynom sei F,(x). Man bilde die Hurwitzschen 
Determinanten für dieses Polynom F;(x) und dividiere 
dann jede ihrer geraden (ungeraden) Zeilen durch diese 
positive Zahlb. Es entstehen die mit einer Potenz der 
Zahl b multiplizierten Hurwitzschen Determinan- 
ten für das Polynom F(x). Ist b > 0, so ist auch 
diese Potenz immer positiv. Notwendig und hin- 
reichend für die Stabilität des Polynoms F(x) ist 
also die Stabilität des Polynoms F;(x). 


Man kann auch zeigen, daß bei dieser Operation 
die Zahl von Wurzeln mit positivem oder negativem 
Realteil unverändert bleibt. Nach Lienard-Chipart 
[2] ist die Anzahl der Wurzeln mit negativem Real- 


Gy dig — Ay Ag 


ad, Ad, — A,d; ad, dg — Ayde 


teil der Anzahl der Vorzeichenwiederholungen in der 
Folge der Quotienten der Hurwitzdeterminanten: 


u, WB 
Ba 
gleich. Durch Multiplizieren mit einer positiven 
Zahl b wird die Anzahl der Vorzeichenwiederholungen 
nicht geändert und damit auch nicht die Anzahl der 
Wurzeln mit negativem Realteil. Dasselbe gilt von 
der Anzahl der Wurzeln mit positivem eil; sie 
ist der Anzahl der Vorzeichenwechsel derselben Folge 
gleich. 

2. Es werde das Polynom F(x) wie in 1., jedoch 
mit > 0, a, >0 angenommen. Man bilde seine 
Hurwitzsche Determinante H„ des n-ten Grades 


dgy H,, 


1} 
‚dı 4 4; To. Agnı 
| % a on E 
| | Een 
er a == 0 für . 
HH=|0 u 9 ana i>n 
Gy n—4 


. . 


0 0 u 
de 


und man multipliziere jede ihrer ungeraden Zeilen 
mit dem Quotienten a,/a, und subtrahiere sie dann 
von jeder darauffolgenden Zeile. Man multipliziere 
alle geraden Zeilen dieser neuen Determinante mit a, 
und entwickle die Determinante nach dem Laplace- 
schen Satz nach den Elementen der ersten Spalte; 
man erhält so die Determinante H#_, der Ordnung 
n—1: 


1} 
- Ay dgn-a— dp Aan-ı | 


4, (Fe Az a Man 

geh HR ge | 
— ae ee 
Hn = az" Hy-ı 1 0 Ym— pl Ha—mG ... lang —Mlgnz 
0 0 a, EI | 
| 
ern, Ns 
I BE für gerade n, = ar er für ungerade n. 


Es folgt a = H$_,: H„ oder bei der Voraussetzung 
a, >00 
sign H„ = sign Hr» 
Die Hauptunterdeterminante Hm der m-ten Ord- 
nung (1 < m <n) der Hurwitzschen Determinante 
H,„ ist 


a, dg — Ag Ag; Ay Ad — Ay, As; Ay Ad —- u; .» 


| a, Az Ay... Agm-ı ß 
pe BE 
m=|% (A (y . dom |» = i > N. 
Hm = az! VE = az! dı dig 


m—2 .. 

A für gerade m, 
Daraus folgt a77= Hm: H%,_, und bei vorausgesetztem 
a, >o0 


sign Hm = sign HH, 5; 


d. h. die Vorzeichen der Determinanten H, und HA, 
und auch ihrer Hauptunterdeterminanten Hy und 
H#_, sind gleich. 

Die Determinante 4%_, ist die Hurwitzsche De- 
terminante für das Polynom f(x) vom Grade n — 1 
mit reellen Koeffizienten: 


fx) = a, "71 + (a, & — ap a5) a? 
+, 77 + (u — na)... 


Positive Vorzeichen der Hurwitzschen Determi- 
nante und aller ihrer Hauptunterdeterminanten sind 
aber notwendige und hinreichende Bedingung für die 
Stabilität eines Polynoms mit positivem erstem 
Koeffizienten. 


Man multipliziere jede ungerade Zeile dieser Deter- 
minante mit dem Quotienten a,/a, und subtrahiere 
sie von ihrer folgenden Zeile. In der neuen Deter- 
minante multipliziere man jede gerade Zeile durch 
a,; es folgt nach der Entwicklung nach der ersten 
Spalte: 


. Ay Ay m. Gy Ay m-ı | 


a; + Aom-ı 


| 
| 
ji m ei für ungerade m. 

Notwendige und hinreichende Bedingung für die 
Stabilität des Polynoms F(&) mit a, >0, >0 
ist also die Stabilität des Polynoms f(«). 

3. Diese Reduktionsformel zwischen F(&) und f(x) 
folgt auch aus dem ersten Stabilitätskriterium von 
J.Schur für Polynome mit reellen Koeffizienten: 
notwendige und hinreichende Bedingung für die 
Stabilität des Polynoms F(x) vom Grade n mit 
Ay > 0, a, > 0 ist die Stabilität des Polynoms f}(x) 
vom Grade n — 1: 


Fıl®) = ay a, — (Q, Ay — ap 45) © + Ay az 2? 
+1 — 5) Hug -:- 


Man dividiere (siehe 1) jedes ungerade Glied des 
Polynoms f(x) durch a, und setze y= x!; man 
erhält dann /(y) = a-f(x), wo a eine positive reelle 
Zahl ist. 


E 
4 
F 
E 
B 
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4. Man multipliziere im ae x) jedes ‚gerade 
Glied mit ay!; es folgt 


(x) = a Wem (a ZQ0g =) am 


Ay 
+ ag 003 Ha a R)anır... 
1 


das ist die bekannte Form des Routhschen Algorith- 
mus. 


5. Notwendige und hinreichende Bedingung für die 
Stabilität des Polynoms F(x) vom Grade n mit 
reellen Koeffizienten mit a, > 0, a, > 0 ist die Stabili- 
tät des Polynoms f(x) vom Grade n— 1. Dies kann 
mit Hilfe von Hurwitzschen Determinanten be- 
wiesen werden, und das Routhsche und auch das 
Schursche Stabilitätskriterium können in die Form 


dieser Reduktionsformel ee werden. Dar- 
aus folgt, daß diese drei Stabilitätskriterien gleich- 
wertig sind und innig zusammenhängen. 
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BUCHBESPRECHUNGEN 


Dr. 6. Hertz (Prof. a. d. Univ. Leipzig), Lehrbuch 
der Kernphysik. Band I. Experimentelle Ver- 
fahren. 227 S.m.139 Abb. Leipzig 1958. B. G. Teub- 
ner Verlagsgesellschaft. Preis geb. 15,— DM. 


Der erste Band behandelt die experimentellen Ver- 
fahren der Kernphysik, insbesondere die Meßverfahren 
einschließlich der Massenspektroskopie und die Ein- 
richtungen zur Herstellung äußerst energiereicher 
Teilchen. Der zweite Band enthält die Physik der 
Atomkerne im engeren Sinne und bildet somit den 
Hauptteil. Der dritte Band umfaßt die angewandte 
Kernphysik, insbesondere Reaktoren, Isotopentren- 
nung, Radiochemie und Anwendung stabiler und 
radioaktiver Isotope. Dabei liegt der Nachdruck auf 
den wissenschaftlichen Grundlagen. Auf die Tech- 
nologie der Verfahren wird, dem Thema entsprechend, 
nicht näher eingegangen. Das Buch soll als Hilfs- 
mittel für die Ausbildung von Kernphysikern dienen, 
es soll außerdem, wie es im Vorwort heißt, den Inge- 
nieuren, welche sich mit ihren auf anderen Gebieten 


gewonnenen Erfahrungen der Kerntechnik zuwenden, ° 


die Möglichkeit geben, sich eine gründliche Kenntnis 
der physikalischen Grundlagen der neuen Technik zu 
erwerben. 

Der erste Band enthält nach einer allgemeinen 
Einführung über Elektronen, Protonen, Atome und 
Atomkerne die Verfahren zur Bestimmung der Grund- 
größen der Atomkerne, Nachweismethoden energie- 
reicher Strahlung und die Darstellung von Teilchen- 
beschleunigern. Dadurch, daß die einzelnen Abschnit- 
te von den entsprechenden Fachexperten geschrieben 
wurden, ist eine Darstellung von hohem Niveau ent- 
standen. Die beteiligten Autoren sind Bernhard, 
Eckardt, Hartmann, Hertz, Lösche, Schintl- 
meister, Weiß. Der Vergleich verschiedener Ab- 
schnitte verrät eine hohe Geschlossenheit der Dar- 
stellung, um die sich der Herausgeber offenbar viel 
Mühe gemacht hat. Damitistein Buch über die Grund- 
lagen der gegenwärtigen Kernphysik entstanden, dem 
man im Interesse dieses Fachgebiets nur eine weite 
Verbreitung wünschen kann. 


Dresden W. Macke 


Dr. F. Regler (o. ö. Prof. a. d. TH Wien), Atont- 
physik. Grundlagen und Anwendungsbeispiele für 
Metallkundler und Ingenieure. Herausgeber: Ver- 
einigte Aluminiumwerke AG — Österreichische 
Metallwerke AG, Ranshofen, Sonderheft 4. 43 S. m. 
39 Abb. u. 7 Tabellen. Salzburg-Großgmain 1956. 
Technische und Fachverlagshandlung F. Wessiak. 
Preis geh. 42,— ö.S. 


Die kleine Schrift bringt in geschickter, leicht- 
verständlicher Darstellung eine gedrängte, zuver- 


lässige Übersicht über diejenigen Gebiete der Atom- 
physik, für die in weiten Kreisen Interesse besteht, 
wobei verständlicherweise die friedliche Anwendung 
der Atomenergie besondere Berücksichtigung fand. 


Dresden Gustav E.R. Schulze 


W. Weizel (o. Prof. d. Physik a. d. Universität 
Bonn), Lehrbuch der theoretischen Physik. 
Zweiter Band: Struktur der Materie. 2. Aufl. XV + 
989 S. m. 268 Abb. a 1958. 
Springer-Verlag. Preis geb. 88,— DM 


Nunmehr liegt auch der zweite Band des bekannten 
Lehrbuches der theoretischen Physik von W. Weizel 
in zweiter Auflage vor, gegenüber der ersten Auflage 
um etwa !/; an Umfang gewachsen. Da der Charakter 
des Werkes unverändert blieb, seien nur die wichtig- 
sten Ergänzungen erwähnt. In einem rund 100 Seiten 
umfassenden Abschnitt ‚‚Feldtheorie der Materie‘ 
werden die Gedankengänge der Theorien der Elemen- 
tarteilchen und ihrer Wechselwirkungen dargestellt. 
Der Abschnitt ‚„‚Kernphysik“ ist etwa auf das Doppelte 
angewachsen. Um das Rüstzeug für diese Abschnitte 
bereitzustellen, wurde der Abschnitt ‚‚Quantentheorie‘“ 
um ein Kapitel ‚‚Stoß- und Streuprozesse‘“ erweitert. 
Ferner entstand ein neuer Abschnitt ‚‚Elektronik‘‘, 
in dem die einschlägigen Kapitel der vorigen Auflage 
mit zwei neu hinzugekommenen (,Relativistische 
Elektronen- und lIonenoptik. Teilchenbeschleuni- 
gung‘ und „Emission, Neutralisation und Absorption 
von Ladungsträgern an Oberflächen‘) vereinigt 
wurden. Im Abschnitt ‚Struktur und Eigenschaften 
der zusammenhängenden Materie‘‘ wurde ein kurzes 
Kapitel über Halbleiter eingefügt. 

Auch die im wesentlichen unveränderten Ab- 
schnitte hat der Verfasser in Einzelheiten über- 
arbeitet und ergänzt, sowohl aus sachlichen als auch 
aus didaktischen Erwägungen. Sein Bemühen um 
ständige Verbesserung der Darstellung zeigt sich z. B. 
darin, daß teilweise auch Paragraphen über bereits 
klassisch gewordene Gegenstände umgeschrieben wur- 
den. So wird jetzt nach dem Hinweis auf das zwangs- 
läufige Versagen der klassischen Physik am Atom- 
modell die Wellenmechanik erst nach dem Besprechen 
der Beugung von Materiestrahlen eingeführt. Oder 
bei den Spektren der Alkalien ist der unveränderte 
Text durch eine sehr schöne Wiedergabe eines 
Natrium-Absorptionsspektrums ergänzt worden. 

Der ‚‚Weizel‘“ bedarf keiner Empfehlung mehr. 
Es ist ihm aber zu wünschen, daß sich der hohe 
Preis — die Ausstattung ist freilich vorzüglich — 
seiner Verbreitung nicht als ernstes Hindernis ent- 
gegenstellt. 


Gustav E. R. Schulze 
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Nicht das Beweisen von anscheinend Selbstver- 


En 


R. Stender, Heft 5). 238. m. 28 Abb. Frankfurt am 
Main/Hamburg 1956. Otto Salle Verlag. Preis geh. 
2,80 DM. 


ständlichem, sondern das Erkennen, daß das an- 
scheinend Selbstverständliche nicht selbstverständ- 
lich ist, das ist die Aufgabe des Unterrichtes. Von 
diesem Grundsatz ausgehend plädiert der Autor dafür, 
die Parallelenlehre, dieses entscheidende Fundament, 
auf dem die euklidische Geometrie basiert, bereits 
auf der Oberstufe, in gewissem Umfange sogar schon 
auf der Mittelstufe zu besprechen. i 

In klarer, didaktisch richtiger Weise wird sodann 
das Parallelenproblem in Euklidscher und Hilbert- 
scher Fassung gebracht. Sodann folgen verschiedene 
andere Arten der Parallelendefinition bis zum heutigen 
Aufbau der Parallelenlehre. Betrachtungen über 
Fernelemente beschließen dieses lesenswerte Büchlein. 


Dresden R. Bereis 


Dr. K. Fladt, Elementarmathematik vom 
höheren Standpunkt aus. Teill: Elementare 
Arithmetik. 80 8. Stuttgart 1957. Ernst Klett Ver- 
lag. Preis kart. 6,20 DM. 


Vorliegendes Büchlein über ‚Elementare Arith- 
metik“ ist der erste Teil der ‚„Elementarmathematik 
von höherem Standpunkt aus‘, eines Werkes, mit 
dem der Autor den Leser mit all den für die Schule 
wichtigen Grundlagen der Algebra und Elementar- 
geometrie bekanntmachen will. Der Autor wendet 
sich daher insbesondere an die Mathematiklehrer an 
den höheren Schulen, sowie an jene Studenten, welche 
Lehrer der Mathematik an höheren Schulen zu werden 
gedenken. Es wird demnach vor allem die Geschlossen- 
heit und innere Folgerichtigkeit des neuen Aufbaues 
der Schulmathematik dargetan und durch Aufnahme 
von modernen mathematischen Begriffen versucht, 
die Kluft zwischen Hochschul- und Schulmathematik 
etwas zu verringern. 

Nach einer kurzen Einführung über den Inhalt 
und die Methode der Elementararithmetik wird der 
Leser in die Axiomatik der elementaren Arithmetik 
eingeführt und sodann unter Verwendung der Begriffe 
von Körper, Ring, Isomorphismus usw. mit der 
Entstehung der rationalen und komplexen Zahlen 
vertraut gemacht. Nach Aufstellung eines Axiomen- 
systems der Arithmetik und Scheidung der einzelnen 
Zahlenarten werden rückwärtsgehend Fragen nach 
dem Zustandekommen der natürlichen Zahlen auf- 
geworfen und in diesem Zusammenhang die fünf 
Peanoschen Axiome besprochen. Die exakte und 
klare Darstellungsart macht die Lektüre dieses 
bemerkenswerten Werkes jedem Freunde der Mathe- 
matik zu einem reinen Vergnügen. 


Dresden R. Bereis 

Dr. K. Fladt, Elementarmathematik vom 
höheren Standpunkt aus. Teil2: Elementar- 
geometrie l. 1548. m. Abb. Stuttgart 1957. Ernst 
Klett Verlag. Preis kart. 13,20 DM. 

Vorliegendes zweites Heft der ‚‚Elementarmathe- 
matik von höherem Standpunkt aus“ bringt die 
elementare Geometrie der Ebene und des Raumes. 
Ausführliche Darlegungen der Axiomensysteme Eu- 
klids und Hilberts, sowie deren Beziehungen zur 
Kleinschen Auffassung der Geometrie als Invarian- 
tentheorie geometrischer Gruppen im Sinne des 
Erlanger Programms bilden den Kern der Betrach- 
tungen. Auf diese Weise wird die elementare Geo- 
metrie zur Geometrie der Hauptgruppe (Bewegungen, 
Umlegungen und Dehnungen). Mit voller Klarheit 
wird die Umwandlung der klassischen Geometrie in 
die moderne Verwandtschaftsgeometrie vollzogen, 


Dr. K.Fladt, Ein Kapitel Axiomatik: Die wobe 
Parallelenlehre. (Schriftenreihe zur Mathematik, le 
herausgeg. v. Prof. Dr. Fr. Drenckhahn u. Dr. 
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cht mehr « 
Vorliegende übersichtliche 


Dr. V. J. Rossow, An Analysis of the Error In- 


volvedin Unrolling the Flow Field in Turbine 


Problems. (Mitteilyngen aus dem Institut für Aero- 


dynamik a. d. ETH in Zürich, Nr. 23.) 528.m. 18 Abb. 
Zürich 1958. Verlag Leemann. Preis brosch. 16,— SFr. 


Die Schaufeln axialer Strömungsmaschinen werden 


üblicherweise so berechnet, daß eine gedachte koaxi- 
ale Schnittebene abgewickelt und die Schaufelumströ- 
mung zum ebenem Problem vereinfacht wird. Faktisch 


handelt es sich in der Strömungsmaschine aber um ein 


dreidimensionales Strömungsproblem. Schon früher 
wurde im gleichen Institut von Ginsburg eine ähnliche 
Aufgabenstellung bearbeitet. (Referat im Band 37 
(1957), S. 315 dieser Zeitschrift). Auch hier sollte der 
Fehler bestimmt werden, der durch Abwickeln der 
Schaufeln in die Ebene entsteht. Während in der frühe- 
ren Arbeit die Störgrößen in unmittelbarer Nähe derWir- 
bel selbst berechnet wurden, wird in der vorliegenden 
Schrift auch die Größe der Störfunktion in einigem 
Abstand von den Wirbeln festgestellt. Zusätzlich wird 
der Einfluß von Quellfäden bestimmt. Die Voraus- 
setzung der Reibungsfreiheit wird auch in dieser Ar- 
beit gemacht unter Hinweis darauf, daß nach der 
Prandtlischen Grenzschichttheorie die Reibung addi- 
tiv zur Potentialströmung ohne nennenswerte Störung 
hinzugefügt werden darf. Das ist richtig für Gebiete 
einer Schaufel, die weit (verglichen mit der Sehnen- 
länge) von den begrenzenden Wänden entfernt sind. 
Nähert sich aber das Nabenverhältnis 1, so wird bei 
ausgeführten Maschinen das Verhältnis Wandabstand 
zur Sehnenlänge klein und die Störwirkung der Rei- 
bung auf die Grundströmung ist erheblich. Das ist zu 
beachten, weil bei Vernachlässigung der Reibung beste 
Übereinstimmung zwischen ebener und räumlicher An- 
ordnung der Schaufeln gerade bei Nabenverhältnis 
= 1 festgestellt wird. 

Die Ableitungen der Arbeit sind soweit vorangetrie- 
ben, daß nun numerische Rechnungen durchgeführt 
werden könnten. Leider vermißt man in dieser Schrift 
Beispielsrechnungen, mit denen eine Abgrenzung von 
Gebieten vorgenommen wurde, in denen die Fehler 
durch das Abwickeln der Schaufeln in die Ebene einen 
bestimmten Betrag nicht überschreiten. 


Dresden W. Albring 


P. Patin (Ingenieur en Chef ä la 8. N.C, F.), Les 
Transmissions de Puissance et la Variation 
de la vitesse. Mit einer Vorrede von A. Caquot. 
338 8. m. 318 Abb. u. 8. Tafeln. Paris 1956. Editions 
Eyrolles. 

Das Buch befaßt sich mit dem Problem der Anpas- 

sung von verschiedenen Antriebsmaschinen durch 
Zwischenschaltung von Getrieben oder Kupplungen 
an Arbeitsmaschinen beliebiger Charakteristik, wobei 
der Steuerung und Regelung besondere Aufmerksam- 
keit geschenkt wird. 
‚ Im ersten Teil werden nach einer kurzen Zusammen- 
fassung der Festigkeitsberechnung die Maschinenele- 
mente, die bei der Kraftübertragung verwendet wer- 
den, z.B. Kupplungen, Zahnräder, Lager, und die wich- 
tigsten Bauteile für die hydraulische und elektronische 
Regelung, z. B. Elektroventile, Trioden und Ampli- 
dyne besprochen. 
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Im zweiten Teil sind die Getriebe bzw. Getriebezüge 
zur Kraftübertragung und ihre Regelung zusammen- 
gestellt. Stufenlos regelbare Reibradgetriebe, viel- 
stufige Zahnrad-Wechselgetriebe, hydrostatische und 
hydraulische Getriebe sind beschrieben und ihre wich- 
tigsten Eigenschaften herausgestellt. Die Probleme 
des elektrischen Antriebes mit Gleich- und Wechsel- 
strom für mechanische und elektronische Regelung 
runden das Bild ab. 

Der Verfasser hat es verstanden, das Wesentliche 
der mechanischen, hydraulischen und elektrischen 
Leistungsübertragung leicht verständlich und kurz 
gefaßt zu bringen und an bewährten Beispielen aus der 
Praxis die Möglichkeiten der Anpassung, Regelung 
und Steuerung zu demonstrieren. 


Dresden St. Fronius 


K. Menninger, Zahlwort und Ziffer, eine Kul- 
turgeschichte der Zahl. Band 2, Zahlschrift und 
Rechnen. XI + 314S. m. 255 Abb. 2. neubearb. u. 
erweiterte Aufl. Göttingen 1958. Vandenhoeck & 
Ruprecht. Preis geb. 24,80 DM. 


Während der 1. Band das Zahlwort und die Zahl- 
sprache behandelt, ist der vorliegende 2. Band der 
Zahlschrift und dem Rechnen gewidmet. Er vermit- 
telt einen reizvollen Einblick in die spannenden Vor- 
gänge, die sich beim Einbruch der indischen Stellen- 
schrift in die abendländische Welt abspielten. Auch 
dem Vertreter der angewandten Mathematik werden 
viele Anregungen geboten. 

Der Verfasser hat’ein gewaltiges Material verarbei- 
tet. Trotzdem verliert der Leser niemals den Über- 
blick über diese Vielfalt, weil jeder Abschnitt mit einer 
klaren Gesamtorientierung beginnt. Was äußerlich das 
Satzbild zur Erleichterung der Übersicht beitragen 
kann, ist wohlabgewogen berücksichtigt. 

Über die Fingerzahlen und die Zahlzeichen des Vol- 
kes (Kerbhölzer, Knoten) führt der Weg zu den Buch- 
stabenziffern, den römischen Zahlzeichen und dem 
Rechenbrett. Sodann wird die Geschichte unserer 
Stellenschrift behandelt. Auf dem Hintergrund un- 
serer abendländischen Entwicklung wird im letzten 
Abschnitt ‚„‚Zahlwort und Ziffer in China und Japan“ 
dargestellt. Es ergibt sich, daß dort Sprache, Schrift 
und Zahlschrift aus ein und demselben Mutterboden 
gewachsen sind, während für uns — wie es der Ver- 
fasser zum Schluß knapp formuliert — gilt: ‚Wir spre- 
chen deutsch, wir schreiben römisch, wir rechnen in- 
disch‘“. 

Viele z. T. wertvolle Abbildungen beleben das Stu- 
dium. Ein reichhaltiges Literaturverzeichnis und ein 
recht gut orientierender Sachweiser sind dem Buche 
beigegeben. 

Potsdam M. Draeger 

Dr. H.R. Voellmy, Experimentelle Untersu- 
suchungen an verschieden stark konvergen- 
ten, schlanken Rotationskörpern bei mäßig 
hohen Überschallgeschwindigkeiten. (Mittei- 
lungen aus dem Institut für Aerodynamik a. d. ETH 
Zürich, Nr. 24.) 828. m. 37 Abb. Zürich 1958. Ver- 
lag Leemann. Preis geh. 17,60 Fr. 

Zur Nachrechnung der aerodynamischen Kräfte an 
schlanken Drehkörpern sind schon früher verschie- 
dene Rechenmethoden entwickelt worden. Um die 
Schwierigkeiten der Integration einer nichtlinearen 
Potentialgleichung für gasdynamische Strömungen zu 
vermeiden, wurden bei einer Gruppe von Lösungsver- 
fahren die Gleichungen linearisiert. Die Körperachse 
wird mit Singularitäten belegt, deren örtliche Ergie- 
bigkeit so bemessen ist, daß die Bedingung tangenti- 
aler Strömung an der Körperoberfläche erfüllt wird. 
Aus den damit berechneten Geschwindigkeiten auf 
dem Körper sind auch Drucke und Kräfte zu bestim- 
men. Im Sinne einer Theorie, die die Reibungswir- 
kung vernachlässigt, würde man genauer arbeiten, 
wenn man mit dem Lösungsverfahren unmittelbar an 


die nichtvereinfachte nichtlineare Potentialgleichung 
anschließt. Nur für den Kegel ist durch Busemann eine 
strenge Lösung bekanntgeworden, für nichtkegelige 
Drehkörper wurden Charakteristikenverfahren aus- 
gearbeitet (Oswatitsch, Sauer, Tollmien und Frankl). 

Es ist nun recht verdienstvoll, an einer Systematik 
von schlanken Drehkörpern mit Parabelkonturen die 
Aussagen der verschiedenen Rechenverfahren experi- 
mentell nachzuprüfen. Die Reynoldszahl — bezogen 
auf die Länge — betrug 1,2. 10%, die Machzahl der 
Anströmung wurde im Bereich zwischen 0,28 bis 1,67 
und der Anstellwinkel zwischen 0° und 12° verändert. 

Die Vergleiche mit der Theorie für reibungsfreie 
Strömung haben zum Teil sehr große Abweichungen 
erkennen lassen. Dabei hat sich gezeigt, daß eine Ver- 
feinerung der linearisierten Methode oder die Benut- 
zung einer Lösung der nichtlinearisierten Potential- 
gleichung keine grundsätzliche Verbesserung erwarten 
läßt. Entscheidend ist die Vernachlässigung der Rei- 
bung. Während bei der Bestimmung der Axialkraft in 
rotationssymmetrischer Strömung eine additive Be- 
rücksichtigung des Reibungswiderstandes statthaft 
ist, liegen die Verhältnisse bei der Normalkraft kom- 
plizierter. Hier spielen Wirbelbildung und Ablösungs- 
erscheinungen an der Leeseite des angestellten Kör- 
pers eine entscheidende Rolle. Der Berichter hat selbst 
schon im Jahre 1942 Windkanalmessungen durchge- 
führt und festgestellt, daß im inkompressiblen Strö- 
mungsbereich für Anstellwinkel zwischen « = 0° und 
& = 90° die Darstellung der Meßwerte in der Form 
C„/sinl® & über der Reynoldszahl (c» - d/v - 1/sin «) 
den besonders bei großen Anstellwinkeln erheblichen 
Einfluß der Reynoldszahl wieder eliminierte. In den 
genannten Koordinaten stellt sich der Normalkraft- 
beiwert ähnlich dar wie die Widerstandsbeiwerte eines 
Kreiszylinders. Deutlich ist ein Unterschied zwischen 
dem laminaren und turbulenten Bereich zu erkennen, 
der Normalkraftbeiwert ändert sich dabei fast um den 
um den Faktor 3! Bei den vorliegenden Messungen 
tritt wegen der noch relativ kleinen Anstellwinkel die- 
ser Effekt zurück, er könnte aber Anschlußwerte zu 
größeren Reynoldszahlen und größeren Anstellwin- 
keln liefern. 


Dresden W. Albring 


Dr.-Ing. W. Kaufmann (em. o. Prof. d. Mechanik a. 
d. TH München), Technische Hydro- und Aero- 
mechanik. Zweite verb. u. ergänzte Aufl. VIII + 
386 S. m. 266 Abb. Berlin/Göttingen/Heidelberg 1958. 
Springer-Verlag. Preis geb. 37,50 DM. 


Man schreibt Prandtl die Äußerung zu, ein Buch 
über die Aeromechanik müßte eigentlich in drei Fas- 
sungen geschrieben werden: für Physiker, für Mathe- 
matiker und für Ingenieure. Beim Lesen des vorlie- 
genden Werkes von Kaufmann hat man den Ein- 
druck, daß in einem gewissen Maße den Belangen die- 
ser drei Sparten gemeinsam Rechnung getragen ist. 
Neben der physikalischen Analyse ist der mathema- 
tischen Beschreibung ein breites Feld eingeräumt, und 
durch manches Beispiel angeregt, werden auch die 
Ingenieure den Weg zur technischen Anwendung fin- 
den. Der Verfasser hat in diesem Buch, das schon vor 
20 Jahren einen Vorläufer in der zweibändigen ‚‚Ange- 
wandten Hydrodynamik“ hatte, den Inhalt und sicher 
auch die pädagogischen Erfahrungen seiner Vorlesun- 
gen über dieses Gebiet verarbeitet. Gegenüber der 
1955 auf S. 198 dieser Zeitschrift besprochenen 1. Auf- 
lage der Neufassung sind nur kleine Änderungen und 
Erweiterungen vorgenommen worden. Bemerkenswert 
ist aber, daß der Verfasser den Empfehlungen über die 
Einführung eines neuen Maßsystemes nur so weit ent- 
gegenkommt, daß er die Bezeichnung [kg], die bisher 
für Masse und Kraft nebeneinander bestand, nun mehr 
für die Kraft nicht mehr benutzt und Kilopond [Kp] 
schreibt, im übrigen aber bei dem wohlbewährten 
technischen Maßsystem bleibt. 
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Jahrbuch 1957 der Wissenschaftlichen Ge- 


sellschaft für Luftfahrt e. V. (WGL). Mit den 
Vorträgen der WGL-DVL-Tagung in Essen vom 9. 
bis 12. April 1957. Herausgegeben von H. Blenk. 
VII + 504 S. m. 623 Abb. u. 23 Tab. Braunschweig 
1958. Friedr. Vieweg & Sohn. Preis geb. 58,— DM. 


Außer gesellschaftsinternen Mitteilungen und einer 
Reihe von Nachrufen enthält das neue Jahrbuch 1957 
der Wissenschaftlichen Gesellschaft für Luftfahrt 
(WGL), das jetzt in der nun schon gewohnten hervor- 
ragenden Ausstattung erschienen ist, den Wortlaut 
des größten Teils der auf der Essener Jahrestagung im 
April 1957 gehaltenen Vorträge. Jedem Vortrag sind 
begrüßenswerterweise Auszüge aus der im Anschluß 
an ihn gehaltenen Diskussion und ausführliche Schrift- 
tumshinweise beigefügt. Es würde den Rahmen dieser 
Besprechung schon sprengen, wollte man auch nur die 
Titel der 49 Vorträge anführen. Es bleibt daher wieder 
nichts anderes übrig, als einige Teilgebiete der viel- 
seitigen und weitgespannten Luftfahrtforschung her- 
auszugreifen und auf die zugehörigen Vorträge hinzu- 
weisen. 

Auf dem Gebiete der Aerodynamik finden sich Be- 
richte über die Einrichtung des neuen Aachener Über- 
druckwindkanals (Naumann. Heyser, Trommsdorff), 
über eine Methode zur vereinfachten Bestimmung der 
Druckverteilung an Flügeln beliebiger Spannweite 
(Oswatitsch, Laschek), über die Strömung an Körpern 
nicht mehr kleiner Streckung in auftriebsloser Unter- 
und Überschallströmung (Keune), über den Auftrieb 
und Widerstand langer Ringflügel in Überschallströ- 
mung (Zierep) und über die Berechnung des induzier- 
ten Geschwindigkeitsfeldes von Tragflügeln (Gersten). 
Mit Untersuchungen zur Aerodynamik von Flügel- 
Rumpf-Anordnungen beschäftigen sich X. Hafer und 
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Verwendung von Freifluganlagen an Stelle von Wind- 


kanälen sind erwähnenswert (Bull, Kutterer), Daß 
darüber hinaus noch andere Gebiete der Luftfahrt- 
forschung zu ihrem Recht gekommen sind, geht aus 
den Vorträgen über Meteorologie, Navigation, biolo- 
gische Strahlenwirkung usw. hervor. 

Alles in allem gibt das neue Jahrbuch wiederum 
einen überaus aufschlußreichen Querschnitt durch das 
gesamte Gebiet der Luftfahrtforschung. 

Dresden.  H. Heinrich 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 
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GAMM-Mitteilung 


Der Vorstand der Gesellschaft für Angewandte 
Mathematik und Mechanik (GAMM) hat auf seiner 
Sitzung am 19. 5. 1959 in Hannover einen Fachaus- 
schuß für Unternehmensforschung (opera- 
tions research) gegründet und die Leitung Herrn 
Prof. Dr. H. Görtler, Freiburg i. B., Sonnhalde 90, 
übertragen. An der Arbeit dieses Fachausschusses 
Interessierte werden gebeten, sich mit Herrn Prof. 
Görtler in Verbindung zu setzen. 
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